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1 Introduction 


The occurrence of vortices is a widespread phenomenon in physics. They are well 
known, for instance, from fluid mechanics of gases and liquids [Bes06}|ZB22], but 
vortices also occur at cryogenic temperatures in superfluids and in the condensate 


of superconducting charge carriers in solids [EH05; |Tin04]. In the latter case, 
vortices of supercurrents exist in Josephson junctions that are superconducting 


leads separated by a tunnel barrier |Jos62;/AR63]. Here, if the junction is tailored 
to an appropriate shape, it allows for spatial and temporal phase distribution with 
changes in integer multiples of 27 [Jos64; |Swi61]. This phenomenon gives rise 


due to a Bohr-Sommerfeld-like quantization of the superconductor’s phase and is 
analogous to the angular momentum quantization in atoms [Lon50;/DN61}|DF61]. 
As a characteristic of superconductors, such phase drops drive supercurrents, which 
are closed rings, namely Josephson vortices. They enclose a quantized magnetic 


flux of Po = h/2e [BL64}|Jos65], often referred to as “fluxon” [BP82} |Lik86]. 


The basic understanding of the phase dynamics of these vortices comes from the 
sine-Gordon model [Kul65} {Sco69}/FD73], which features a nonlinear wave equation 
for the Josephson medium. Its vortex solution is of particular interest [Ust98], as 
these excitations are pinned at the tunnel barrier and may propagate along the 
nonlinear Josephson medium as a wave packet that maintains its shape, i.e., as 
a soliton [Swi61]. Fluxons in long Josephson junctions have been a rich research 
field for decades since they serve as an excellent physical system to study highly 
nonlinear electrodynamics and soliton dynamics on a microscopic scale 
Ped83]. Thereby, classical effects such as high-frequency cavity resonances 


FD73; |FD74b Kul65} |Kul67], but also relativistic mass increase and Lorentz 
contraction |MS78} Lau+95] have been studied. Among others, applications for 
microwave generation |Nag+83 and amplification [Nag+85 based 


on flux-flow oscillators have been developed. 


Although Josephson vortices host quanta of magnetic flux, only their phase is 
quantized. However, they do not behave like quantum particles because both 
phase and charge of the vortex commute and hence do not feature Heisenberg’s 
uncertainty relation. 


The key question is whether and how the Josephson medium can be modified so 
that Josephson vortices behave like quantum particles. 
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Experimentally realized Josephson quantum particles are interesting, as there are 
theoretical proposals to explore relativistic quantum sine-Gordon physics [Roy+21], 
and quantum coherent excitations in spatially periodic potentials simulating solid- 
state models. Moreover, topological solitons would theoretically be auspicious 
candidates for superconducting quantum bits 
with good energy relaxation and coherence properties. These characteristic 
time scales are expected to be in the order of hundreds of microseconds 
since Josephson vortices are well decoupled from other electromagnetic excitations. 
This weak interaction is due to the vortices’ small magnetic dipole moment and 
the large capacitance of their hosting long junction that serves as an intrinsic filter 
for low-frequency charge noise below the Josephson plasma frequency. Unlike 
established superconducting quantum bits, which are all compositions of Josephson 


junctions, capacitors, and inductors |VD17 , these elements are intrinsically 


included in Josephson media, enabling very compact circuits. 


All these considerations are theoretical, and the experimental implementation 
requires junctions that allow for quantum vortices. The system’s quantumness is 
governed in essence by the junction impedance, i.e., depends on the inductance of 
the electrodes, the capacitance of the junction, and its width. So far, entering the 
quantum regime of Josephson vortices has been experimentally very difficult to 
reach and the coherent quantum regime of Josephson vortices has not 
been realized yet. 


The goal of this thesis is the experimental realization of Josephson vortices with 
a quantum nature in a distributed junction. To increase the long junction’s quan- 
tumness, one needs to enhance its impedance. We use the approach of very narrow 
junctions and replace the junction’s usual bulk electrodes with thin films of a 
high-kinetic inductance superconductor. With this additional parameter of kinetic 
inductance, which exceeds the magnetic contribution by orders of magnitude, the 
impedance can be increased beyond the purely geometrical limit. 


The thesis continues this introduction with a theoretical background on the prin- 
ciples of quantum electrodynamics in superconducting circuits, with particular 
attention to long Josephson junctions. Following this, experimental methods for 
cooling to cryogenic temperatures, the measurement setups, and the fabrication 
techniques developed for impedance-tailored long Josephson junctions are speci- 
fied. 


In the first part of the experimental results, these high-impedance long junctions 
are characterized by transport measurements to observe Josephson vortices in the 
presence of high-kinetic inductance electrodes and to check whether the impedance 
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is increased as intended. In the second results chapter, we investigate a suitable 
manipulation and readout scheme to explore quantum coherence, where a classical 
long junction is embedded in a microwave resonator. In the third experimental 
chapter, we experimentally study whether we can indeed observe Josephson vor- 
tices with properties of quantum particles with the chosen approach. Finally, I 
summarize the main results and conclude with possible future experiments that 
could follow this work. 


2 The principles of quantum 
electrodynamics in long Josephson 
junctions 


This chapter summarizes the theoretical background that supports the understand- 
ing of later presented experiments on high-impedance long Josephson junctions 
(LJJs). As these tunnel junctions are weak link superconductors, basics of the 
superconducting phase are introduced, where I devote special attention to the 
kinetic inductance of circuits composed of granular superconductors that is uti- 
lized to increase the junction’s impedance. After expounding the fundamentals of 
lumped element Josephson junctions that remain valid for distributed versions, 
one-dimensional long junctions are elaborated, as experimentally studied through- 
out this work. Moreover, superconducting microwave circuits are treated, used as a 
novel tool to read out long junctions. 


2.1 Basics of superconductivity 


The successful liquefaction of helium via the Hampson-Linde cycle by Heike 
Kamerlingh Onnes in 1908 opened the research field of low temperatures 
solid state physics. This technology paved the way for studying metals’ electrical 
resistivity near zero temperature. In 1911, Kamerlingh Onnes discovered the re- 
sistance of mercury to drop abruptly to zero below a transition temperature, an 
effect he named “superconductivity” [Kam11]. Except for exceeding this critical 
temperature Te, superconductivity also vanishes above a critical magnetic field He, 
which is also linked with a critical current Ke, generating an intrinsic magnetic field 
(Silsbee’s rule). 


Among many pure metals, superconductivity occurs also in certain alloys, inter- 
metallic compounds, and doped semiconductors and the material specific criti- 
cal temperatures Te range from below 1mK for elementary rhodium over 1.2K 
and 9.2K for aluminum and niobium, respectively, up to 139K for the cuprate 


Hgo Tlo 2BazCaCu3Os_341|[Kit05]. 


1 The given Te values are no physical limit, but merely the extremes yet known at atmospheric pressure. 
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Beyond the superconductor’s zero DC resistance, Walther Meissner and Robert 
Ochsenfeld found an even more fundamental property: Up to a critical magnetic 
field, superconductors expel external magnetic fields from their interior regardless 
of whether the magnetic field is applied before or after the transition to super- 
conductivity [MO33]. This behavior is independent of the path that leads to it, 
consequently distinguishing superconductors essentially from ideal conductors and 
characterizing them as ideal diamagnets. Explaining these results was elementary 
for subsequent phenomenological theories, such as the London theory and the 
Ginzburg-Landau theory, and suggested already to explain superconductivity as a 
macroscopic quantum effect. 


2.1.1 Phenomenological description 


For a phenomenological description of the superconductor’s unique electromag- 
netic behavior of zero DC resistance and ideal diamagnetism, the brothers Fritz 
and Heinz London extended the Drude model of normal metals to su- 
perconductors within the scope of classical electrodynamics. London’s theory 
considers the superconducting state at finite temperatures, consisting of two sepa- 
rate, noninteracting fluids, the normal (index j = n) and superconducting (index 
j = s) charge carriers with either mass mj, charge qj and mean scattering time fo ;. 
Electrodynamics of both fluids is described by the Drude model each 


d 1 m; d 1 : 
E=m,; v= ! j (2.1) 
fj ] (3 =) J nigj E =) Jj 


and the normal conducting static case results in Ohm’s law E = onjn with the 


2 
normal conductivity on = Patia, The superconducting state, however, is charac- 
terized by zero DC resistance, and thus tos —> œ, which is why, the classical Ohm’s 
law is replaced by the first London equation 
djs = ns45 p 
ot Ms 
Equation is Newton’s second law for superconducting charge carriers and 


(2.2a) 


means that electric fields vanish from the superconductor’s interior in the stationary 
state. To examine the magnetic behavior, this first London equation is inserted 
in Faraday’s law VxE = — ƏB and by taking magnetic flux expulsion from 
the inside, according to the Meissner-Ochsenfeld effect, into account, the second 


London equation is found to be 


2 
vxj=- sp. (2.2b) 


Ms 
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These two London equations and can be combined using the electro- 
magnetic vector potential A and the London gauge and result in 
__ nsqsħ nsq2 A. 
Ms Ms 
The additional term (x V@) in Eq. originates from London’s later derivation 
|Lon48] that is based on the postulate of superconductivity as a macroscopic quan- 


(2.2c) 


tum state, described by one single wave function encompassing all superconducting 
charge carriers 


P(r) = polr) ® = y/ng (re, (2.3) 


The Meissner-Ochsenfeld effect is comprehended by these basic assumptions: By 
applying the Maxwell equations to the London equation (2.2ch, the magnetic flux 
density and current density have to satisfy the differential equations 


2 

V2B = Kos, (2.4a) 
S 
2 

V?j= mj, (2.4b) 
S 


respectively. Within a bulk superconductor, the solution of this second-order dif- 
ferential equation describes screening currents underneath the superconductor’s 
surface that decay exponential on the length scale of the so-called London penetra- 
tion depth 


Honsgs 


In the interior of the bulk superconductor, these screening currents compensate an 


applied magnetic field Ha, manifesting ideal diamagnetism with u, = 0. 


Apart from explaining zero direct current (DC) resistance and the Meissner- 
Ochsenfeld effect of superconductors, the London theory can explain the su- 
perconductor’s response to alternating (AC) fields Eac = Ee!®! and consequently 
jac = je’. Within the scope of the two-fluid model, the normal conducting and 
the superconducting charge carriers exhibit a complex net conductivity 


2 2942 
159; toj i nq; tg 
242 2 42 
mj (1 +w D mj (1 +w 5) 


where the real part 7, denotes ohmic dissipation in the normal state, whereas the 


j = jn +js — = (oj — ion j) E, (2.6) 


imaginary part 72 accounts for the superconducting charge carrier’s inertia. 
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2.1.2 Fluxoid and magnetic flux quantization 


Due to their macroscopic wave function, which has to be single-valued for in- 
tegration along a closed path ƏS [Lon50J, Fritz London predicted a “fluxoid” 
quantization in multiple connected superconductors. The wave function of Eq. 
is invariant for phase changes in integer multiples of 27r, such that the closed line 
integral of the phase gradient holds 


rı 


r2 
$ Vedl = lim J vedi = lim [6(12,t) — 6(r,t)] =27n, neZ. (2.7) 
os nor m7 

rı 


Thus, the closed line integral of the London equation yields the fluxoid 


h 
®' = | Bas ag ee 2.8 
4 + HoAL ee! n q (2.8) 


S 
which is Bohr-Sommerfeld-like quantized in multiples of the superconducting 
flux quantum o = h/gs. In case of bulk superconductors, well inside the super- 
conductor, screening currents vanish according to the Meissner-Ochsenfeld effect 
and the fluxoid quantization of Eq. results in a magnetic flux quantization 
[BdS = = n®o. This means that multiple connected superconductors can 
merely enclose integer multiples of the superconducting flux quantum ®o. 


Fritz London assumed erroneously qs = -e, but experiments, e.g., by Robert 
Doll and Martin Näbauer as well as Bascom Sine Deaver Jr. and William 
Martin Fairbank in 1961, showed the superconducting charge quantum 
to be qs = —2e. This supports the idea of Cooper pairs retrospectively, on which 
the microscopic BCS theory bases. The superconducting magnetic flux quantum’s 
current value in the literature is 


h 
Do = 5, = 2.067 833 848 x 1071 Vs. (2.9) 


2.1.3 Ginzburg-Landau theory 


The previously mentioned observation, that the Meissner-Ochsenfeld effect is ther- 
modynamically reversible, reveals that superconductivity is a separate thermody- 
namic state. This motivated Witali Lasarewitsch Ginzburg and Lew Dawidowitsch 
Landau to describe the superconductor’s second order phase transition by Landau’s 
theory of phase transitions [Lan37], extended by the superconductor’s response 
to external magnetic fields [LG50;|GLO9}. The Landau theory relates the system’s 
changing properties with its breaking symmetry that is characterized by the oc- 
currence of a complex, spatially varying order parameter y(r). As per definition 
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of second order phase transitions, this order parameter exists exclusively in the 
less symmetric phase, approaches zero continuously at the phase transition and 
its square turns out to be the superconducting charge carrier density ns. As the 
order parameter is small near the phase transition, the free energy without external 
magnetic field can be expanded into a symmetric power series 


i= frav = [m + acip + PCL iy + o(ipf)av. (2.10) 


Here, fn is the free energy density of the normal conducting state that is assumed to 
be constant. The occurrence of superconductivity as the energetically favorable state 
below Te necessitates acy < 0 and Bor > 0 in order to fulfill the energy lowering 
fs < fn at a finite nonzero order nn (see Fig. 2.1a)). Then, the equilibrium 


=0, yield Bot = — TA with the 


states at the local minima, which hold 


Hily 
order parameter of the spatially piesa equilibrium Poo: 


As superconductors expel external magnetic fields from their interior, their Gibbs 
free energy density gs is increased by the magnetic work — f MdB, in the presence 
of external magnetic fields. In case of perfect diamagnetism (Xm = —1), the mag- 
netization is M = —H = —B,/yo and the field expulsion work becomes B2/2p19 
[GC34]. induced currents are further associated with the charge carriers’ kinetic 
energy rn , where ĵ is the canonical momentum —ifV — qsA derived from the 
principle of minimal coupling. All in all, for superconductors the Gibbs free energy 
density functional reads 


1 2 
geld, A] = gn + acrlpl? + ÉSE yt 4 -|(-inv gA) p? + SE. (2.11) 


2 2m 


The order parameter in the equilibrium is obtained by minimizing the Gibbs 
free energy both with respect to variations in the order parameter dp (while y* 
is regarded independently) and the vector potential dA and results in the two 
Ginzburg-Landau equations 


0 = acLY + bcLlyl y + = — —qsA)* Y, (2.12a) 
je = ine "(pvp yy") — Typ 
= T Re(y* (-ihV — qsA) 2 (2.12b) 


For homogeneous superconductors, the whole ensemble of superconducting elec- 
trons can be treated as macroscopic quantum state with one single wave function 
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Figure 2.1: Ginzburg-Landau free energy and field dependence of the emerging order parameter. 
a) Symmetry breaking at the normal to superconducting phase transition. For acy, > 0 the equilibrium 
state is characterized by the absence of the order parameter = 0, which is the case for T > Te, while 
the occurrence of a finite order parameter Y = Yo below T < Te requires ag < 0. The resulting 
energy lowering (difference of the super- and normal conducting free energy densities fs — fn) can 
for instance be suppressed by applied magnetic fields, and thus overcome by exceeding the critical 
agg? 
HoßGL 
nonlinear supercurrent density js with a maximal (critical) value je. 


y(r) = |e!) according to Fritz London’s postulate |Lon48) and the supercur- 
rent density of Eq. (2.12b) results in 


js = Epl? (HV 0 — gs) = sl pl?vs. (2.13) 


value He = . b) The order parameter is suppressed in strong electric fields, which results in a 


The comparison with the London equation identifies the order parameter 
to be given by the superconducting charge carrier density Ip? = ns. For weak 
perturbation fields, the last two terms of gs in Eq. are negligibly small, so 
that the thermodynamic equilibrium state emerges in the potential minima at 
|p| = |Poo| = \/—acL/BcL, as illustrated in Fig. In case of strong fields, 
however, the superfluid’s kinetic energy Insmg|vs| can no longer be neglected. 
Since the kinetic properties of the superconducting charge carriers will be of 
particular interest for the experiments in this work, they will be considered in 
more detail below. Optimizing the homogeneous order parameter of the Ginzburg- 
Landau Gibbs free energy 


1 B2 
glp A] = gn + ac)? + ÉSE lyt + Fms Pp + 28 (2.14) 


2 This assumption was justified by means of microscopic BCS theory by Lev Petrovich Gor’kov ; 
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acı] 


rent density js = qs|p|vs by inserting in Eq. (2.13). Figure|2.1 b)|shows both the 


reduced charge carrier densities with increased superfluid velocities vs. Depairing 


2 
yields Iyl? z {theo 7 (1 = no and consequently a nonlinear supercur- 


thermodynamically is thus favored, and the maximum supercurrent density is 
Imsv2 


reached for |&gı| = "7. This means that there is a maximal (critical) current 
density 

: 2 

je = 34s|tPool ve (2.15) 


above which superconductivity collapses. With the corresponding critical velocity 
Vc, the supercurrent density results in 


js = [feel govs h — 302 ) ; (2.16) 
Cc 


As mentioned above, the second Ginzburg-Landau equation (2.12b) corresponds 
to the London equation with the characteristic London penetration depth 
in the case of a homogeneous superconducting charge carrier density, and thus 


order parameter. Furthermore, the Ginzburg-Landau theory exhibits a second 
characteristic length scale of superconductors, on which y(r) may vary spatially. 


In the absence of magnetic fields and for small order parameters, e.g., near Te, 
the O (|p|?) term is negligible and the first Ginzburg-Landau equation (2.12a) 
becomes V7) = arlicl, This means that small perturbations dy decay exponen- 
tially on the characteristic length scale, the so-called Ginzburg-Landau coherence 
length 


Pr) 
== 4) 2.17 
EGL TAE (2.17) 


These two length scales, the coherence length Zcı and the London penetration 
depth ^, are main results of the Ginzburg-Landau theory and important quan- 
tities to describe the different properties of superconductors. For instance, the 
ratio ÀL/ČGL defines whether superconductivity collapses abruptly above one 
critical magnetic field H. (type-I superconductor with AL /éc. < 1/V2) or the 
superconductor allows magnetic flux to penetrate partially above the first critical 
field H.ı and becomes completely normal conducting not until the second critical 
field Ho > Ha (type-II superconductor with A/c > 1/V2). In the so-called 
Abrikosov or Shubnikov phase between Ha and Hez of the latter case, flux quanta 
with normal conducting cores get in the superconductor, and thus make for imper- 
fect diamagnetism (—1 < Xm < 0) [DP19]. This second critical field can exceed the 
first critical field by orders of magnitudes [MW18]. 
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The phenomenological Ginzburg-Landau theory covers the superconductor’s 
unique behavior, namely electrodynamics and thermodynamics, using a macro- 
scopic wave function of conduction electrons. However, it provides neither a mi- 
croscopic explanation of their origin nor of the occurrence of superconductivity in 
general. 


2.1.4 Microscopic Bardeen-Cooper-Schrieffer theory 


The microscopic theory of superconductivity by John Bardeen, Leon Neil Cooper 
and John Robert Schrieffer (BCS) heralded a breakthrough in basic understanding 
of superconductivity. This theory is founded on a boson-like ground state of con- 
densed, paired conduction electrons originating from an attractive electron-electron 
interaction, as it would not be possible for single electrons according to the Pauli 
exclusion principle [Tin04}|Sch19}[DP19]. Despite the repulsive Coulomb interaction 
between two isolated electrons, a virtual phonon can mediate a net attractive inter- 
action in the solid body of conventional superconductorsP] This counterintuitive 
exchange mechanism is inspired by evidencing the isotope effect of superconduc- 
tors [Max50}|Rey+50} |Fr650} |Bar50], whose transition temperature depends on the 
isotope mass as Te x 1/,/my and therefore on the phonon frequency. Since phonon 
energies are limited to the Debye frequency fwp, the interaction happens between 
the Fermi energy Er and Er + hwp and is most effective for electrons with opposite 
momentum and spin |k |, —k *). |The energy of these two electrons, a so-called 


Cooper pair, is lowered by approximately [Tin04: 
ee 
AE x —2hwpe "EPY (2.18) 


compared to two “free” electrons of energy 2Er. The Hamiltonian of the whole 
ensemble of superconducting electrons in the volume V holds 


at a» Skk st sto oa A 
HBgcs = L erôkoĉko + L V Cyr pO |C—k,L ke t (2.19) 
ko k,k’ 


where ĉt} „ and ĉk, denote electron’s creation and annihilation operators, respec- 


h? k? 
2m 


tively} The first term represents the electron’s kinetic energy £k = — u with 


It should be noted that this net attractive interaction does not necessarily need to be mediated by a 
virtual phonon. This is only the case for “conventional” superconductors, but for instance Te 2 90K 
of high-temperature superconductors cannot be explained by a phononic interaction ; 

These singlet s-wave Cooper pairs with zero total spin and orbital angular momentum S = L = 0 are 
the simplest case, but for instance singlet d-wave Cooper pairs characterize high Te superconductors 


papag 


ese fermionic creation and annihilation operators fulfill the typical commutation rules 


êkas tte a} = OK K deo! and {é a} = {ko wo} =0. 
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the chemical potential u. The Cooper pairing mechanism in the second term is 
described by the attractive interaction 9,4 of two electrons that are scattered 
from momentum +k to +k’. These Cooper pairs typically extend over a range of 
100nm-1000 nm, overlapping with ~ 10° Cooper pairs. Therefore, the interaction 


oe mee . can Be (Eade aE sachen with a Cooper pair mean-field 


Moreover, the .. wave function on the superconducting system 
simplifies in the mean-field approximation to |Tin04}/Sch19 


(Focs) = TT (un + nek seh) 190) = [I (lel + [orle ek etn) [¥o) - (2.20) 
k k 


Here, lok]? denotes the probability of the singlet pair |k |, —k f) being occupied, 
|u|? = 1 — |v,|* the probability being unoccupied, respectively. The coefficients 
UK, Vk € C can be chosen as ug = |ux| and vk = |v, |e, since the creation operator 
Cepek |, creates one Cooper pair with a certain phase ei? from the vacuum state 
|Yo). This reveals that the BCS ground state is a coherent many-particle state 
in which all Cooper pairs are in a single common quantum state. Due to the 
long-range coherence of the BCS ground state, superconductors can effectively be 
described by the macroscopic wave function of Eq. (2.3), as postulated by London 
and assumed in the Ginzburg-Landau theory (see Sec. [2.1.3). 


As can be seen in the BCS-Hamiltonian of Eq. (2.19), electrons are created and 
annihilate in pairs exclusively, so that their creation and annihilation operators 
need to be decoupled from each other by means of a linear transformation, the 
so-called Bogolyubov-Valatin transformation. The new fermionic operators 


Val58] are given by 


Vt \ _ [uk -% Kt Kt \ _ Hr Uk Vest 221 
vt Not uk Fad êt ~ (ož u e tea) 
Y-kL k k C k} -k,| k "k nue Kl 


In case of phonon mediated singlet s-wave Cooper pairs, the attractive electron 
interaction 9,’ is bound to the Debye energy around the Fermi surface and is 


further isotropic [Col15], such that 


=g, lex, \ew| < hwp A, lex| <hwp 
Sue = $ Ir | \ Ak = bier (2.22) 
0, otherwise 0, otherwise 


The diagonalized effective BCS-Hamiltonian finally holds 
vr 
Bes = = PEK Tk, otk +} (En — Ex) 4 =. (2.23) 
k 


ko 
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The second and third term in Eq. (2.23) denote the constant energy gained due to 
the creation of Cooper pairs (condensation energy). The introduced energy 


Ex = y & + |A]? (2.24) 


emphasizes that A manifests as a gap in the superconductor’s energy spectrum 
even at the Fermi level, where e, = 0. This energy gap and the macroscopic wave 
function are the reason for the superconductor’s vanishing DC resistance up to the 
critical current, since the absence of free states does not allow for scattering with 
an equivalent energy transfer less than 2A. The size of the superconducting energy 
gap varies among different materials and at absolute zero temperature is given by 


h oa 
MD Ohwpe MERE = ne YkgT, ~ 1.764kpT. (2.25) 


^o = — ~~ 
; 1 
sinh (veg) 
with Euler’s constant y ~ 0.5772. Equation (2.25) explains the superconducting 
isotope effect since the superconductor’s gap energy is both directly proportional 


to its critical temperature T. and the Debye frequency. 


Excitations that exceed the minimum energy of 2|A| break a Cooper pair into 
two Bogoliubov quasiparticles, as described in the first term of Eq. by the 
operators Mol] The quasiparticles’ density of states, as visualized in Fig. 
results in 


V dS, Ex 
Nop (Ex) = 53 J = = Na (Er) ———— 
(2m) E=const IVREx| VER -4 


with the normal conducting density of states at the Fermi level Mn (Er) = 


@(E, — A) (2.26) 


m|ke| 
. . . . 272 he ` 
From the quasiparticles’ density of states, a band diagram can be deduced that fea- 


tures no contribution of quasiparticles to the conductance at zero temperature and 
a thermally populated conductance band above the gap for finite temperatures|’| 


6 It is conspicuous that the operators couple electrons and holes, which is why these new fermions are 
called quasiparticles. 

7 Cooper pairs at the Fermi energy Er might be confusing since their energy is lowered by the 
superconducting gap according to Eq. 2.18). In the picture of band diagrams that illustrate fermionic 
densities of states and possible gaps, the boson-like Cooper pairs that condensate in a common 
ground state can not be included in the same manner. To illustrate transport properties, however, a 
combined band diagram of Cooper pairs and their quasiparticle excitations is useful. First, the Fermi 
energy in superconductors at zero temperature is above a fully occupied band and below an empty 
band because there is no fermionic conductance band. Small changes within this gap do not change 
the superconducting behavior, so they can be considered fermionic insulators. However, they show 
perfect conductivity but carried by boson-like Cooper pairs. Second, the energy 2A is required to 
excite two quasiparticles, which were not in any other fermion states of the density of states before. 
Thus, Cooper pairs are below the lower edge of the upper unoccupied valence band by A. 
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Figure 2.2: Quasiparticle density of states and the deduced superconductor’s band structure. a) The 
quasiparticle density of states shows a superconducting gap of size 2A around the Fermi level, where 
there are no states available. However, near above the gap, far more states are possible compared to 
the normal conducting case, whose value N„(Er) is approached far above the gap. b) Against this 
backdrop, a band diagram can be designed for superconductors, which is a useful concept to describe 


the transport behavior of coupled superconductors such as the Josephson effect | NSS60}|/KW19]. 
P P P P 


2.2 Granular high-kinetic inductance 
superconductors 


For the design of superconducting circuits, the basic components, namely resistors, 
capacitors, and inductors, need to be known. As superconductors stand out be- 
cause of their zero DC resistance, resistors rather appear from normal conducting 
quasiparticles. Capacitors are dominated by their dielectric, and thus are not influ- 
enced by superconducting electrodes. However, the inductance of superconducting 
devices may feature additional terms compared to common normal conductors, 
which are addressed in this section. 


2.2.1 Kinetic inductance 


Like any other current-carrying conductor, superconducting strip lines exhibit a lin- 
ear geometric inductance Lg, which depends on their actual shape and dimensions 
and stores energy in the surrounding B-field Tia f B’dv = ILP. This magnetic 
energy, however, is not the entire energy that is carried by the electric current since 
moving charge carriers themselves naturally hold kinetic energy 5 f nsms|vs| dV. 
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As current I and charge carrier velocity are proportionally to the first order, the 
kinetic part contributes to the total inductance, too. In normal conductors, the 
magnetic due exceeds the kinetic due by orders of magnitudes, except for very 
high frequencies f > 10! Hz [SMU97]. In thin superconducting strips or films, 
however, charge carriers can gain kinetic energies comparable or exceeding the 
magnetic energy even for f — 0 because of the substantially higher charge carrier 


mobility [Ann10]. 


The fact that the energy, associated with a kinetic inductance, is spatially confined 
to the superconductor’s interior leads to a material specific constant inductivity 
and a constant inductance per unit square for thin homogeneous films of fixed 
thickness d. In the latter case, the element’s total inductance kinetic inductance 


Ly = L,/N™ scales with its number of squares N- = !/w, where I and w denote 
its length and width, respectively. In contrast, geometric inductances couple to 
the surrounding magnetic fields and depend on the element’s actual shape and 
dimensions but not on its material. This is why the geometric inductance is usually 
referred to the constant inductance in terms per unit length L, = Lg/l. 


By definition, the inductance relates the change in current through and the voltage 


drop across an element as El = Ly dls Ly L wd djs, The electric field E, in turn, 
accelerates the superfluid qsE = dps = ms Is to increases its momentum ps = 


MsVs. This leads to the kinetic square inductance of a superconductor |MT69 
CHS71;|ASB89] 


L [E| /dljsl = _ ms dlvs| (dljs u _ Ms ( djs = 2.27 
N ee ee a. (2.27) 
As already discussed as part of the Ginzburg-Landau theory in Sec. the 
homogeneous supercurrent density is described by the London theory js = nsqsVs 
with ns = Iyl? = Ivo]? in the simplest case of thin superconducting films in weak 
perturbating fields and therefore the kinetic square inductance results in 

ms 1 A2 

neqed Mo 
In case of strong perturbations, however, the js(vs) dependence becomes nonlinear, 
as given in Eq. and calculating the kinetic square inductance of Eq. 
results in 

2 


= 
Lov, Me 150 (1 = =) . 2.29 
x (v ) poo ad ( u 2) k ( ) v2 ( ) 


To use the more measurable quantities, supercurrent and critical current density, 
js(vs) needs to be inverted along its stable branch. The Eqs. (2.15) and (2.16) relate 


(2.28) 


Ly (0) = 
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2.2 Granular high-kinetic inductance superconductors 


Table 2.1: Properties of several high-kinetic inductance superconductors. The kinetic sheet inductance 
iB is both affected by critical temperature Te as a measure of the superconducting gap energy and the 
normal state resistivity pn according to Eq. (2.32}. These quantities are reported from several circuits 
made from different high-kinetic inductance superconductors with comparable thicknesses. 


material Pn Te iD references 
(uQcm) (K)  (pHsq"') 
NbN ~1000 n72 80 

TiN ~450 ~3.0 ~235 Vis+10] 
NbTiN ~200 v9.3 ~35 
MoGe ~220 ~6.0 ~45 
AIO, 6000 ~1.6 ~2000 
js — 3Vs v3 


mom and for small bias currents, where only the first order term is 
considered, the kinetic square inductance is approximated to [ASB89 


Roa ye HPO ( | EY l (E) o(£)); (2.30) 


Apart from this bias dependence, LH (0) « nz! strongly depends on the temperature 
since the superfluid density vanishes at the phase transition at T — Te. The exact 
temperature dependence is individual for different superconductors |Tin04], but 


4 
can be approximated by the Gorter-Casimir two-fluid model at = =1- (x) ; 
Finally, the kinetic square inductance is given by 


ie js) = Lk (0)- G (: | (2) | (2), | o(é)). (2.31) 
(= 


In any superconductor, the kinetic inductance diverges close to Te, where order 


parameter approaches zero, and plays an important role for any superconductor. At 
low temperatures far below Te, kinetic inductance is only relevant in superconduc- 
tors with small superfluid densities. This can be explained descriptively as follows: 
By driving a given bias current, less participating superconducting charge carri- 
ers are accelerated to antiproportionally faster velocities and their kinetic energy 
thus increases quadratically, such that the kinetic energy of the whole superfluid 
increases antiproportionally to ns. 


8 ata film thickness of d = 20nm 
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2.2.2 Mattis-Bardeen theory of granular superconductors 


Due to their suppressed superfluid density, granular superconductor devices are 
propitious to achieve large values of kinetic inductance. In order to describe the ki- 
netic inductance of circuits made from granular superconductors by means of their 
normal conducting properties, it is worth to look at kinetic inductance from another 
perspective of AC circuit analysis. The superconductor’s AC response exhibits 
a negative susceptance Im(Z~') « o (see Sec. [2.1.1}, which corresponds to an 
inductance that accounts for the superfluid’s inertia. Hence, the superconductor’s 


š = Ber à š á wd _ 1 = 1 
complex admittance gs = 01 — ion defines its impedance by 09,77 Z Kar 


Considering the superconductor’s vanishing DC resistance, its impedance is exclu- 
sively given by its kinetic inductance |Z| = wLx, such that this kinetic inductance 


results in Lg = @Im(o,)d' 


In the local (lo < AL), dirty (lo < &) and low-frequency limit (hf < ksT) this 


imaginary superconductivity Im(os) is related to the real normal state conductivit 
On by the Mattis-Bardeen theory as min) = =e tanh (a4). [MB58 Tino4]. 


Thus, the kinetic sheet inductance 


1 hR A(T)\ 70 AR 
L = 2 th n 2.32 
k = wim )d mAT) (3) 7tA(0) oe 
is directly linked to the normal sheet resistance m = pn! = (Rndw/ u E 


(Ra d) 1 To reach large kinetic inductance values, superconductor with small 
superconducting gaps (x Te) and large normal state resistances are preferred. An 
incomplete compilation of different popular granular superconductors as well as 
their kinetic inductance values is given in Tab. 


2.3 Superconducting microwave circuits 


Nowadays, superconducting quantum circuits are commonly manipulated and 
read out in the microwave regime [Gu+17]. Unlike the DC measurement scheme, 
RF coupling can prevent galvanic contacts to the environment that would be 
a substantial decoherence channel. Thereby, microwave resonators are used as 
suitable detectors for coupled quantum systems since their resonant frequency 
shifts dispersively depending on the quantum state of the coupled system 


waron). 
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2.3 Superconducting microwave circuits 


2.3.1 Scattering matrix of a microwave network 


A common technique of electrical engineering to describe a microwave signal’s 
interaction with any circuit, representing an n-port network, is an n x n complex 
scattering matrix, which relates the complex output and input voltage at porti and 


port j by [Poz11] 


yout 
Sij = 


in (2.33) 


J 


V=0,k#j 
This means that S;; denotes the reflection coefficient of the port i, whereas Sj; is 


the transmission coefficient from port j to i. Throughout this work, microwave 
spectroscopy is done at a two-port network, as illustrated in Fig. given by?| 


yp“ Su 52\ (vin 
ts 7 = S22 vin ee =n 


In the present experiments, both ports are connected with a continuous transmis- 
sion line, whose characteristic impedance is typically terminated to Zo = 50 O, and 
this transmission line is shunted to ground by a load with an arbitrary impedance 
Zı (see Fig. P-3). For this so-called notch-type configuration of a parallel scatterer, 


Kirchhoff’s laws yield the scattering matrix |Poz11 


wi (G8) +2), 
gnote +3)" (1422) . (2.35) 


To give a vivid description, the abrupt impedance mismatch in the transmission 


line due to the external load yields a part of the incident signal being reflected 
and the other part being transmitted. Thereby, complete reflection is obtained for 
Zı — 0, complete transmission for Z} — œ, whereas the signal is split in equal 
parts in case of Zı = Zo/2. 


2.3.2 Superconducting microwave resonators 


One of the simplest loads featuring quantized excitations are quantum harmonic 
oscillators in the form of RLC-circuits. In this work, such RLC-resonators are 
coupled to a microwave feedline on the one hand and to long Josephson junctions 


? The matrices ST = S of the circuits studied in this work are symmetric, which means that interchang- 
ing the ports 1 and 2 will not vary the system’s behavior. The whole measurement setup where these 
samples are embedded, however, include active components (amplifiers) and ferrites (isolators) that 
break the time reversal symmetry so that the S-matrix is no longer symmetric. 
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a) 


Figure 2.3: Schematic of a half-wavelength notch-type resonator embedded in a two port-network. 
a) The particular notch-type configuration consists of an LC transmission line with specific impedance 
Zo = VL/C that is shunted to ground via a load with impedance Zı. In the present case, this load is a 
terminated transmission line whose two open ends define a distributed A/2 resonator. The two-port 
network of this assembly is characterized by relating the incident and scattered microwaves at both 
ports of the transmission line. The coupling between feedline and resonator is assumed to be capacitive 
with loss rate x, and the resonator itself has an internal linewidth xi. b) As a schematic circuit, the 
infinite transmission line is modeled as an inductive waveguide that is capacitively shunted to ground. 
The distributed RLC-resonator has a well-defined length, whose impedance mismatch at the open ends 
forms a resonating cavity. 


on the other hand to investigate single photon excitations in these junctions. In 
order to probe these extended junctions, open-ended distributed A/2 resonators 
are suitable since their current distribution is defined by the excited mode. In the 
center part of the resonator, A/2 resonators exhibit current nodes and antinodes for 
odd and even harmonics, respectively, that probe the junctions differently. A A/2 
resonator itself is formed from a waveguide of a certain length 1, which constrains 
resonating wavelengths to integer multiples of A/2. Neglecting the frequency shift 
induced by the coupling capacitors, the resonant frequency of the n!* harmonic of 
a A/2-resonator holds 
1 
DIYTICH 


The fundamental frequency wł depends on the inductance per unit length L’ and 


(2.36) 


D i Lo 
wy =N: wj =n 


capacitance per unit length C’. Since this quantity does not account for dissipation, 
quality factors are introduced, describing the ratio between the total energy stored 
in the resonating system Etot and the energy dissipated per cycle Egiss = Pioss/w"- 


r 
ee (2.37) 
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Here, x denotes the inverse photon lifetime in a linear resonator, meaning that the 
excited resonator decays into its ground state at a rate x/2. In the frequency domain, 
« corresponds to the full width at half maximum of a resonator’s Lorentzian shaped 
amplitude signal. For the whole load, there are two main loss channels: an internal 
due originating in ohmic dissipation and dielectric losses and a coupling due to 
other systems as e.g., the transmission line. For that reason, the total quality factor 
of the whole load Q; composes of an internal quality factor Q; and a coupling 
quality factor Qe. Since decay rates x add linearly, the loaded quality factor is given 


by 
ed KI = Ki + Ke- (2.38) 
In the notch-type configuration, the load’s frequency-dependent input impedance 


can be expressed in terms of the internal and coupling quality factor as well as the 
frequency detuning from the resonant frequency 
r 
plas) = Zp (2 -i E ) , 
The S-parameters of ideal notch-type resonators can be found by inserting Eq. 
into Eq. (2.35). Since the scattered signal propagates to both ports, S11 + S21 = 1 
holds. In the following only 521 is treated since this is the later measured matrix 


(2.39) 


entry and as it implies S11 = 1 — S21. The complex signal is usually given in the 
polar form with amplitude and phase 


gnoteh _ (1 Qı/Qe | 1-5, (2.40a) 
1+ 2iQ 2 
2 — (Q; — Q.) / Q2 

gpoteh]? = q _ 1- Q-Q) /Qe (2.40b) 

14407 (5) 

= 2010. w- Wy 

arg (spotch) = — arctan (2 = 2) — arctan ( ZSIS ji 2.40c 
(s2) a, OO a nn 


Figure [2.4/shows both the complex plane, the amplitude, and the phase of S21 of a 
notch-type resonator. In the complex plane, the S21-parameter manifests as a circle 
with diameter Q;/Q, around the center, located at 1 — Q);/2Q,. The Lorentzian 
shaped squared amplitudes |S; (w) j? correspond to the Fourier transform of an 
exponential decay in the time domain and exhibit a full width at half maximum 
of x). As can be seen in Fig. [2.4| the interplay of Qi and Qe plays a key role in the 
response of resonating systems that are usually categorized in three cases: 


e under-coupled (Qi < Qe): Here, most of the energy is lost internally resulting 
in a weak signal both in amplitude and phase, and thus in a circle with a 
small diameter and a large distance from the origin 
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Figure 2.4: S21-parameter of a notch-type resonator for different coupling regimes. a) In the complex 
plane, frequency sweeps around the resonant frequency feature circles, starting and ending at the 
off-resonant point S21 = 1. The the diameter of the circle increases with the ratio of Q;/Q. and tends to 
1 for Qi /Q. — œ (dashed). b) The squared amplitude shows Lorentzian dips, the deeper the larger 
Qi/ Qe and c) the phase signal maximally rolls off from —7t/2 to 7/2 for QGi/ Qe > ov. 


e critically coupled (Qi = Qe): Here, half the energy is lost internally and the 
other half by coupling such that the amplitude on resonance equals |S2;| = 0.5 
(—6 dB) 


e over-coupled (Qi > Qe): Here, most of the energy is lost during coupling, 
resulting in the strongest signals with maximal diameter 1, zero transmission 
amplitude and a phase roll-off from — 7/2 to 7/2. 


The complex-valued S2; data are experimentally acquired with the two-port vector 
network analyzer (VNA) and analyzed with a circle-fit routine (see Fig. [2.5), which 
provides values for the quality factors Qi, Qi, and Qe [Pro+15]. For this purpose, 
Eq. for an ideal notch-type resonator is adapted to the experimental setup. 
The environment is considered by an overall factor ae'*'e”""e that includes at- 
tenuation and amplification of the whole signal a, a global phase offset a9, and 
a cable delay Tę denoting the signal’s traveling time through the setup. Possible 
impedance mismatch is further taken into account by an asymmetry angle &z that 
rotates the circle’s center around the origin. This means that the resonance point is 
no longer on the real axis of the complex plane and the squared amplitude does 
not follow a clean Lorentzian. Having regard to these corrections, the transmission 


coefficient holds [Kha+12 Sch20] 


gnotch — goitogiwta [1 _ Q Qe \ (2.41) 
1+ 210158 


22 


2.3 Superconducting microwave circuits 


6.46450 6.46475 6.46500 6.46525 
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Figure 2.5: Measured and fitted S>ı-parameter of a notch-type resonator. a) In the complex plane, 
the circle is rotated around the origin by the angle ¢, which accounts for impedance mismatches. 
b) This further manifests in the asymmetrically distorted Lorentzian amplitude, whose right shoulder 
exceeds the left one. All the amplitudes are offset corrected by a global factor including attenuation 
and amplification of the measurement setup. c) The phase signal is shifted by the electrical delay and a 
global phase offset of the setup. The circle fit routine yields the quality factors Q} ~ 51300, 
Qe ~ 64400, and consequently Q; ~ 250500. 


Power dependence of high-kinetic inductance microwave resonators 


Microwave resonators with a substantial kinetic inductance feature significant 
power dependence. On the one hand, the resonant frequency, given in Eq. 
is dominated by the nonlinear kinetic inductance of Eq. and consequently 
power-dependent, too. Thereby, the relevant scale is the induced current flowing 
through the resonators with respect to the superconductor’s critical current Te. This 
nonlinear resonant frequency also affects the whole S21 spectra for high powers on 
the other hand: While sweeping the probe frequency close to and even through 
the initial resonance, the resonator becomes more and more excited, which in turn 
increases the kinetic inductance, and thus lowers the effective resonant frequency 
[Swe+13]. This behavior is often called “bifurcation” and exhibits interrupted 
circles in the complex plane and abrupt jumps in amplitude and phase, as shown in 
Fig.|2.6] As a consequence, the S21-parameter that is analytically given in Eq. 
is only valid for small probe powers, meaning few photons in the resonator. 


The average photon number is a suitable power scale, especially with regard to 
quantum systems, interacting with single photons. Therefore, the total energy 
in the definition of the internal quality factor Q; of Eq. is considered as 
Etot = (photon) hw". In the notch-type configuration, the power lost in the system 


23 


2 The principles of quantum electrodynamics in long Josephson junctions 


= jö P (dBm) 
—-108.9 
0.0 0.00 ==. 95.9 
92,9 
2 025 — 90.9 
= 02 £ 
a > 
E | A —0.50 
—0.4 bo 
0 = —0.75 
—0.6 —1.00 
6.4645 6.4650 
f (GHz) 
0 
6.4650 
a 
= 
T 6.4649 = X 
= N 
a an 
Ā 6.4648 E 
—10 5 
6.4647 E 
6.4645 6.4650 —110 —100 —90 
f (GHz) probe power P (dBm) 


Figure 2.6: Power dependence of high-kinetic inductance microwave resonators. The resonance fre- 
quency of devices with dominating nonlinear kinetic inductance strongly depends on the drive power 
until the resonance vanishes. This entails bifurcating S21 spectra, featuring a discrete jump both in 
amplitude and phase as well as a discontinuous circle in the complex plane. The given powers bear on 
the feedline to which the notch-type resonator couples. 


can further be related to the input power and the reflection and transmission 
coefficients Pjoss = Pin — Pout = Pin (1 — |S]? — Sl). On resonance (w = w”), 
the S-parameters of Eq. (2.40) become S$ 1(w") = 1- Sıı(w') = 1 — St and in 
consequence, the average photon number in a notch-type resonator holds |Sch20]| 


2Pin Q? 
(photon) = har? Os (2.42) 


2.3.3 Coupled resonant systems 


The load, coupled to the transmission line in notch-type configuration, as illus- 
trated in Fig. naturally can not only consist of one microwave resonator, but 
also of two or even more coupled resonating systems. Such arrangements are 
of particular interest, if the individual resonant frequencies may be tuned by an 
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Figure 2.7: Avoided level crossings of two coupled resonant systems with constant coupling strength, 
but different dispersion relations. If the two individual systems are strongly detuned, they are mostly 
unaffected by the other, such that their resonant frequencies are shifted only weakly. The closer to bare 
eigenfrequencies (dashed gray lines) come to each other, the larger the new eigenfrequencies (solid 
lines) repel each other. Where both bare frequencies would coincide, the frequency difference w4 — w_ 
takes its minimal value of 2g. 


external parameter, as for instance by magnetic or electric fields, such that their 
bare resonant frequencies would cross each other. However, if the systems can 
exchange energy at a constant rate g, their oscillation modes will hybridize when 
wt z wir 


their resonance frequencies are close, ~ g. The interplay between these 


two coupled systems, each oscillating with a;(t) = aje t, is described by the 


eigenvalue equation [Nol18 


Ww wt 


8 ay(t)\ _ 
> en) (2.43) 


and results in the new eigenfrequencies 


rl r,2 wtl — w2)? + 4g2 

pa „ví Bia (2.44) 
2 2 

These solutions are schematically depicted in Fig. 2.7Jand reveal an avoided level 

crossing, also called “anticrossing”. 


w 


The approach of Eq. provides merely information about the eigenfrequencies, 
but not on the whole spectra including all specific linewidths. Within the scope 
of the input-output theory [GC08], the S21-parameter of a resonant system with 
bare frequency w"? and internal linewidth «;2, coupled to a notch-type resonator 
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Figure 2.8: Calculated S21 spectra of avoided level crossings at different internal quality factors. The 
first resonator is chosen to be critically coupled, its linewidths are kept fixed at x; ı/2rr = Ke/27 = 
10 MHz, and the coupling strength between both systems remains constant at g/27 = 50 MHz, too. 
From left to right, both amplitude and phase signal blur out, during the transition from strong to weak 
coupling by increasing the internal linewidth of the second resonator xj. 


(coupling constant 2g) with its resonant frequency w"!, internal, and coupling 


linewidth x; 1 and xqMlis given by zan 


5 =1 A 2.45 
21 (w) T E aD | ag2 (2.45) 


2i(w—w"??) +1; 9 

Both amplitude and phase of the S2|-parameter are plotted in Fig. [2.8| depending 
on the internal linewidth of the second resonator xj 7. In the case of strong coupling 
Ki2 < g (left), the two new eigenfrequencies are clearly visible and correspond 
to Eq. (2.44), but the avoided level crossing becomes more and more indistinct 
with increasing losses in the second system. In this case of weak coupling, the 
second resonator is nothing but an additional loss channel for the first one since 
its excitation decays faster than it can couple back to the first resonator. This 
mechanism further accounts for the line widening of the first resonator in case of 
weak coupling (right). 


10 Note that all linewidths «x are defined as full width at half maximum, whereas the coupling constant 
g is the half width at half maximum for reasons of convention. 
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2.3.4 Coherence and relaxation properties of a two-level 
system 


In general, not only two harmonic systems, but any kind of resonant system can 
be coupled to a harmonic oscillator. In particular, superconducting resonators are 
widely used as detectors for coupled quantum systems and to read out their state. 
In the simplest case, such a circuit is an effective two-level system comparable 
to a spin 1/2 particle, having a ground state |0) = (0,1)! and an excited state 
|1) = (1,0)' orthogonal to |0). Thus, every possible state of this “quantum bit”, 
often abbreviated as ”qubit”, can be expressed as a superposition thereof 


|%) =a|0) +b]1). (2.46) 


Since |) is normalized to (¥ | Y} = 1, the complex amplitudes must hold la]? + 
|b|? = 1, where |a|* and |b|? describe the probabilities that the system is in the 
ground and excited state, respectively. The global phase of |Y) has no physical 
meaning, so a € R can be chosen to be real-valued. In spherical coordinates with 


the azimuth angle 0 and the polar angle ¢, an arbitrary state reads [Kra+19}|Sch20 
0 i 0 
F) = cos( 5) |0) +etsin(5) [Iy (2.47) 


and is on the surface of a sphere, the so-called Bloch sphere [Blo46]. For historical 
reasons, |0} is aligned in positive and |1) in the negative z-direction. In this picture, 
any qubit operation can be described as linear combination of the Pauli matrices dy, 
dy, and 6, that correspond to rotations of 7t around the respective axis. For example, 
creation and annihilation operators that excite the ground state ô+ |0) = |1) and 
deexcite the excited state ö_ |1) = |0) can be defined as 

ĝ = ; (ôx + iö,) . (2.48) 
In terms of the qubit transition energy Eo, = hwo, the qubit Hamiltonian is given 


by 
hwo. 
Haubit = T Oz. (2.49) 


The time evolution of the pure qubit state given in Eq. (2.47) is governed by 


i 


| (t)) =U(E, to) [F (to = 0)) = exp (= Haus (t- to) ) [F (to = 0)) 


= cos (5) |0) + eilP wort) sin($) hoe (2.50) 


using the time evolution operator U and the qubit Hamiltonian of Eq. (2.49). This 
time evolution corresponds to a continuous rotation around the z-axis with the 
transition frequency wy, i.e., a Lamor precession. 
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Figure 2.9: Drive frequency and amplitude dependence of Rabi oscillations. The qubit population 
oscillates with the pulse length of a photon field, where the frequency strongly depends on the drive 
frequency and amplitude (see Eq. (2.52). a) This Rabi frequency is maximal for zero detuning Aw, 
meaning the drive frequency wa exactly matches the qubit’s transition frequency woı. Small detunings 
give rise to typical Chevron patterns and the Rabi frequency approaches Aw for large detunings. b) 
For quantum mechanical Rabi oscillations, their frequency increases linearly with the amplitude of the 
drive tone. 


Rabi oscillations 


In order to drive a qubit transition, a microwave tone of amplitude V and frequency 
wa is applied, which couples transversally, and thus in terms of the Pauli matrix 


ôy. Therefore, the Hamiltonian of a driven qubit holds [SN17 Sch20 


h 
H= — Oz +hV cos(wat)dr. (2.51a) 


As the microwave drive itself rotates with the wg around the z-axis, the qubit is 
more intuitively described within the rotating frame of this electromagnetic field. 
Using the unitary operation U = exp(iwatöz/2), the driven qubit Hamiltonian of 
Eq. can be transformed from the reference to the rotating frame and results 
in 


Ei + in (2.51b) 
Here, only rotations of the frequency detuning Aw = woı — wg remain. The di- 
agonalized eigenfrequency corresponds to the generalized Rabi frequency Qr = 
V V2 + Aw? and depends on the frequency detuning Aw = wo, — wa and the drive 


amplitude V. 


A = 
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Analysis of the time evolution in the new basis yields the qubit’s excitation proba- 
bility Pi), Or rather the expectation value of ô; [SN17] 


1 Aw \?* 
Pity = 5 (: = (=) ) (1—cos(Qrt)), (&) = 2ppy — 1. (2.52a) 


Dispersive readout 


Nowadays, qubits are often manipulated and read out via the resonator to which 
they are coupled. In analogy to an atom in a cavity as in optics, the resonator- 
qubit interaction can be described by the Jaynes-Cummings model. There, the total 
Hamiltonian Hjc composes of parts of the bare harmonic oscillator Hr, the pure 


qubit Hqubit and their interaction Hint and reads in the rotating frame [JC63; HR06 
WMoO8]. 


1\ ch 
Hyo = hh (aa ! ) | f 6, the (+4 + o_a*) (2.53) 
2 una oS 
Hr Hqubit Hint 


Here, ât and â denote the photon creation and annihilation operator, respectively, 
wr the resonator’s harmonic eigenfrequency, and g the coupling strength between 
resonator and qubit. 


In general, qubit and resonator exchange energy, which is why the eigenstates of 
the coupled system are neither eigenstates of the pure resonator nor of the qubit. 
This means that a part of the resonator excitation is in the qubit and vice versa, 
enabling manipulation, readout, and decay of the qubit via the resonator. 


If the detuning of the resonator and the qubit is large compared to their coupling 
Aw = wr — Wo; > g, the eigenstates are similar to those of the pure resonator and 
qubit. In this so-called dispersive limit [Bla+04], the system’s eigenstates can be 
approximated as product states of resonator and qubit, while the coupling is only 
a small perturbation. Perturbation analysis of the Jaynes-Cummings Hamiltonian 


(2.53) yields 


= 1 ` 2 
Hyc = ħ (wr + xô) ata 7 (wo ! x) Oz H o( £) (2.54) 


with the dispersive shift x = g?/ Aw for true two-level systems without higher har- 
monics. Since the Jaynes-Cummings Hamiltonian in the dispersive limit (Eq. (2.54)) 
is diagonalized, the effective qubit Čo = wo, +x and resonator frequencies 
Õr = Wr + XÖz equal the diagonal elements. The latter’s dependence on 6; means 
that the qubit state can be determined by measuring the resonator frequency. This, 
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however, is a projective measurement, where the qubit’s wave function (2.47) col- 
lapses either into |0) or |1) with the probabilities cos?(0/2) and sin?(@/2), so only 
Õr = wr + x can be measured. Repeating the same measurement many times to per- 


form statistics yields the expectation value of the qubit state (@,) that corresponds 


to the qubit’s excitation probability as [Sch20 
(6) = (F |ô- |¥) = —a|? + |b|? = — cose. (2.55) 


2.4 Dynamics of long Josephson junctions 


2.4.1 Basics of the Josephson effect 


Key elements of manifold superconducting devices, such as superconducting quan- 


tum bits INPT99;|Wal+00a;|CW08], but also superconducting quantum interference 
devices (SQUIDs) |[Jak+64} |ZTH70 Cla66], and voltage standard circuits 
[Tay+67} |FFT73;|Ham00], are Josephson tunnel junctions as they feature unique 


electromagnetic characteristics. The Josephson effect was predicted by Brian David 
Josephson in 1962 [Jos62] and experimentally proven for the first time by Philip 
Warren Anderson and John Martin Rowell in 1963 [AR63]. 


It describes coherent Cooper pair tunneling between weak link superconduc- 
tors, separated by a thin insulating or normal conducting barrier or constric- 
tions, point contacts, etc. Due to the good controllability of electronic properties, 
superconductor-insulator-superconductor (SIS) junctions find widespread usage 
and are therefore utilized in this work. It is remarkable for tunneling processes 
between superconductors that the tunneling probability of a whole Cooper pair 
equals the one of single electrons, and thus is of measurable order of magnitude. 
This originates in the superconductors’ macroscopic quantum states, whose wave 
functions are reduced from the bulk values inside the tunnel barrier but overlap 
nonvanishing. In a perturbation theory treatment, this manifests by a constant 
coupling energy g on the off-diagonal of the time-dependent Schrödinger equation 


|FLS64a] 


Hi\ [H g pı 
am () = (* 8) (2). m 


Thereby, both superconductors hold the macroscopic wave function Y;(t) = 


J ng (te), as given in Eq. that solves the individual Schrödinger equation 
ihdsp; = Hp; = Ejj with the ground state eigenenergies Ej. Assuming weak 
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tunneling between identical superconductors with ngj = ns2 and fg, = Mia, the 
real part of Eq. (2.56) results in the so-called first Josephson equation 
iy = jesin (p). (2.57) 


Instead of a DC voltage, as expected in case of normal state tunneling, Eq. (2.57) 
reveals that DC supercurrents across a tunnel barrier need to be driven by the 
gauge-invariant phase difference 


2 
T E = | Adl (2.58) 


up to the junction’s critical current density across the barrier je = 2gns/ħ. This 
quantity is a substantial parameter and includes information about barrier proper- 
ties, temperature dependence, etc. 


Within the scope of the BCS theory, Vinay Ambegaokar and Alexis Baratoff gave a 
description of coherent Cooper pair tunneling between two, in general different, 
superconductors [AB63} |Wal76}/BP82} [W0111]. They found the dependency of the 
critical current of lumped element junctions on the normal state tunnel resistance 
Ry and the gap energies A; and A» as 


_ Aı(T) Ai(T)\* R 
k R. K| 4/1 (F) . (2.59a) 


Here, K(x) denotes the complete elliptic integral of the first kind and in case of 
identical superconductors, Eq. (2.59a) simplifies to the well-known Ambegaokar- 
Baratoff relation 
nA(T) A(T) \ T> A(T = 0) sah Te 
= tanh 2.388 x 10°" VK" - —. (2.59b 
em ea (Ee er“ i nn 


Hence, the critical current of the Josephson junction can be estimated by the normal 


conducting tunnel resistance Rn, the critical temperature Te or the gap energy A and 
allows experimentalists to characterize the tunnel junctions at room temperature 
without the elaborate cooling far below T.. From the imaginary part of Eq. 
follows another characteristic dependency, commonly referred to as the second 
Josephson equation. If the eigenenergies of the two weak link superconductors 
are shifted against each other by E2 — Eı = 2eV, this difference in the electrical 
potentials is given by 
_hdg ode 
2e dt 2ndt 
According to the two Josephson equations, applied constant voltages V < 2A/e 


V(t) (2.60) 


imply a rotating phase difference, and thus an oscillating current across the barrier, 
often referred to as the AC Josephson effect. Moreover, the second Josephson 
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equation resembles Faraday’s law of inductance with the phase-flux relation @ = 
we. However, the voltage does not originate in the magnetic energy but in the 
Cooper pair’s kinetic energy instead. In analog to Sec. [2.2.1] this goes along with 
an inherent kinetic-type inductance 


V V V Do Ly 
L(p) = u (2.61) 
ary) yee le cos(p) V 2rtlecos(p) cos(p) 


Depending on the phase difference applied, and thus the current flowing across the 
tunnel barrier, Josephson junctions exhibit a nonlinear inductance with the mini- 
mum absolute value L; = Sr = seh, the so-called Josephson inductance [Wai76]. 
This unique nonlinear inductance makes Josephson junctions key elements for most 
superconducting quantum circuits since it enables, for instance, anharmonic LC 
oscillators with distinct artificial two-level systems utilized as quantum bits. The 
AC Josephson effect, originating in the nonlinear Josephson inductance implies, 
that energy is stored in the lumped element junction. When increasing a bias cur- 
rent from zero (where ọ = 0, too) to a finite value I < Ie, the phase difference has 
to change consistent with the first Josephson equation (2.57). Although the voltage 
will be zero, when reaching the current I and Cooper pairs can tunnel across the 
junction without any dissipation, while this phase change yields a nonzero voltage 
during the current increase. As a consequence, an external current source has to 
do work on the junction and the free energy stored in the kinetic energy of the 


moving superfluid is given by |OD91;\|GMD16 


(2.62) 


From another perspective, this energy resembles molecular binding energy be- 
tween atoms, originating in the overlap of the macroscopic wave functions and the 


increased accessible space for Cooper pairs. Therefore, Ey = Pole is referred to as 


TE 
Josephson coupling energy |Jos65; |BP82]. 


Resistively and capacitively shunted junction model for current biased junctions 


In order to gain an insight into the phase dynamics of lumped element Josephson 
junctions, a useful model was established by Wilbert C. Stewart and Dean E. 
McCumber [Ste68}|McC68]. On the one hand, the Josephson junction’s weak link 
can be considered as two conductors, separated by a dielectric layer, like in a 
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capacitor|!"| On the other hand, the weak link also provides tunneling of both 
superconducting Cooper pairs and (in case of finite temperature or magnetic 
field) normal conducting quasiparticles like a tunnel junction. As a consequence, a 
Josephson junction can be considered as a parallel connection of a resistor with 
current Ip = F, accounting for ohmic quasiparticle tunneling, a capacitor with its 
displacement current Ic = q = cg and the Josephson junction with I; = ksin ọ, 
representing Cooper pair tunneling. For such resistively and capacitively shunted 


junctions (RCSJ) [McC68], [Ste68], the net current reads 
$e aw ng EA DcCdp Do dp 


di += g tsin on d2 + IAR dt ` ksin @. (2.63a) 
In e units Eq. (2.63a) can be rewritten as 
d?o 1 dọ 
= +4 7 +sing, 2.63b 
y IR Jp: dt + sin Q (2.63b) 
where y = Í denotes the normalized bias current and the time is normalized 


to F = wpt in nie of the inverse plasma frequency of the junction w p|] This 
characteristic frequency describes the frequency of small phase oscillations at zero 
bias current and equals the resonance frequency of the LC-circuit formed by the 
Josephson inductance Ly and the junction capacitance C 


2rtl. 1 
Wp = = : 2.64 
p Jz YLc ee) 


Moreover, the so-called Stewart-McCumber parameter be = z R2CI. is introduced 


as a measure for damping. The dissipation of a Josephson junction is caused by 
ohmic quasiparticle tunneling, whose corresponding time constant equals Trc = 
RC. Hence, the quality factor of a Josephson junction is given by |Tin04]| 


Q = wptrc = Be: (2.65) 


Equation (2.63) can be considered as a classical equation of motion for a fictitious 
phase particle of mass mg = C(®/27)? ne 2 well comparable with a damped 


gravity pendulum, described by më — (x = —S- U In analogue, the potential energy 
of a current biased Josephson junction in the RCS] model can be identified as 
U(p) = -Ej (vp + cos ẹ), (2.66) 


11 The existence of a junction capacitance is in general not necessarily given (for instance in 
superconductor-normalconductor-superconductor junctions), but the resistively and capacitively 
shunted junction model is an adequate description of dielectric tunnel barriers such as in SIS junc- 
tions used throughout this work. 

12 The term “plasma” frequency is due to its dispersion relation, which equals plasma oscillations in 
metals, as can be seen later in Sec.P.4.3] 
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which is often called “tilted washboard potential” and is illustrated in Fig. [2.10 a)| 
Without any bias current y, the potential energy is cosinusoidal, but applying 
bias currents tilts the potential proportionally to y. This yields local maxima at 
Pp = 7 —arcsin(y) and local minima at go = arcsin(7), giving rise to metastable 
states of the phase particle. The barrier height in between is 


Uo = U(gp) — U( go) = 2E (v -%7 — varccos(1) ) (2.67a) 
z I (1-7)? 4 Sr (1-7)2+0 (a 7°) . (2.67b) 


The typical experimental situation E; > kgT implies that thermal excitation be- 
comes relevant for y — 1, the center of the Taylor expansion, given in Eq. (2.67b). 


Due to its bias current dependence, the potential wells are in general anharmonic 
and current-dependent, too. By Taylor expanding the potential around the local 


Up) 


curvature in the well 


ag 


minima, the frequency of small amplitude oscillations is related to the potential 
a and the oscillating particle’s mass m, by [Dah+68 
=P0 


= = wp (1 = ry (2.68) 


Beyond these energetic and dynamic considerations, the RCSJ model provides 
information about Josephson junction’s time averaged IV characteristics (IVC). 
Within the scope of the tilted washboard potential, the IVC shows a vertical slope 
in the superconducting branch below I. due to (¢) = (V) = 0 (zero-voltage state) 
and a finite slope Ry 1 in the normal conducting branch above Ie, where the current 
cannot be carried completely by the supercurrent through the junction but also by 
dissipative quasiparticles. 


The IVCs vary for different Stewart-McCumber parameters, as can be seen in 
Fig. In the overdamped case ße < 1, even though the phase particle can 
move from one valley to another for I, Z Ic, it does not achieve the running state. 
This is because the released energy is dissipated, and the phase particle needs to 
be activated in every single valley. Hence, the rate ad is small and the resulting 
voltage does not jump directly to the gap voltage Va = #4. In the underdamped 
case b. > 1, however, the phase particle can run down the full potential after 
overcoming the first potential well since the gained energy (minus the dissipated 
energy that is small because of low damping) suffices to overcome the next potential 
well. As a consequence, this running state maintains, even if the bias current I, 


34 


2.4 Dynamics of long Josephson junctions 


a) b) 104 
1.5 
10 & 
= > 8 
> 5 1.0 10? = 
b 8 it § 
E 3 5 
T N 0.5 10° S 
Š F = 
a E | 
£ 10 5 
z 
u 
0.0 Ù 
107? 
0.0 0.5 1.0 1.5 
phase difference normalized voltage V / kR 


Figure 2.10: Potential energy for varying bias currents and numerically simulated IV-characteristic 
for different damping values within the scope of the RCSJ model. a) For zero bias current, the 
phase particle is trapped in a minimum of the potential well with potential height E; and oscillates 
inside with the specific plasma frequency wp. Since wp is typically in the order of tens to hundreds of 
gigahertz, in transport measurements, only time averaged values are measured. These mean values 
are time-independent ($) = 0 and according to the second Josephson equation no voltage drops 
across the junction. However, an applied bias current tilts the washboard potential and if it exceeds 
the critical current I, > I, the phase particle can overcome the potential well with a finite phase 
velocity (p) # 0 which results in a voltage drop across the junction. b) A damping-independent 
superconducting behavior characterizes the up-sweep of the bias current I, up to the critical current Ic, 
where the voltage drops abruptly to the superconducting gap voltage V4 in case of low damping (dark 
line color). Beyond an underdamped junction shows an ohmic dependency due to ohmic quasiparticle 
tunneling. With increasing damping (brighter line color), this mentioned voltage drop decreases, and 
normal conducting behavior approaches the linear ohmic slope slowly and more slowly. Furthermore, 
the hysteretic behavior between up- and down-sweep of the bias current increases with decreasing 
damping (from dark blue to yellow), meaning the IVC follows the ohmic straight below I, down to a 
retrapping current I, where it becomes superconducting again. 


is decreased again down to a retrapping current I, < Ie. The resulting hysteresis 
and, in particular, the ratio of retrapping to critical current is a measure of the 
Josephson junction’s quality factor, given in Eq. (2.65). 


Although the RCSJ model describes significant parts of the junction’s IV charac- 
teristics correctly, such as the occurrence of a zero-voltage state and a damping- 
dependent hysteresis, it is essential to recognize that real Josephson junctions do 
not follow the ohmic branch during the down-sweep, but remain almost constant at 
the gap voltage V4 (see Fig. 2.11 a)). At this branch, the RCS] model is limited by its 
assumption of an ohmic, voltage-independent quasiparticle tunneling resistance. 
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Figure 2.11: IV-characteristic of a current biased and underdamped Josephson junction. a) IVC of a 
110 um x 5 um wide underdamped Josephson junction with electrodes made from AlO, and pure Al 
measured at T ~ 20 mK. b) The simulated IVC by means of the nRCSJ model in Eq. exhibits an 
almost vertical slope (conductance) around the gap voltage at the down-sweep between the critical and 
the retrapping current. c) Various conductance on different branches of a real junction’s IVC within the 
scope of band diagrams of two weak link superconductors [Poo+14]. (i) Below the critical current, 
the junction remains in the zero-voltage state, where superconducting (s) Cooper pairs may tunnel 
coherently and without any dissipation (s — s). (ii) At the critical current, the junction jumps to the 
gap voltage 2A/e, while Cooper pairs are broken in two quasiparticles (q) each. These fermions may 
tunnel through the barrier dissipatively with a nonlinear resistance due to the quasiparticles’ increased 
density of states (s — q). (iii) Far above the critical current, and thus way beyond the gap voltage, 
lots of normal free electrons (n) tunnel yielding the ohmic normal state resistance (n — n). When 
decreasing the bias below the critical current underdamped junctions remain in the running state at the 
gap voltage down to a retrapping current (analog to (ii)). (iv) Close to the retrapping current, and 
thus below the gap voltage, a small number of not yet recondensed quasiparticles may tunnel with 
dissipation (q — q), which is often referred to as Stewart-McCumber branch. 
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In real Josephson junctions, however, this resistance is much smaller for voltages 
close to the gap voltage Vı than for those far above. This nonlinear resistance 
around the gap originates in the drastically increased quasi-particle density (see 
Fig. that are available for only few thermally activated quasiparticles 
(PS73}|BP82}/KW19]. This behavior is covered by a nonlinear resistive and capaci- 
tively shunted junction (nRCSJ) and can be modeled by the power law expression 


[MK83] 


n>1 (2.69) 


and results in a nearly constant gap voltage in the down-sweep, similar to experi- 


mental data (compare Fig. |2.11 a)}and|2.11b)). 


2.4.2 Magnetic field and size effects of distributed junctions 


For all the former mentioned phenomena of Josephson junctions, they were so far 
taken as point like “lumped elements” with spatially homogeneously distributed 
supercurrent densities js(x) = js and phase differences p(x) = g. Real SIS junc- 
tions, however, involve a spatial extension along x and y of length £ and width w, 
respectively, which may feature spatially varying js(x,y) and (x,y). 


Low-inductance limit of bulk electrodes 


In order to work out the phase difference introduced between two positions, say, 
along the x-coordinate, a cross-section of a distributed junction with bulk electrodes, 
as sketched in Fig. [2.12} is considered. In this scheme, two closely spaced pairs of 
points in the vicinity of the junction are regarded and the closed contour C through 
these four points can be split into the two paths across the junction and these along 
the electrodes 


= $ Adl = ie Adl + fae Adl + A Adl + ik Adl. (2.70a) 
eee eee” 
ae a et Fe. el. junction 2 top el. 


Well inside bulk electrodes d > AL screening supercurrents j are negligible small 
and the first London equation reduces to A = Su v9. Consequently, the phase 
drops along the bottom and top electrode are described by 


2a pri by „t by „t 
m. In. Adl = 0(r8\r}) —6(rb\7). (2.70b) 
12 


37 


2 The principles of quantum electrodynamics in long Josephson junctions 


a) 2 b) 27 


pat 
i] 
I 
Basen 
(x) 
© a 


| 
T 
B(x)/ ( 


AL —1 


oe 


AL 


@------- 
255252555 


Figure 2.12: Spatial phase variation in distributed junctions and the Josephson vortex solution. a) The 
Josephson junction is supposed to consist of both a bulk bottom and top electrode that are separated 
along the xy-plane by a tunnel barrier layer of thickness dj. An external magnetic field is applied 
along the y-axis and penetrates the bulk electrodes on the scale of the London penetration depth AL. 
b) Gauge-invariant phase difference, magnetic flux density, and current density of the Josephson vortex 
solution for the Ferrell-Prange equation. The magnetic field and the Josephson current density decay 
with the characteristic length scale Aj. The current’s opposite signs indicate the ring current associated 
with 277 phase kink, whereas the net current is zero. The integration of the magnetic flux density B 
across the vortex’s cross-section equals a total flux of the superconducting magnetic flux quantum ©. 


With the definition of the Josephson junction’s gauge-invariant phase difference in 
Eq. (2.58), the phase drops across the junctions result in 

Qn pN Bed a oe 

= I. Adl = e(r? \r5) — o(ri\r) F 9(x1\x2), (2.70c) 
depending on the direction of walking the contour. The enclosed magnetic flux 
through this the selected area is ® = uoHy - (x2 — x1) A with the magnetic thickness 
A = M1 + Aı,2 +4 [Jos64]. For an infinitesimal small spacing dx in the x-direction, 
the magnetic field results in 

Do am PR) Pi) _ Po dpa) „AP 

= l = ‘ 2.71 

(x) 2rruoA Ka X2 — X1 2rupA dx i dx G7) 


The supercurrent along z, driven by this phase gradient, is on the one hand linked 


with the magnetic field by Ampere’s law jz = oo and on the other hand defined 


by the first Josephson equation (2.57) jz(x) = jesin(p(x)). In combination, the 
phase holds a stationary, unperturbed sine-Gordon equation that is often referred 


to as Ferrell-Prange equation |FP63 


Do do d?o i 
TTA =F = sin(@). (2.72) 
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With the electrode’s sheet inductance Lg = po^ [Swi61 Lik86], the junction’s 


characteristic length scale is given by 


Do Do 
À = = 2.73 
J i mat Te ( ) 


that is often called as Josephson penetration depth. In dimensionless units, the 


length x is normalized to Aj so that ¥ = x/Aj. The general solution of this Ferrell- 


Prange equation is for instance given in |O567}|/LS67]. From a physical perspective, 
the most interesting particular solution in linear geometry junctions in the limit 


£ — œ are Josephson vorticed!5| 


gr (x) = +4arctan(e* *°) +2rın. (2.74a) 


whose center of mass is located at Xo. As can be seen in Fig. this solution is 
associated with a change of the phase difference across the junction of 27t, which 
is why it is often referred to as “27t-kink solution”. This twist extends over the 
characteristic length scale of the Josephson penetration length and is centered 
around Xp. 


From Eq. (2.71), the magnetic flux density that goes along to this solution results 
in 
Po Opp OX Do 1 


B x) = =x ; 2.74b 
rlt) = an OF dx AA; cosh(¥ — žo) N 
The corresponding magnetic flux equals the superconducting magnetic flux quan- 
co A/2 = a ` ` K . 
tum Dp = [_ „dx Ar dzBr,(X) = Bo (see Fig. |2.12b)), which is why this solution 


is also often called “fluxon”. 


According to Ampére’s law, the associated supercurrent 
Se 1 OBry ax o — sinh(x — Xo) . sinh(x — Xo) 
je,(X) = 32 Ze 2 7a a = 2je 2/2 = 
uo OX 9x TAAS cosh“ (X — Xo) cosh* (* — Xo) 


comes as a ring current around the tunnel barrier, which is why the denotation 
“Josephson vortex” solution. It should be noted that in contrast to Abrikosov vor- 


(2.74c) 


tices in type-II superconductors, Josephson vortices feature no normal conducting 
core. As superconductivity is weakest in the junction region, the fluxon’s energy is 
the lowest around the junction, which is why they enter there and are pinned to 
the tunnel barrier (see Fig. (2.16). 


13 Richard A. Ferrell, and Richard E. Prange quoted the solution p(£) = 2 arcsin[sech(% — %)]| that is 
an equivalent representation for a semi-infinite Josephson junction (0 < ¥ < œ) which fulfills the 


boundary condition g(co) = 0 : 
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High-kinetic inductance limit of thin-film electrodes 


In Josephson junctions with electrodes made from high-kinetic inductance super- 
conductors, the magnetic penetration depth AL ~ 3m exceeds the typical film 
thickness d ~ 20nm, which yields incomplete screening of applied magnetic fields. 
Therefore, the assumption of vanishing screening currents well inside the electrode, 
made in Eq. for the derivation of the Ferrell-Prange equation, is no longer 
valid, and this case needs to be discussed separately. 


By solving the Poisson equation of screening currents (see London theory in 
Sec. under the boundary conditions of an extended Josephson junction by 
means of Fourier transformation, Alfimov and Popkov worked out that the phase 
difference across a junction with thin-film electrodes is given by 


. 1 se x —u|\ qd? 

sin(p) = z. ca(! 7 a (2.75a) 
es tanh (dv1 +F) i 

Ga) = = i i a a (2.75b) 


Here, d = $ denotes the normalized electrode thickness with respect to the 
London penetration depth and A = a the ratio of London penetration depth 
and the bulk Josephson length. In Be limit of large Ay, meaning A < 1, as it 
is the case throughout this work, the kernel of the integral operator, given in 
Eq. simplifies to lim;_,o 16,( 54) = tanh(d)d( — u) with the Dirac 
delta distribution ô and Eq. result in the spatial phase variation 


œ 2 2 
sin(p) = / tanh(d)5(z— Lau = tan). (2.76) 
This result is a local sine-Gordon equation with the effective Josephson length 
d 
Agi = ap i tanh (+) . (2.77) 


The equal result is obtained from the Josephson length’s definition in Eq. (2.73) by 
inserting an effective magnetic thickness in junctions with thin-film electrodes, as 


derived earlier by Weihnacht [Wei69] 
m | = + dy (2.78) 
tanh( 7) tanh ( #2) 


with identical bottom and top electrode and negligible tunnel barrier thickness dy. 


Aeff _ 


The Josephson vortex solution (2.74) remains valid for Eq. (2.76), though the char- 
acteristic length is Agi instead of AP |AP95]. It should be noted, that this solution 
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still features a phase winding (fluxoid) quantization in 271, which, however, does 
not necessitate quantized magnetic flux anymore. As the major part of the vortex’s 
total 277 phase winding drops at the dominating kinetic inductance that does not 
couple to magnetic fields and implies incomplete screening, the magnetic flux 
possessed by a Josephson vortex ® is thus significantly smaller than the magnetic 
flux quantum ®o. Therefore, this kind of vortex should be more correctly denoted 


as “fluxoid” instead of “fluxon” |Wil+22]. 
Magnetic diffraction pattern of the critical current 


The occurrence of Josephson supercurrents is an interference effect of two macro- 
scopic wave functions that may vary both in time and space. As already mentioned 
in Sec. 2.4.1]and what is known as AC Josephson effect, applied electric fields (DC 
voltages) make for a temporally rolling phase difference, according to the second 
Josephson equation (2.60)). This yields (oscillating) supercurrents that are periodic 
in time. Applied magnetic fields, however, cause a spatially rolling gauge-invariant 
phase difference, as described by Eq. (2.71). Therefore, if the self-field is small 
compared to the external field, which is the case for short junctions with £ < Aj, 
the Josephson current across a rectangular junction of length £ and width w results 
in 


er 2rtugHA sin( 3° ) 
D(H) = J dy J dxje sin ( HA, + oo) = T.(0) te sin(@o). 
-w/2 -4/2 Po 


(2.79a) 


As can be seen in Eq. (2.79a), the Josephson current vs. magnetic field distribution 
generally equals the Fourier transform of the junction’s cross-section |Tin04]. In 
particular, the magnetic diffraction pattern of a rectangular junction cross-section 
gives rise to a sinc-shaped pattern function. Maximizing the current with respect 


to po, results in the critical current [Jos64 


I, = I,(0) (2.79b) 
which resembles the Fraunhofer pattern in optical single-slit experiments, where 
two wave functions also interfere with varying path differences. In distributed 
Josephson junctions, this path difference originates in the magnetic field depen- 
dence of the gauge-invariant phase difference and modifies the interference ampli- 
tude of the two superconducting wave functions. Figure 2.13|shows both simulated 
and experimentally acquired magnetic diffraction pattern of short junction’s critical 
currents. 
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Figure 2.13: Fraunhofer-like magnetic diffraction pattern of a short junction’s critical current. a) Mea- 
sured diffraction pattern of an 80 um x 5m junction in inline geometry with electrodes made from 
high-kinetic inductance AlO, and pure aluminum. The slight asymmetry arises because of the inhomo- 
geneous bias current distribution in the inline bias scheme in combination with different inductances 


in the bottom and top electrode |Sch70}|BJV75;|Mon+13]. These experimental data resemble b) the 
P P 


analytically calculated sinc-shaped diffraction Fraunhofer pattern, which, however, features dissonantly 
sharp turns at integer multiples of o. More advanced simulations |PRS91], assuming a normalized 
junction length ¢/Aj = 0.3, cover this detail splendidly. 


The simulated data are generated with the approach of Pagano ef al. that is 
based on evaluating the phase space configuration Px(P) of the perturbed Ferrell- 
Prange equation. In doing so, stable points can be found at a given applied field 
and for different fluxon branches, which yield extremal junction length depending 
on the applied bias current. This dependency is numerically inverted to find a 
critical current I. at a fixed junction length £. 


In case of the above analytical derivation, the self-field, generated by the induced 
screening currents is completely neglected, such that the magnetic field inside the 
tunnel barrier equals the applied magnetic flux density B = 9H. This assumption 
is valid as long as the applied field penetrates the entire junction, which is the case, 
while the junction length £ exceeds the twice Josephson length Aj, the so-called 


“short junction” limit £ S 2af] [0567]. 


In the opposite regime of long junctions ¢ > Aj, however, self-fields become sig- 
nificant since screening currents repel magnetic fields from the junction’s interior, 
in the same way as described by the Meissner-Ochsenfeld effect of bulk materials 


14 The factor 2 arises since the field penetrates the junction from both sides, as illustrated in Fig.|2.16 a) 


42 


2.4 Dynamics of long Josephson junctions 


Ly AAN l ANA l A 


345 45012 3 4 


normalized applied magnetic flux ®/ ®o 


1.0 
0.8 
0.6 


0.4 


0.2 


normalized critical current /./ I.(0) 


0.0 


Figure 2.14: Magnetic diffraction pattern of the critical current for different junction lengths. The 
simulations follow the approach of Pagano et al. and assume a homogeneous bias current 
distribution, such as is the case for the overlap geometry. The initially rounded lobes of the Fraunhofer 
diffraction pattern of short junctions (2 = 0) become more and more linear with increasing junction 
lengths. Thereby, the main maximum widens and eventually covers the former adjacent higher-order 
side maxima. The different lobes correspond to the number of inherent flux quanta, while it should be 
stressed that junctions enter the state that may carry the maximal possible current to minimize their 
free energy, which is observed in experiments, too. 


(see Sec. [2.1) Similar to type-II superconductors, in Josephson junctions, screening 
currents increase with the applied magnetic field to maintain perfect diamagnetism, 
which collapses above a first critical field Ha. There, entering of a magnetic flux 
quantum in the junction is energetically favorable, and since screening currents 
only need to repel parts of the applied flux, they decrease. As a consequence, 
the junction can carry larger bias currents, and a side maximum occurs in the 
magnetic diffraction pattern. Increasing the magnetic field further lets the same 
effect start all over again, and each lobe corresponds to the number of inherent 
magnetic flux quanta, starting from zero in the case of the main maxima [OS67]. 
The Meissner-like linear increase in screening currents manifests in long junc- 
tions’ magnetic diffraction pattern of the critical current as linear slopes of the 
individual lobes. As can be seen in Fig. there is no sharp transition from short 
to long junctions with any hard criterion, but typically junctions are denoted as 
one-dimensional “long junction” for normalized junction lengths 0 = ¢/A, = 7, 
whereas short widths w < Ay do not allow for spatial phase variations along the 


y-direction [PRS91]. 
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Aside from their normalized length, long Josephson junctions’ exact magnetic 
diffraction pattern is defined by the junction geometry as this goes along with dif- 
ferent current distribution [0S67}|Sch70} |BJV75} /BP82]. There are three widespread 
main types of linear electrode configurations, whose schematic and length- 
dependent maximal critical current are shown in Fig. 


e overlap geometry: In this configuration, as shown in Fig. |2.15a) (i); bias 


currents distribute homogeneously across the whole length, such that the 
junction’s maximum Josephson current increases proportionally to the total 
junction area as 


I 
Dt (2.802) 
2wAjjc 247 
e inline geometry: In the symmetric inline arrangement, as depicted in Fig.|2.15 a) 
bias currents are injected and ejected close to the two edges over a length 
of Ay, which is why the normalized maximum Josephson current holds 
inh( +- 
D 2 sinh (5) (2.806) 
2wajj A ’ 
He cosh (54; +1 


and thus converges towards 2 for long junctions. An asymmetric inline 


geometry, as illustrated in Fig.|2.15 a) (iii)|is also possible, where only the one 


edge close to the bias leads is stressed, and thus the normalized maximum 
Josephson current is given by 


D £ 
aA tanh 6 ) i (2.80c) 


saturating at 1 for long junctions. 


e crossover geometry: This configuration, as schematically presented in Fig.]2.15h) 
is a particular combination of inline and overlap configuration, where 
bias currents are injected along the short side (as in the case of inline ge- 
ometry) and ejected along the long side (like in overlap junctions). If the 
Josephson inductance dominates the low lead inductance, the normalized 
maximum Josephson current is 1 < lim% Beh < 2, depending on the 
fraction of bias current flowing in each edge. In the particular case of high- 
and low-inductive leads in the y- and x-direction, this yields a homogeneous 
current distribution, which is why it is referred to as quasi-overlap geometry 


[Sar+91). 
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Figure 2.15: Schematic layout and maximum Josephson current for different junction geometries. 
a) Popular one-dimensional linear long junction geometries are referred to as (i) symmetric inline, 
(ii) asymmetric inline, (iii) overlap, and (iv) crossover. b) These configurations differ in their current 
distribution, for what reason, their maximum Josephson current varies significantly. Merely the overlap 
geometry provides a uniform current distribution, and thus its maximal Josephson current does not 
saturate in the long junction limit, as otherwise. This figure is adapted from Ref. [BJV75}. 


2.4.3 The sine-Gordon model of long junctions 


To analyze phase dynamics, covering both spatial and temporal phase variations 
p(x, t), the junction capacitance needs to be considered similarly to the RCSJ model, 
but for distributed junctions with inductive electrodes M] Within this scope, one- 
dimensional junctions can be modeled by N — œ infinitesimally small lumped 
elements of resistively and capacitively shunted junctions in the z-direction, which 
are extended along the x-axis, and thus connected via inductive and resistive leads. 
The finite resistance in parallel to the electrode’s inductance is due to surface losses. 
As shown in Appendix {A.1] the inductances of both electrodes can be mapped 
to a single effective inductance, resulting in the equivalent circuit illustrated in 


Fig.|2.16 b)| Using common circuit quantization techniques |VD17] with N node 
fluxes ®;, the system holds the Lagrangian 


2 
N-1 G ®_1 Dr t oe) 


N-1 
cIwAx 2 
Per ce) Sr: (2.81a) 
N-1 n; N-1 
PojewAx 27T i : 
L aa (1 cos( z 2) L hwaâr®; 
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a) 


Figure 2.16: Schematic screening currents and equivalent circuit of long Josephson junctions. a) Ac- 
cording to the Meissner-Ochsenfeld effect, screening currents (black solid lines) are induced by an 
applied magnetic field, which penetrate the bulk electrodes on the scale of the London penetration 
depth ÀL. By comparison, superconductivity is suppressed in the junction region, for what reason 
the penetration depth there, the Josephson length Aj, exceeds the A, significantly. Screening currents 
penetrate vertically on the scale of the magnetic thickness A that composes of the tunnel barrier’s 
thickness dy and the À, of each electrode. b) The long junction’s equivalent circuit is represented by 
lumped elements of resistively (ri) and capacitively (ci) shunted junctions (Jc) in the z-direction, 


which are extended along the x-axis, and thus connected via inductive (Lj) and resistive (Rj) leads. 


£ 2 : 
_ cy 2 (0x2 = A) Poje 27T RE 
| I (9) m ae (cos g) ) + jo dr. 


(2.81b) 


Here, cy = C}*/w£ denotes the specific junction capacitance, Lg = Li w/l the 
electrode’s sheet inductance that sums up both the geometric and kinetic parts, 
je = ['/wé, and jp = Ift/wé describe the critical and bias current density, 
respectively. The quoted bias current is commonly driven by external current 


Hy (x) 
x 


c : ; d 
sources, but magnetic field induced currents jz = —j,— can be treated analogously. 


During the transition to the continuous limit of distributed long junctions (N — ©, 
Ax — 0), the sum is replaced by an integral LE a . AX > i, ...dx, the difference 
PP; _ d(x) 


quotient becomes a derivative lim, 0 ~y ~ = -q and the external flux is 
an Dt arte) 
transferred to the electromagnetic vector potential lim,,;,9 —y — = —g = 


A(x). In order to describe the junction’s phase dynamics, node fluxes are translated 
to node phases via the phase-flux relation ®; = Loj 


With the assumption of static external magnetic flux 3p®*t = 0, the phase dif- 
ference can be redefined, shifted by electromagnetic vector potential A, so that 
9 9+ oe J Adz = 9 — o™ with the externally applied phase gt. Finally, the 
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long junction’s phase dynamics are described by the sine-Gordon Lagrangian in 
normalized units 


[L= i (dp)? 5 (0x9) (1 —cos( 0) +79| dz. (2.81c) 


Here, the Josephson plasma frequency wp, the Josephson penetration depth Aj, 
and € are the characteristic inverse time, length, and energy scales in the long 
Josephson junction, such that these quantities are normalized to 


= 27Tje 
t= ht, Wp = j 2.82a 
p eee (2.82a) 
Po 
t=x/r = 2.82b 
¥=x/Ajy, J mja (2.82b) 
Ae _ DojcAyw 
Y = jo/ jc- (2.82d) 
Applying the Euler-Lagrange equation [BP82}|Leg87; CF91] yields the equation of 
motion [Jos65| 
u d ƏL | = ƏL a (2.83a) 
did; dxdzp dp 
y = 079 — IP + sing. (2.83b) 


This so-called sine-Gordon equation is a nonlinear wave equation, which describes 
dynamic properties in a Josephson transmission line. The characteristic prop- 
agation velocity inside this Josephson medium is denoted as Swihart velocity 
Tit and is typically in the order of ~10~* co. The corresponding 


junction impedance Zj = 4 y E & ©! is consequently in the order of a few ohms. 


C= WpAy = 


For the purpose of Lagrangian formalism, dissipative perturbations that are un- 
preventable in real long junctions are usually included by a Rayleigh dissipation 


functional R = f % (9)? dx [JV09}|EGK15]. Introducing both the ohmic quasi- 
particle tunneling resistance Ry = pydj/wAx = r;/wAx and the lateral electrode 


resistance due to surface losses in the superconductor Ro = ppAx/wd = Rp Ax/w, 


as depicted in Fig. |2.16b)} results in 


N-1 N-1 
wAx 2 w 2 
R= Y — (,0,)*4 Bd; — dD; 2.84 
D Foray (i) L zra P! t®j—1) (2.84a) 
J J 
-wf | er ©)? | dx (2.84b) 
2m‘ 2RD S l 
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With the typical normalization of time, length, and energy, given in Eq. (2.82), the 
dimensionless Rayleigh dissipation functional results in |MS78 


TI L f’ (kan? + B@x9)?) ar (2.840) 


with the dimensionless damping coefficients covering quasiparticle tunneling 
across the tunnel barrier « and electrode’s surface losses 6 


= = 671/2 2.85 
7 Zntjerfe) Be" en) 
és 2 
2rie (15) (2.85b) 
> l 
Docy (R ) 


The junction’s Hamiltonian, governing its total energy, can be obtained by a classical 
Legendre transformation and is commonly split in a sine-Gordon summand ASS 
and a perturbation summand HP including bias and dissipation 


ei (2.86a) 
2 
= 1 1 
HS = I E (Ip) 5 (9:9) + (1 — cos(p — Pext))| Ax, (2.86b) 
0 
dH? 2 2] 45 
0 


The Euler-Lagrange equation including dissipation by means of the Rayleigh 


functional reads |GPS14: |Nol14] 


teen) oe ee Ge) 1 OR (2.87) 
dt \ d (drp) ] dx (3 (9:p) 9P Wp dx \Ə (ðP) wp 0 (079) 
and results in the perturbed sine-Gordon equation|MS78] 

Ip — Oxep + sin p = y — RIP + PderrP. (2.88) 


In order to describe the phase dynamics of long junctions with finite length prop- 
erly, magnetic fields, influencing the spatial phase dependence (x), need to be 
analyzed. The part of externally applied fields in the y-direction Ha always affects 
p(z) according to Eq. (2.71), whereas the part Hp, induced by a bias current Ip, 


strongly depends on the junction’s geometry |MS78; \LPS83; |Ols+86]. 
sly dep J 8 y 


In case of the overlap configuration (see Fig. the bias current both in 
the top and bottom electrode is applied in the y-direction, hence, the induced 
Hp is bound to the xz-plane. Since, it is perpendicular to the junction’s area, the 
z-component does not influence the gauge-invariant phase difference at all. Due to 
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the small junction width w < Aj, the magnetic flux through the yz-plane, generated 
by the x-component of the induced magnetic field, is negligible small for phase 
variation on the scale of Aj. That is, the bias current induced magnetic fields do 
not affect the phase junction’s dynamics that is fully governed by the externally 
applied field according to Eq. 271), such that the boundary conditions of overlap 
junctions hold 


> 2710 AA] Be 


(pz + Beez) pz + Por) |; Dy 1. (2.89) 


The situation changes for the inline geometry (see Fig.|2.15 a) (ii)), where the bias 
currents flow in parallel to the junction’s long side along x and its induced magnetic 


a 


fields are in the yz-plane. As for overlap junctions, the z-component does not affect 
p(X), but the y-component contributes in the same manner as the Hay does. Inside 
the junction (0 < x < £), Hy of the bottom and top electrode cancel to zero, but 
at one junction end, the Hp above the bottom electrode and at the other end Hp 
below the top electrode contribute with different signs. Therefore, the boundary 
conditions of symmetric inline junctions are given by 


(2+ Bex)|,_5=—@, #4 awe F FARI 
2rrug AA] I, $ 

z oT := 4 — .90b 

(pz PPri) pp © 2jeApw 7 Yz (2 90 ) 


or vice versa, depending on whether both dues sum up constructively or destruc- 
tively. 


Small-amplitude plasma wave solution 


The generally nonlinear sine-Gordon equation (2.83b) of an infinite length Joseph- 
son junction can be linearized for small amplitudes  < 1 and becomes a Klein- 
Gordon equation 


Prz — PH = P, (2.91) 
whose solution are plane waves with normalized wave number k = k - Ay, frequency 


@ = w/wp, and small amplitude po < 1 


Pp(&,F) = poè 2) + 2700. (2.92) 


The propagation of these small amplitude excitations in space and time is illustrated 
in Fig. [2.17 a)| Inserting this plane wave solution of Eq. in the Klein-Gordon 
equation yields the dispersion relation @ = v1 + K (see Fig. 2.17b)), which 
resembles plasma oscillations since it features a gap in the spectrum at a small k of 
1 (of wp in SI units). This is the reason why such a linear excitation is commonly 
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Figure 2.17: Small amplitude plane wave solution of the linearized sine-Gordon equation. a) Plane 
waves, solving the linearized sine-Gordon equation for small amplitudes po = 0.1 < 1, propagate both 
in space with the wave number k = 1 and in time with the frequency © = 1. b) Their dispersion relation 
exhibits a gap for small wave numbers k and approaches a linear course for k > 1. 


called “Josephson plasma oscillation” or simply “plasmon” [OD91]| and the former 
introduced oscillation frequency in Josephson junctions is referred to as “Josephson 
plasma frequency”. 


Josephson vortex solution 


The solution of infinite length junctions, which attracted most attention in the last 
decades, is the Josephson vortex solution. Neglecting perturbations on the right 
side of Eq. (2.88), this sine-Gordon equation is solved by 


pr(X,F) = 4arctan (op E l 7 ) ) + 27tn, (2.93) 
-ö 


which is nothing but the dynamic version of Josephson vortices of Eq. 
discussed in Sec. however, traveling along the junction with velocity 6 = v/c. 
This means that the sign of the exponent of Eq. determines the direction of 
propagation, manifesting as vortex or antivortex with opposite ring current. 


These Josephson vortices behave absolutely as solitons [Ust98], which are self- 
reinforcing wave packages with permanent shape due to balanced perturbations 
while propagating a medium. This requires nonlinear effects that cancel dispersive 
effects in a manner that the wave package’s faster frequency components are 
converted into slower ones and vise versa, resulting in the dynamic equilibrium 
of the soliton. The propagation velocity is limited to the Swihart velocity č, the 
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Figure 2.18: Josephson vortex solution in long Josephson junctions. a) The 27 phase kink solution 
(red line) spans over two neighboring valleys in the tilted washboard potential of a spatially extended 
long junction. b) At relativistic velocities, the Josephson vortex reveals Lorentz contraction such that the 
27-kink becomes sharper and sharper (from dark to bright). 


speed of light inside the Josephson medium. As can be seen from the Lorentz factor 
(1-3?) -12 in Eq. (2.93), these solitons are invariant with respect to the Lorentz 
transformation 


a a. (2.94a) 
= 0) 

= E, (2.94b) 
— Od 


since both the d’Alembert operator and sin(p) in the unperturbed sine-Gordon 
equation are Lorentz invariant. As a consequence, the solitons behave like 
relativistic particles, exhibiting Lorentz contraction (see Fig. for velocities 
approaching the Swihart velocity. Furthermore, Josephson vortices’ total energy 
are given by 

- = 8 

Er = H( 9p) = Aa (2.95) 
which features a relativistic mass increase and the normalized rest mass Mr = 8 
[MS78]. The soliton propagation can be driven up to relativistic velocities by 
external forces, such as a current bias that acts via a Lorentz-Magnus force on 
the vortex’s participating charge carriers. This acceleration is accompanied by 
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dissipation that yields an equilibrium velocity ö%, which can be determined by 
equating the time derivatives of the sine-Gordon and perturbation Hamiltonians 


di 8 dö 
dt han (1 = ps? dt’ (2.96a) 
=P u n 
i - - ° any, (2.96b) 


dt = = (1— #2)? 3P (1-22)? 
and results in the first-order ordinary differential equation for ö(f) [MS78] 


T- ao (1 a) „Bo li ey (2.97) 


The situation, where the input power supplied by the external current source 


equals the dissipated power, is characterized by de 
ö 


= 0 with the power 
=y 


balance velocity ö% and yields 


_ 4 [Goo P 
: rine ts) =R 
2.4.4 Quantumness of long Josephson junctions 


From a more quantum mechanical perspective, the long junction’s Lagrangian is 


conveniently written in the form [HSB94; |SBM97] 


he 1 1 1 
Be A (9:9)? - 5 xp - A) 7 (1 —cos(@)) dx. (2.99) 
0 
The overall coefficient 
_ hw Arch L 
iy = —P = Z . 
we T? / awk Zi (2.100) 


defines the system’s energy scale with respect to its quantized eigenenergies. In 
doing so, hj can be interpreted as a normalized Planck constant in the Josephson 
medium that is a measure for the quantumness of LJJs. It is zero in the 
classical limit, whereas the larger fiy, the smaller the system’s energy scale, and thus 
the more quantum the junction is. For a typical 1 pm wide Al/AlO,/Al-junction, 
this normalized Planck constant is about hy x 2x107 ÍKI96], which means that 
energy scale of inherent excitations such as Josephson vortices exceed the quantized 
eigenenergies by 1/hj. In consequence, such a conventional Josephson medium 
and its excitations are quite classical according to the Ehrenfest theorem. 


The quantum regime can be approached by increasing the junction impedance 
by either narrowing the sample, decreasing the junction’s specific capacitance, or 
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increasing the electrode’s inductance. In previous works, the junction widths w 
were reduced down to ~ 200nm-300 nm, yielding normalized Planck constants 
in the range fiy = 7 x 1073-9 x 107° [Kem06]. The further narrowing is 
experimentally challenging since the state-of-the-art fabrication technology limits 
the feasible structure sizes. Moreover, the widths cannot be reduced arbitrarily to 
ensure well-defined phases in the superconducting electrodes by avoiding phase 


slips |Vos+21]. The electrode’s geometric inductance L, is mostly defined by 
the actual shape and dimensions and is difficult to increase further. The specific 
capacitance is bound to suitable tunnel barrier materials that also control the critical 
current density and quasiparticle dissipation so that these three parameters are 
mainly exhausted. 


However, the advent of high-kinetic inductance materials, exceeding the geometric 
inductance by orders of magnitude, is a promising, additional parameter for fy 
that relaxes the requirement on the lithographic fabrication resolution. The ap- 
proach, taken throughout this work, to increase the junction’s quantumness is the 
replacement of the junction’s bulk electrodes with thin films of very high-kinetic 
inductance superconductors. 


53 


3 Experimental methods 


The replacement of bulk long junction electrodes with a high-kinetic inductance 
superconductor is a promising approach for increasing the junction’s impedance 
and its concomitant quantumness (see Sec. 2.4.4). In the following, the experimental 
realization chosen for this work is reported. There, the choice of granular aluminum 
as a high-kinetic inductance superconductor is motivated by elaborating the main 
properties of this material. After that, details of fabrication techniques and the 
measurement framework are specified. The latter consists of a refrigerator to cool 
the sample to cryogenic temperatures with inherent low-noise DC and RF wiring 
that ensure an appropriate environment. Finally, measurement electronics for 
transport or microwave measurements outside the cryostat are sketched out. 


3.1 Fabrication techniques for high-impedance long 
junctions 


Since most superconducting parameters are material-specific, the junction’s elec- 
trode material, providing high kinetic inductance, plays a crucial role in the ex- 
perimental realization. As seen in Tab. thin films of granular aluminum oxide 
(AlO,) are promising candidates since their large kinetic inductance exceeds the 
geometric due by orders of magnitudes. This is mandatory to increase the junc- 
tion’s impedance substantially. The combination of the material’s adjustable large 
kinetic inductances with good microwave coherence of circuits made from AlO; 
[Rot+16} |Grii+18} |Grii+19} Sch+20], are the main arguments for the usage of granu- 
lar aluminum oxide for high-quality and high-impedance Josephson junctions. This 
material further turns out to be very suitable for making long Josephson junctions 
because of the good compatibility with common AlO, tunnel barriers that are most 
widely used for quantum junctions due to performance and the straightforward 
native oxidation. In addition, AlO; is resilient to magnetic in-plane fields 
(Chu+81b} /Bor+20], which, for instance, are required to inject Josephson vortices 
in long junctions. Moreover, the critical temperature Te ~ 1.6 K-2.0K well exceeds 
typical dilution refrigerators’ base temperatures of ~10mK, avoiding thermal 
quasiparticle excitations. 
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Figure 3.1: Transmission electron microscopy photograph of granular aluminum oxide. The high- 
lighted areas with parallel stripes in the diffraction pattern are identified as crystalline pure aluminum 
grains with an approximated diameter of 3nm-4nm. These nanocrystals are randomly embedded in 
an insulating amorphous aluminum oxide matrix, such that this granular material can be regarded as 
disordered network of (super-) conducting islands separated by tunnel barriers. 


3.1.1 Thin films of granular aluminum oxide 


Initiated by the superconducting isotope effect and suggestions of electron-lattice 
interactions that are affecting the superconducting behavior, superconductors have 
been studied in the dirty limit for decades. In doing so, Buckel and Hilsch found, 
pure aluminum, deposited on a 4K cold quartz substrate, to form a granular 
structure and to reveal a increased transition temperature up to Te = 2.59 K |BH54]. 
The grain size can be reduced, while increasing the critical temperature to Te ~ 3K 
by depositing aluminum inside an oxygen atmosphere [ACC66]. 


Such granular aluminum oxide is used throughout this work and consists of 
pure crystalline aluminum grains embedded in insulating amorphous aluminum- 
oxide, as depicted in Fig. Both the normal and the superconducting transport 
properties are dominated by these intrinsic nanoscopic tunnel barriers [Deu+73a]. 
In the superconducting phase, this granular structure can be treated as a network 


of nanoscopic Josephson junctions [Deu+73a]| providing a kinetic-type Josephson 
inductance [Jos64], which is related to the normal state tunnel resistance Rn and 


the superconducting gap A by Ly = hRn/ nA |MB58}|GT57 (see 


Eq. (2.32)). The AlO;’s conductivity and inductance can range over five orders of 


magnitude [Deu+73b; |Chu+8la] depending on the oxygen concentration in the 
nanoscopic tunnel barriers (TBs) |ZHB78], which is controlled by the oxygen partial 
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pressure during the reactive sputtering process [Rot+16]. This enormous versatility 
enables us to use AlO, in different regimes for depositing junction stacks with 
high kinetic inductance as bottom electrode, an insulating tunnel barrier, and a top 
electrode made from AlO, and pure aluminum. 


3.1.2 Trilayer sputter deposition 


The junction trilayers are deposited on a c-plane sapphire substrate using a reactive 
DC magnetron sputter process in a homemade sputter tool. This device consists of 
an ultrahigh vacuum chamber to avoid contamination and a sputter gun formed by 
a ring-shaped anode close to the aluminum sputter target that serves as the cathode. 
Applying sufficiently high voltages ~ 400 V enables a dielectric breakdown of the 
argon working gas, ignites an argon plasma that bombards the sputter target, and 
thus eject clusters of aluminum atoms. The sputtering rate is increased via a radial 
magnetic field that deflects charge carriers on an elongated spiral track, leading to 


an intensified plasma |DMO8]. 


To grow AlO,, oxygen is added as a reactive gas that oxidizes aluminum naturally 
and yields the desired granularity. The oxygen concentration, thus the normal 


sheet resistance Rg, and the concomitant kinetic inductance per unit square Lg, are 
defined by the oxygen partial pressure during the deposition [Rot+16]. Especially 
for sheet resistances beyond ~ 1 KO, the resistivity is so sensitive to slight variations 
of the oxygen partial pressure that the reproducibility suffers from this exponential 
dependence. 


To gain necessary control in this regime, both film thickness d and resistance Ry 
are monitored during the film deposition]|| The film thickness is measured using 
a quartz oscillator, whose resonant frequency shifts with the mass of deposited 
material. This quartz oscillator is placed close to the sample, and the frequency 
shift to thickness conversion is calibrated using a profilometer. 


Conventional four-point probes of the film resistance are falsified by large currents 
~10mA that flow from the plasma to the sample. To overcome this impact, large 
measurement signals, exceeding this plasma current by orders of magnitudes, 
cannot be used since they heat the film beyond its melting point. Such heating 
results in an inhomogeneous resistivity with undefined granularity and even 
bubbles on the film surface. Therefore, a counter voltage ~ 200V between the 


1 This in situ R(d) measurement enables the estimation of the bottom electrode’s sheet inductance, 
which would not be possible for the whole trilayer. 
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insulated sample holder and the grounded chamber is applied to decelerate the 
argon plasma, and thus avoid targeting the film and falsifying the measurement. 
The uncharged sputter material is unaffected by this electrical field, and the sputter 
rate remains the same. 


This in situ R(d) measurement enables us to fit the specific conductance o(d) = 


= (Raa) T and to estimate the final specific resistance pno already during 
the deposition [Wil+22]. The granular aluminum oxide’s thickness-dependent 
specific conductance follows the model for fine-grained polycrystalline thin films 


by Mayadas et al. |MSJ69] 


0 _ Pn0 1 a 2 3 1 
= = H In(1+-=-)|. a 
= an 3 377 a“ — a n( +2) (3.1) 


Equation is the solution of a linearized Boltzmann equation including ordinary 
scattering mechanisms as in bulk materials and superimposed scattering at grain 
boundaries. There, vo is the intrinsic thickness-independent conductivity from the 
film interior and a := ee denotes the ratio between the background mean free 
path Io and the film thickness d, reduced by a scattering reflection coefficient r. 


For the in situ fit, as shown in Fig. an offset thickness has to be included, 
above which the model holds. The fit performed for trilayer A yields the resistivity 
po = (70.7 + 0.2) Q cm and zur = (5.04 + 0.17) A. The product of resistivity and 
mean free path is a material specific measure of its Fermi surface Sp and can be 


assumed to be polo = 127°h/e*Sp = 1.6 x 1071! Q cm~? for granular aluminum 
oxide [CA68]. This yields the mean free path Io = (22.62 + 0.07) A in the order of a 
few nanometers and the substantial reflection coefficient r = (18.2 + 0.5) %, which 


confirms diffusive transport in the granular material. The discrepancy between 
measurement and theoretical model for small film thicknesses originates in an 
inhomogeneous film thickness and the conductivity of the not perfectly repelled 
argon plasma that contributes especially for small film thicknesses, where the film 
is not entirely connected. 


As can be seen in Fig. 3.2] this in situ R(d) measurement is a powerful tool to achieve 
the aimed kinetic inductance value. On the one hand, the deposition can be stopped 
when reaching the desired sheet resistance, subjecting the exact film thickness to 
slight variations. On the other hand, if the fitted final resistivity deviates from the 
goal, the oxygen partial pressure can be readjusted, visible as knees in the R(d) 
curves for trilayers B and C in Fig. Consequently, the desired sheet resistance 
RẸ at the fixed film thickness can be achieved with a deviation of less than 10 % 
[Wil+22], which is great progress in controlling the sheet resistance of granular 
aluminum oxide. 
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Figure 3.2: In situ resistance measurement during reactive sputter deposition of different trilayers. 
a) On two opposing sides of the substrate, there are previously deposited silver strips to ensure that 
the resistance of the squared film equals the sheet resistance, which is connected to a four-point 
measurement. b) The measured thickness-dependent normal resistance Rn (d) is fitted to the model of 
fine-grained polycrystalline thin films by Mayadas et al. in order to estimate the final resistance. 
If needed, the oxygen partial pressure is readjusted during the deposition to achieve the desired 
resistance, as is the case for trilayers B and C. The inset points out the tunnel barrier creation, where 
low-ohmic pure aluminum is deposited (a to b) and afterward statically oxidized (b to c). (This Figure 
is also published in Ref. .) 


The trilayers’ tunnel barriers are formed by depositing a thin film of pure alu- 
minum, that is later oxidized statically. There, the critical current density is con- 
trolled by the product of oxidation pressure pox and time tox [KMM95 


jc(Pox, tox) = Aox (Pox : tox) : (3.2) 


Fitting results of the later transport measurements yield the constants aox = (3.36 + 
0.12) x 10* and box = —0.666 + 0.016, when using the SI units ampere, meter, bar 
and second. For details, see Fig. [B.2| given in Appendix|B.2| 


The barrier thickness can be estimated from R(d) in situ measurement. The pure 
aluminum layer as a low-ohmic shunt first decreases the total sheet resistance (from 
a to b in the inset of Fig. before coming back to the initial resistance at the 
end of the static oxidation of the whole aluminum layer (from b to c in the inset of 
Fig. [3.2). A complete oxidation is achieved for ~ 1.2nm-1.5nm thick aluminum 
layers, but slightly depends on the barrier’s transparency. 


Finally, a top electrode is deposited similar to the bottom electrode but eventually 
with different oxygen concentrations. The characteristic parameters of the trilayers 
used in this work for long junction experiments are given in Tab. 
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Table 3.1: Trilayer fabrication parameters. During the deposition of the bottom electrode, the sheet 
resistance R is measured by four-point probes, and the film thickness dı is measured with a quartz 
oscillator. To create the tunnel barrier, a layer of pure aluminum of thickness dy is first sputtered and 
statically oxidized in an oxygen partial pressure pox for a time tox. The layer of pure aluminum is then 
deposited on the bottom electrode. If necessary, the top electrode of thickness dz is continued with 
another layer of thickness d3. 


trilayer Ra Pox tox dı dy da d3 
(Q) (bar) (min:s) (nm) (nm) (nm) (nm) 


38 21.0 20:00 20.4 2.6 20.6 20.0 

75 051 4:00 20.4 2.1 21.1 20.2 

778 51 4:00 18.6 1.5 20.4 20.1 
2295 0.025 4:00 34.6 1.2 30.0 
947 0.286 4:00 46.6 1.2 30.0 


HoOnuwu» 


Long Josephson junctions are patterned from these stacks using e-beam and pho- 
tolithography techniques in combination with inductively coupled plasma etching, 
anodic passivation and thermal evaporation of leads to the junction’s top electrode. 
As the exact process differs for junctions read out in the DC or the RF scheme, this 


is discussed in detail in Sec. [4.2} Sec. 5.2] and Sec. 


3.2 Measurement setup 


To study the physics of long Josephson junction devices, an appropriate framework 
of low-noise DC and RF connections as well as proper static magnetic fields are 
required. Therefore, the long junction samples are hosted in a gold-plated oxygen- 
free copper sample holder that provides two coaxial SMP jacks, which are soldered 
to 50 Q matched microstrip transmission lines on microwave-suitable printed circuit 
boards (PCB). Aluminum wire bonds both connect these PCBs with the on-chip 
feedline and the DC microstructures with another PCB that supplies 24 DC lines 
from a micro D metal (MDM) connector. Moreover, a homemade superconducting 
niobium coil (Te 7 9.2 K) is placed around the sample holder to apply an in-plane 
magnetic field. To shield external magnetic fields, the whole box is surrounded 
by a Cryoperm shield [] The assembly used throughout this work is depicted in 


Fig. 


2 Cryoperm is low-temperature suitable permalloy made from a soft magnetic nickel-copper- 
molybdenum-iron alloy. Due to its high permeability u, ~ 70000 at 4K and the entailing low 
reluctance, magnetic flux prefers being inside the shield instead of in the vacuum at its inner. 
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Figure 3.3: Sample holder with DC and RF connectors, surrounded by a magnetic coil and a permalloy 
shield. a) The sample chip, containing long junction experiments, is embedded in a sample holder, 
where aluminum wire bonds contact both junction’s leads and microwave feedline to appropriate 
PCBs installed in the sample holder. b) The DC lines are connected to an MDM connector, whereas 
the RF-suitable PCBs are soldered to SMP jacks (accessible from the sample holder’s bottom side, and 
thus not visible on this photograph). The hosting sample holder and its lid are made from copper and 
gold-plated to provide a good thermal connection to the cryostat. c) This closed box is surrounded by 
a superconducting coil to apply well-defined magnetic fields and mounted deep inside a Cryoperm 
shield, which screens magnetic fields from the sample. 


Naturally, Josephson junction devices, made from aluminum (oxide) trilayers, ne- 
cessitate cooling down to millikelvin temperatures. First superconducting phenom- 
ena, such as zero DC resistance, ideal diamagnetism and even coherent Josephson 
tunneling, already occur just below Te. However, avoiding thermally populated 
dissipative quasiparticles and thermally activated tunnel processes that would 
destroy quantum coherence require way colder temperatures. 


3He/*He dilution cryostats meet the demanding requirements and provide contin- 
uously sufficient cooling power to cool the solid bodies of samples and necessary 
measurement equipment. Therefore, the whole sample holder as portrayed in 
Fig. B.3lis thermally connected to the base plate of a ?He/*He dilution cryostat. 


The following describes how such refrigerators operate in principle and how they 
must be equipped to suit measurements on Josephson junction devices. In this 
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work, these samples are probed both by DC and RF measurement techniques, 
whose essentials are presented after that. 


3.2.1 Pulse tube precooled *He/*He dilution cryostat 


3He/*He dilution cryostats are multistage cooling units whose precooling to a 
few kelvins can either be implemented by baths of liquid nitrogen and helium or 
mechanical cryocoolers like a closed cycle pulse tube refrigerator. The latter are 
very popular as they circumvent moving parts in their low-temperature part, and 
thus avoid vibrations. As shown in Fig. [3.4| two of such pulse tube refrigerators 
are thermally connected to the first two stages and precool the ?He/?He mixture 
down to liquefaction. This is cooled further by means of a Joule-Thomson valve, 
where the temperature decreases under isenthalpic expansion [ADK18]. 


Now, below ~ 870 mK, the 7He/#He mixture undergoes a spontaneous phase 
transition and separates in two phases: a He rich phase (concentrated phase) 
that consists almost exclusively of ?He and a ?He poor phase (dilute phase) that 
consists in its thermodynamic equilibrium at 20 mK of ~ 93.5% *He and © 6.5% 
3He [EH05]. Due to their different densities, the concentrated phase in the mixing 
chamber floats above the dilute phase, which is connected with a reservoir at still. 
There, the higher temperature and the decreased pressure favor the evaporation 
of almost exclusively “He because of the isotopes’ different vapor pressures. As 
a result, there is a nonequilibrium with a concentration gradient that goes along 
with an osmotic pressure between the dilute phases at the base and still stage. 
To compensate for this and maintain the thermodynamic equilibrium, *He of the 
concentrated phase must mix itself with the dilute phase. The required entropy 
of mixing is extracted from the system in the form of heat [EH05]. To ensure 
continuous cooling, the evaporated He is sent to the cryostat again, cooled by 
heat exchangers on the condenser side, and fed into the concentrated phase in the 
mixing chamber. 


Even though this cooling principle enables theoretically arbitrary low temperatures, 
its finite cooling power ~ 10 pW limits the base temperature of 7He/*He dilution 
cryostats to ~ 10 mK, because of thermal connections between the stages, heat 
input through circulating mixture, thermal radiation, and unavoidable Joule heat 
input of electronic devices. Nontransparent metallic shields at each stage suppress 
the latter, where the warmest ones have the greatest bearing according to the 
Stefan-Boltzmann law P « T*. For this reason and for convenience, the used 
Bluefors LD250 cryostat abstains from a base shield. Thermal connections between 
the stages as well as their thermal contact with the environment are reduced by 
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Figure 3.4: Cryogenic measurement setup. a) The cryostat is a multistep refrigerator, whose precooling 
unit is a two-stage pulse tube generating about 40K and 3K, respectively, that is thermally coupled 
by heat exchangers to the primary dilution refrigerator. By means of a Joule-Thomson valve in this 
3He/*He cycle, temperatures at still of about 700 mK are reached, which suffices to obtain a separation 
into a concentrated and a dilute phase of the 7He/*He and run the dilution refrigeration cycle. Thereby, 
3He is pumped down at still from the dilute phase and recondensed to ensure continuous cooling. 
b) The photograph of a Bluefors LD250 refrigerator shows the implementation used in this work, 
including the wiring. c) The measurement framework consists of low-pass filtered DC wires, coaxial RF 
lines, and high-current supply lines for magnetic coils. 
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pumping the cryostat’s outer vacuum chamber to ~ 1076 mbar and by an elaborate 
choice of materials used for the cooling unit and measurement wiring. 


The later presented experiments were run in the commercial ?He/*He dilution 
cryostat Bluefors LD250, as depicted by the photograph in Fig. [3.4b)| whose mea- 
surement setup is configured to the special needs of Josephson devices. For this 
purpose, DC and RF connections, as well as high-current supply lines, are mounted, 
where the heat input to the individual stages due to thermal conductivity is kept 
minimal by the careful choice of suitable materials and wire cross-sections. Great 
care is also taken to reduce the thermal noise in the measurement lines by means 
of filters at different temperatures. 


Commercial microwave electronics generate signals at room temperatures with 
order of magnitudes too high powers for Josephson devices that require a few to a 
few thousands of photons. This is the reason why thin stainless steel coaxial cables 
with comparable large damping (23.5 dB m! at 300K (compare to Fig. in 
Appendix[B.1) and 14.7dBm 1! at 4K, both at 1 GHz) can be used as input lines, 
to reduce the thermal conductance between the individual stages. The stainless 
steel’s comparably high electrical resistivity originates in a strong electron-phonon 
interaction that is beneficial to thermalize the electrons. For experiments with a 
few photons, the input signal’s noise temperature, which is approximately room 
temperature, is of particular interest and needs to be reduced. Therefore, it is 
advantageous to thermalize attenuators at the lowest possible stages that provide 
sufficient cooling power since each resistor contributes with thermal Johnson- 
Nyquist noise. The signal’s final noise temperature passing an attenuation chain of 
N stages, each with attenuation factor aj > 1 and temperature Tj (with Ti ~ 10mK 


at base and room temperature Ty+1 = 300 K) results in [Fri44] [Sch20] 


N T 
T=T™+)—*. (3.3) 
i=1 l lj=1 fj 
As the base temperature T; contributes without any attenuation, this should be 
kept as low as possible, and due to the limited cooling power, the attenuation 
cannot be arbitrarily large there. For this reason, and to decouple the cold stages 
thermally, attenuators are typically distributed over the three coldest stages, as 
illustrated in Fig. 3.4 c)|[Kri+19]. The total attenuation, which is usually desired to 
be in the order of ~ 70 dB, sums up by about 40 dB at 10 GHz in the stainless steel 


cables and of several attenuators with in total 30 dB in the present setup. 


To analyze the signal’s few photons, scattered by the sample, at room temperature, 
resistive coaxial lines are not applicable, by which ~ 99.99 % of the signal would 
be lost. This requires amplification at cold temperatures that is realized by high 
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electron mobility transistors (HEMTs) at the two pulse tube stages|>| The sample 
is optimally connected with the first HEMT by superconducting niobium coaxial 
lines due to their perfect electric conductance and the coexistent thermal isolation. 
After an amplification of ~ 35 dB, normal conducting coaxial lines can further be 
used. In the present setup, the signal is led through silver-plated stainless steel 
cables that combine weak thermal conductivity in the outer and the middle of 
the center conductor and good electric conductivity underneath the surface of the 
center conductors. This silver layer is exactly the area where most of the microwave 
signals are located according to the skin effect. 


In close vicinity to the sample, both input and output signals are low-pass filtered 
with a cutoff frequency of 18 GHz and infrared filters based on Mie scattering on 
sapphire powders of different diameters. The latter filter is essential as the coaxial 
line’s and attenuator’s typical dielectric PTFE, commonly known by its brand name 
Teflon, is transparent to infrared radiation. Moreover, to avoid incoming (reflected) 
signals at the sample’s second port, a circulator with a 50 Q termination is installed 
that only allows for transmission in one direction and otherwise attenuates by 


~ 20 dB (see Fig. |3.4c)). 


Conventional source measure units for transport measurements provide insufficient 
signal-to-noise ratios, so these signals must also be filtered. In doing so, the DC 
lines are made from twisted pairs of constantan wires with a rather high resistance 
of 66 Q m71, all embedded in epoxy glue with copper powder acting as infrared 
filter [LUO8], and shielded by a stainless steel mesh. The signals are further low- 
pass filtered at room temperature by a CLC filter in 7 topology with a cutoff 
frequency of about ~ 300 kHz, and at the second pulse tube stage by an RCR filter 
in T topology with a cutoff frequency of about = 1kHz. The latter RCR filter 
is installed at ~ 3K in order to reduce the resistors’ thermal Johnson-Nyquist 
noise voltage Vn = J 4kgTRAf, while having sufficient cooling power available. 
More details such as the transmission spectrum of such DC measurement lines are 
specified in Fig. [B.1b)|in Appendix {B.1] 

Supply lines for magnetic coils that need to carry high currents consist of massive 
copper wires between room temperature and the first pulse tube stage, whose 
thermal load is increased in this manner. Due to its good thermal conductance 
and the lower cooling power at further stages, such copper wires with large cross- 
sections are not possible to use there while maintaining millikelvin temperatures 


3 A semiconductor heterostructure provides a two-dimensional electron gas with high electron mobility, 
such that these amplifiers can operate at a few kelvins, where most other semiconducting amplifiers 
would not be able to operate. 
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at base. This is why high-temperature superconductors, featuring zero resistance 
and good thermal isolation, are used between both pulse tube stages. Below this 
stage, with about ~ 3K, conventional metallic superconductors are available that 
are way cheaper and much easier to handle. For this purpose, wires are used where 
multiple filaments of NbTi are embedded in a CuNi matrix. With this combination 
of materials, avoiding plug-connections and using only clamp connections and 
superconducting solder joints, the wire resistance is minimized to about 0.3.0, 
while enabling coil currents of a few amperes with negligible Joule heat input to 
the cryostat. 


As ?He/*He dilution refrigerators are conceptualized to operate at millikelvin 
temperatures, but not around 1K, an additional stage is thermally weakly coupled 
to the mixing chamber by thin-walled stainless steel tubes if needed for such tem- 
peratures. This extra stage has an individual heater and a separate resistive RuO2 
thermometer, thus allowing for precise temperature control and fast temperature 
sweeps. This gadget benefits from its small thermal mass and weak thermal con- 
nection, so it avoids heating the mixing chamber beyond the critical temperature, 
where the mixture would vaporize. Thus, temperatures around 5K can be achieved 
easily, with base not exceeding 100 mK. 


3.2.2 Transport measurement setup 


To take IV-characteristics of current biased Josephson junctions with critical cur- 
rents down to a few nanoamperes, high signal-to-noise ratios are required even at 
current resolutions in the order of tens of picoamperes. The bias current is gener- 
ated by means of a homemade voltage biased current source, providing different 
amplification factors (2 x 107? A V7! to 2 x 1076 A V-!). Its core is a reference 
amplifier consisting of two high-precision monolithic dielectrically isolated FETs 
OPA2111. This current source is supplied by a source-measure unit (SMU), whose 
output signal is prefiltered by the low-noise preamplifier Stanford SR560. 


The measured voltage drop across the Josephson junction, in turn, is amplified 
by a low-noise instrumentation amplifier INA 105KP and measured by the SMU. 
IV-characteristics are recorded by linear staircase sweeps, provided by the SMU, 
which is in this work either the sourcemeter Keithley 2636A or the data acquisition 
system NIDAQ USB-6259 BNC, if the real-time resolution is needed. 


A method to take data of transport measurements with customized sequences of 
sweeps is integrated into the Qkit measurement framework [gki] that is developed 
in Alexey Ustinov’s group at KIT. 
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3.2.3 Microwave measurement setup 
Spectroscopic measurements 


Scatter parameters of 2 port networks can be acquired during frequency sweeps 
for each frequency point employing vector network analyzers (VNA) such as the 
model Keysight VNA E5080B, which was used in this work. In doing so, the complex 
scattered signal, detected at port 2, is related to the complex signal applied at port 1 
to obey the definition of the transmission coefficient S21 in Eq. (2.33). The obtained 
microwave spectrum contains the whole cryogenic RF setup in general but is 
dominated by the Josephson device in case of proper wiring, as described above. 
As the receiver’s (at port 2) typical intermediate frequency (IF) bandwidth is in the 
order of 1Hz-10* Hz, such measurements analyze the system in the steady-state, 
which is sufficient for basic characterizations such as resonances. 


As theoretically treated in Sec. coupled systems are typically characterized via 
two-tone spectroscopy. Thereby, the second system is driven by an additional tone 
with power Py and frequency fg that is swept and in the vicinity of the resonant 
frequency, the second system is populated. This results in a dispersive shift of 
the first resonator monitored by a VNA. The second tone is generated by a local 
oscillator (LO), such as the model Keysight N5173B in this work, and is added to the 
VNA’s probe signal by means of a directional coupler. The setup is schematically 


shown in Fig. 


Time resolved measurements 


If steady-state measurements are insufficient and time-resolved measurements 
are needed, a time-domain setup, as schematically shown in Fig. 3.5 b)|is used. 
Commercially available electronics, such as arbitrary waveform generators (AWGs) 
or digital-to-analog converters (DACs), usually operate with a few gigasamples 
per second, which is why the required pulses cannot be generated directly with 
carrier frequencies in the gigahertz range. Instead, only the pulses’ envelopes are 
generated with such devices at a maximum intermediate frequency of 500 MHz 
and mixed with another carrier frequency generated with an LO, which is at most 
20 GHz. 


To obtain only a single sideband with this up-conversion, two 90° phase-shifted 
envelopes are generated and multiplied together with the carrier frequency by IQ 
mixers to form a complex signal. This way, both manipulation and readout pulses 
are generated, which are then added with a directional coupler before being sent 
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a) 21d tone spectroscopy b) manipulation readout 
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Figure 3.5: Microwave setup for steady-state two-tone spectroscopy and time-resolved measurements. 
a) In a two-tone spectroscopy microwave setup, a probe tone of a vector network analyzer (VNA) 
with frequency fyna and power Pyna is added by means of a directional coupler to a drive tone of 
another local oscillator with frequency fa and power Py. The scattered signal is detected at the VNA’s 
receiver and related to the probe signal to determine the scatter parameter S21. b) For time-resolved 
measurements, the pulses’ envelopes of both I and Q are generated by a digital-to-analog converter 
(DAC) and multiplied by IQ mixers with a carrier frequency of a local oscillator (LO) for single-sideband 
mixing. Both manipulation and readout pulses are amplified with variable gain, and their signals are 
added by a directional coupler. For the readout signal, the microwave is previously split into equal 
parts for the up and down conversion to ensure a common reference for both frequency and phase. 
The signal scattered by the sample is down-converted, low-pass filtered, amplified and recorded using 
an analog-to-digital converter (ADC). This schematic of the time domain setup is reduced to the basic 
components. For details, see the thesis of my colleague Richard Gebauer [Geb22], in whose context this 
setup was developed. 


to the sample. The pulses scattered by the sample are again down-converted with 
IQ mixers using the same LO signal as for the previous up-conversion to ensure a 
common frequency and phase reference. Before these I and Q measurement signals 
are recorded with analog-to-digital converters (ADCs), they are low-pass filtered 
to get rid of high-frequency mixing products and amplified to match the ADC’s 
input levels. 


During this work, a time-domain setup was used, which was developed in cooper- 
ation with the Institut fiir Prozessdatenverarbeitung und Elektronik (IPE) at the 
KIT. The setup combines central processing units (CPUs), field-programmable gate 
arrays (FPGAs), DACs and ADCs to generate and detect the signals 
Geb+20b]. An in-depth description of the system can be found in [Geb22]. 
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4 Fluxons in high-impedance long 
Josephson junctions 


Long Josephson junctions (LJJs), whose usual bulk electrodes are replaced by thin 
films of granular aluminum oxide, must be tested experimentally for their suitability 
for quantum circuits. For this purpose, we experimentally check whether such high- 
impedance junctions are still governed by the sine-Gordon model, which means 
that vortices can be excited in this Josephson medium. The existence of vortices 
in the LJJs is mandatory, which means that the junctions are describable by the 
sine-Gordon model. Moreover, it needs to be checked whether the additional kinetic 
inductance can really increase the junction impedance, and thus the quantumness. 
Furthermore, high-quality tunneling barriers are required to observe quantum 
coherence, which entails as few dissipative quasiparticle excitations as possible. 
All these properties are studied in transport measurements by analyzing the LJJs’ 
IV characteristics, their magnetic field dependence, their temperature dependence, 
and their response to microwave irradiation. The main results of the experiments 
presented in this chapter have already been published in Ref. [Wil+22]. 


4.1 Sample considerations and design 


In order to study the influence of the kinetic inductance in the junctions’ electrodes 
on the characteristic junction parameters, namely magnetic thickness A, Josephson 
length Aj, plasma frequency Wp and Swihart velocity ¢, the kinetic inductance of the 
bottom electrode is varied. For this purpose, the partial oxygen pressure, and thus 
the sheet resistance, is varied during trilayer deposition and the critical current 
density is adjusted during static oxidation of the tunnel barrier so that the theo- 
retically predicted Josephson length is a few tens of micrometers. In this manner, 
junctions can be patterned from the trilayers A, B, and C by photolithography with 
the same optical masks, varying the width w = 2ym-5 pm < Aj and the length 
£ = 20 um-120 pm > Aj (for stack properties see Tab.|3.1). 


Moreover, since the spatial distribution of the phase difference and current density 
is defined by the junction geometry (see Sec. [2.4.2), three different types of them 
are studied: In the quasi-overlap geometry, the high-kinetic inductance bottom 
electrode is contacted along the entire length of the junction (along x) on both sides 
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Figure 4.1: Schematic patterning of junction stacks for transport measurements. The upper row 
illustrates the structure’s top view and the lower row its side view along the dashed line after each 
fabrication step. Starting from the deposited junction stack (i) the areas of junctions and DC leads to the 
bottom electrode are patterned by an argon-chlorine-based inductively coupled plasma (ICP) etching 
process. The junctions are defined from the remaining stacks (ii) by removing the top electrode in the 
vicinity of the active area by ICP etching. The bare surfaces are further anodic oxidized to insulate the 
bottom electrode. Onto the resulting active (iii) and close parts of the passive region, DC leads are 
finally contacted to the top electrode by thermal evaporation of pure aluminum, where the galvanic 
contact is ensured by the previous removal of the native oxide at the surface by argon milling. The 
quasi-overlap geometry of the final long junction (iv), illustrated here, is characterized by two DC leads 
to both the bottom and top electrode that allows four-point probes. 


with leads that themselves consist only of the AlO, bottom layer. Their high kinetic 
inductance compared to the Josephson inductance of the LJJ leads to a homoge- 
neous current distribution since the top electrode is shunted with low-inductive 
pure aluminum, which is led out at the two short sides of the junctions (along 
y). This geometry allows four-point probes, even for room temperature prechar- 
acterizations. In the inline geometry, bottom and top electrodes are contacted at 
both ends of the LJJs to allow four-point probes. In this geometry, the bias current 
typically peaks at the junction ends, and the top electrode’s low-inductance shunt 
increases the current distribution’s resulting inhomogeneity since it also becomes 
asymmetric. In addition, squared short junctions of different sizes are fabricated 
with a current distribution as homogeneous as possible to determine the critical 
current density. 


4.2 Sample fabrication details 


Starting from trilayers A, B, and C that are deposited as described in Sec. long 
junctions that provide galvanically connected leads for probing the sample with 
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microwave 
coupler 


Figure 4.2: Micrograph of a 100 um x 5 ym long Josephson junction device in quasi-overlap geometry. 
To enable four-probe measurements at the junction, both its bottom and top electrode are galvanically 
connected with two bond pads each. One line offers an L-shaped coupler close to a microwave feedline 
that allows for RF irradiation as required for measuring the plasma frequency. In the inset, the bright 
field image is superimposed by a dark field micrograph to pronounce edges and highlight the junction 
underneath the top electrode. 


transport measurements, are patterned in three steps: First, the Josephson junctions 
and their leads to the bottom electrodes are patterned using optical lithography and 
an argon-chlorine-based inductively coupled plasma (ICP) reactive ion etching (RIE) 
process. The latter technique offers high etching rates, high selectivity and low- 
damage processing by combining a high chlorine radical density for the aluminum’s 
chemical etching and a high argon-ion density for physical removal, similar to the 
sputtering process. In doing so, all parts that are not protected by the photoresist 
are etched down to the substrate, as illustrated in the top view of Fig. 


The junctions’ exact dimensions are defined in the second step, where a photoresist 
covers the active regions and distant parts of the bottom electrode. At the bare areas, 
only the trilayers’ top electrodes are removed by ICP etching, and the resulting 
surfaces are oxidized using an electrolyte of ammonium pentaborate and ethylene 
glycol for anodic oxidatior]"| (see Fig. 4.1 (iii). This passivation is mandatory for the 
subsequent contacting of the top electrodes without shorting the junction stacks by 
the deposited superconducting layers in the third step. 


1 For the anodic oxidation, all junctions have to be connected to a common pad, where the voltage is 
applied. These joints between individual junctions are cut during dicing the chip into 5mm x 5mm 
pieces that are installed in the sample holder. 
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Table 4.1: Properties of the fabricated high-impedance trilayers. The normal conducting sheet resis- 
tance of the bottom electrode RẸ is extracted from the film deposition. The critical current densities je 
are determined from squared junctions, and Ambegaokar-Baratoff estimations coincide to switching 
current measurements. The superconducting gap energies are fitted to temperature-dependent gap 
voltages. Magnetic thickness AP and Josephson penetration depth A; are derived from the critical 
current’s magnetic field dependence of inline junctions at T ~ 300 mK, from which the geometry factor 
g(r) also follows. From overlap junctions, the Swihart velocity ¢ and the impedance Z is acquired 
from the periodicity of zero-field and Fiske steps at around 1K. (These data are also published in 


Ref. [Wil+22].) 


trilayer Ra Je At Ay 
(0)  (Acm?) (neV) (neV) 
A 38 0.32+0.03 227.1+0.8 278.6 + 0.6 
B 75 125 +03 1932+04 267.5 0.6 
C 778 1.90 0.01 201.5+0.7 278.1 0.7 
APP Ay g(r) Czts/Co Crs /Co Z 
(nm) (um) (107?) (107°) (107°) (0) 


A 69+2 — LJJ limit not reached for @<120pm 6.56 +0.03 2.79 0.04 
B 72+9 17.242.2 15+4 4.27 + 0.06 3.64+0.02 4.11 0.06 
C 9449 1951.8 1343 3.37 0.08 3.22 +0.03 14.0 0.4 


To ensure galvanic contact in the active region, the native oxides at the surfaces 
of the trilayers’ former deposited top electrodes are removed by argon milling, 
whereas the passive region remains insulating. Without breaking the vacuum, 
the final top electrode and efferent leads are deposited thermally immediately 
after (see Fig. [4.1 (iv)), which is experimentally much easier to realize with pure 
aluminum instead of granular aluminum oxide. See Appendix |[B.2| for details of 


parameters concerning optical lithography (Tabs. and B-4), ICP etching 
(Tab.|B.5), and anodic oxidation (Tab.|B.6). 


Figure |4.2|shows an optical microscopy image of a typical fabrication results of 
this process. The sample includes an LJJ in quasi-overlap geometry (see inset) and 
four DC lines, where one of them is attached to an L-shaped radio frequency (RF) 
coupler in the vicinity of a microwave feedline. 
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4.3 IV-characteristics of high-impedance Josephson 
junctions 


First, IV-characteristics of high-impedance Josephson junctions of different sizes 
and geometries are taken at low temperatures of about 20mK. As can be seen 
exemplarily in Fig. such junctions exhibit a large hysteresis between the 
switching and retrapping current that confirms a high quality factor originating 
in only little quasiparticle dissipation. The switching currents themselves are 
stochastic, which is why their distribution is measured statistically, as depicted 
in the inset of Fig. The probability that the junction escapes from the zero- 
voltage to the nonzero-voltage state within the interval dI around the current I is 


given by 
f -j pes(f) a) (4.1a) 


pe (DdI = T(t ) = 
This probability depends on the inverse bias sweeping rate dt/dl,, the activation 
rate (I) and the probability that the junction was not already switched at a lower 


dhl 


current. The integral equation for the escape probability, given in Eq. (4.1a), can be 
solved for the escape probability density (epd) [FD74a 


=1 -1 
ae (- oe / rena) (4.1b) 
0 


ecs — pecs web 
pe(1) = ree(1)|F 
and further for the escape rate |FD74a] 


fed 

ecs = dh I 

T°®(T) Ji In | — (4.1c) 
Í pess(T)dI 
I+AI 


In the experiment, the bias current is swept with a constant rate dI„/dt using 
a 16-bit digital to analog converter (DAC) and the switching current from the 
zero-voltage to the nonzero-voltage state I. is recorded for 10000 events with a 
sample rate of 1.25MS/s by a 16-bit analog to digital converter (ADC) of the data 
acquisition system NIDAQ USB-6259 BNC. A histogram is calculated from these 
switching currents, by counting the number of switching events noS within the 
je bin of the width AI, as done in Fig. The escape probability density is 
obtained by normalizing the individual number of events n®“® to unity both by the 
total number of switching events N and the bin width AI, so that 
nes 


pss(1) me (4.2) 
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Figure 4.3: Critical currents for different junction sizes and geometries, measured with various 
techniques. a) The IV-characteristic of a 10 um x 101m junction of trilayer B shows a large hysteresis 
between critical and retrapping current, indicating low damping, and thus a high junction quality factor. 
The inset points out the statistics of measured switching currents that differ every time, ramping up the 
current as each escape from the zero voltage to the nonzero voltage state is a stochastic process. (This 
figure is also published in Ref. [wil+22].) b) Different methods to determine the critical current are 
compared exemplarily for square-shaped junctions made from the same trilayer B with various junction 
areas. The critical current, estimated from the normal state resistance by the Ambeagaokar-Baratoff 
relation (AB), conforms with the most likely switching current (ML), which is slightly less than the 
fitted result of switching current measurements (SCM) due to thermal excitation. The area variation of 
junctions in quasi-overlap and inline geometry is exclusively covered by different lengths £, whereas the 
width is kept constant at w = 5 pm. The critical currents are linear to the junction area for the quasi- 
overlap and squared geometry, which goes along with homogeneous current distribution. In contrast, 
the critical current of inline junctions saturates for large junction lengths due to the inhomogeneous 
current distribution. 


As mentioned in Eq. (4.1c), the escape probability determines the junction’s activa- 
tion rate, whose discrete version is calculated from the experimentally acquired 


values Pl ) (see Fig. 4.4) 


au. Lee 
dt j>k 
res (Ip) = a Z pes (4.3) 
jok+1 ] 


Within the scope of the RCSJ model (see Sec. (2.4.1), underdamped junctions switch 
to the nonzero voltage state when overcoming the potential barrier of the junction’s 
tilted washboard potential Up(I) given by Eq. (2.67). In general, this process can 
either be because of macroscopic quantum tunneling of the phase through the 
potential barrier or because of thermal activation over the barrier (see 


Fig. (4.4). 
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Figure 4.4: Switching current measurements of squared high-impedance junctions. a) To evaluate 
the distribution of switching currents from the zero-voltage to the nonzero voltage state with 20000 
events, the activation rate is determined from the escape probability distribution (epd), whereby the 
linear course indicates that this is thermally dominated. A fit of this thermal activation rate as described 
by Eq. provides the critical current I. and escape temperature T° as free fit parameters. b) The 
escape of the phase from the junction’s potential well can be of thermal and quantum mechanical origin, 
both depending on the bias-dependent barrier height and the thermal population of higher energy 
levels. 


The crossover temperature, below which the escape rate transitions from thermal 
to quantum mechanical origin (where it becomes temperature-independent) is 
given by T* = hw ,/27kg for underdamped junctions |DMCS3l. This is 
not a sharp transition, so even just below the theoretical transition temperature, a 
relevant number of higher levels are still thermally occupied due to the long tail 
of the thermal occupancy distribution. Furthermore, it is challenging to bring the 
electrons in the measurement lines down to the bath temperature so that the real 
system temperature is higher than the bath temperature [Gia+06]. 


Accordingly, in the present measurements of the escape from zero voltage to the 
nonzero voltage state thermal activation must be considered, whose rate is given 
by the attempt rate « wp as well as the ratio of barrier height Uo and effective 


escape temperature T°: [BHL83}/DMC85| 


V2 3/2 
(1) = SP ( aie a! Ey? ( a) (4.4) 


ex 
th T P kgTesc 
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The experimentally determined activation rate of Eq. is fitted to this thermal 
activation rate to determine the fitting parameters, namely escape temperature T* 
and critical current ko. 


Using square contacts as examples, critical fitted current values Ieo (SCM) and 
most likely switching currents Ie (ML), both determined by switching current 
measurements, are compared with those estimated by normal resistance according 
to the Ambegaokar-Baratoff relation of Eq. [AB63]. As shown in Fig. 
the Ambegaokar-Baratoff (AB) estimations are in good agreement with the most 
likely switching currents Ie, which are systematically just below the fitted critical 
currents leo. 


Figure (4.3 a)]also shows the critical currents for different junction geometries, each 
with different junction areas £w, only dependent on £ if w is kept constant for long 
junction geometries. It can be seen that the critical currents of the quasi-overlap 
junctions increase linearly with the junction’s length, coincide with the squared 
junctions, and thus exhibits a homogeneous current distribution as desired. In 
contrast, the critical currents of the inline junctions saturate already above a length 
of 40 um (corresponding to an area of 200 pm? at a fixed width of 5pm), which 
allows a Josephson length in the order of 20 um to be roughly estimated. The 
critical current densities listed in Tab. /4.1|are determined for squared junctions. 


In addition to the critical current, the voltage at which the junction jumps is also an 
important parameter for IV characteristics of Josephson junctions. After the escape 
from the zero-voltage to the nonzero voltage state, the phase particle runs down the 
washboard potential, whereby an equilibrium velocity of energy gain due to the tilt 
and loss due to dissipation emerges. These losses increase dramatically when the 
voltage is reached, which is just enough to break Cooper pairs into two dissipative 
quasiparticles each, so that the equilibrium velocity is approximately equal to 
this voltage. As shown in Fig. [2.11 c)} this voltage corresponds to V, = 2A/e 
for SIS junctions with identical superconductors with gap A. In case of S1IS2 
junctions with different superconductors for the bottom and top electrode with 
gaps A, and Aj, respectively, the gap voltage holds V4 = (A, + At) /e. In the 
latter case, a second singularity is possible at V = |A, — 4| /e BP82]. To 
observe these two singularities, the IV characteristics are recorded up to the critical 
temperature, the two effective gaps are extracted, and if both are measurable the 
two superconductors’ gaps are calculated by 


f |Ap — Atl 
AT = a ed Ay = As + AS (4 5) 
S l 
neff _ [Ab za Ay = AFE — age. 
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Figure 4.5: Temperature-dependent IV-characteristic and extracted critical current and gap energies. 
a) Since superconductivity is suppressed with increasing temperature, both the gap voltage, which 
is a measure for the Cooper pair density, and superconductor’s overlap, the critical current, decrease 
and finally disappear completely for T + Te. The increasingly thermally excited quasiparticles thereby 
increase the dissipation in the system, which b) translates into a higher Stewart-McCumber parameter, 
thus a larger retrapping current, and c) in a reduced superconducting gap. 


The resulting gaps of the bottom electrode Ay ~ 265 eV-280 eV and the top electrode 
At © 190 eV-230 eV of the three different trilayers A, B, and C corresponds to that of 
granular aluminum oxide and that of by AlO, proximitized aluminum, respectively. 
The temperature dependence of both superconductors can each be described with 
the Gorter-Casimir two-fluid model 


(4.6) 


The gap values for the different trilayers listed in Tab. |4.1]are extrapolated to T = 0 
from such fits of A;(T), as for example depicted in Fig. 


4.4 Magnetic diffraction patterns 


The high-impedance junction’s characteristic lengths along and perpendicular to 
the tunnel barrier, the Josephson length A; and the magnetic thickness A, respec- 
tively, are experimentally determined from the critical current’s dependence on 
in-plane magnetic fields. As explained in Sec.|2.4.2| externally applied magnetic flux 
causes a spatial variation of the junction’s phase difference along its x-direction, 
which leads to a diffraction pattern of the critical current intensity. The linear 
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Figure 4.6: Magnetic diffraction pattern of a 110 um x 541m long junction of trilayer B in inline 
geometry at different temperatures. a) The linear decrease of the Meissner phase confirms reaching 
the LJJ limit and the extrapolated root of the main lobe corresponds to +H.ı. The asymmetric lobes 
of the magnetic field-dependent critical current arise because of inhomogeneously distributed bias 
currents in inline junctions in combination with different inductances in the two electrodes 
(This figure is also published in Ref. [Wil+22].) b) The maximum of the main lobe 
decreases for higher temperatures according to the Ambegaokar-Baratoff model of S,IS2 junctions 
with different superconducting gaps A;(T), whereas the Josephson length increases because of the 
temperature-dependent kinetic inductance. 


slopes of the lobes in such patterns, as exemplified by Fig. confirm that the 
junction can be considered as long, meaning that self-field effects are not negligible. 
Accordingly, Meissner-Ochsenfeld-like screening currents can shield the junction’s 
interior from small fields. Since these screening currents penetrate deeper and 
deeper along the tunnel barrier from the edges, only smaller and smaller additional 
bias currents can be carried so that the Meissner phase decreases linearly when the 
magnetic field is increased from zero. A long junction maintains this behavior up 
to its first critical field 
Po 


Ha =—, 4.7 
cl THAN ( ) 


at which a quantum of flux penetrates the junction by overcoming its edge barrier. 
Far above, in large in-plane fields, the junction is considered to be completely 
penetrated by magnetic flux, such that flux quanta are laterally lined up close 
together along the tunnel barrier. The flux-penetrated extent in the z-direction, the 
magnetic thickness is accordingly given as 

Po 


u = 4. 
TAB,’ (4.8) 
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where AB, is the critical current’s periodicity of higher side maxima. As can be seen 
in Tab. Lx affects A, since the proximitized top electrode’s London penetration 
depth A, enlarges with increasing Lx, whereas the bottom electrode of each stack 
is without exception in the thin-film limit d4 < A11, and thus contributes to A with 


dı /2 [Wei69]. 


The first critical field H.ı can be determined by the extrapolated zeros of the linear 

slopes of the main maximum. By inserting this value and the already determined 

magnetic thickness A in Eq. (4.7), the Josephson length can finally be calculated as 
Po 

T oHa l 

Using the definition of the Josephson length from Eq. (2.73), the kinetic inductance, 

contributing locally to Aj, can be determined. The comparison of this value g(r)L; 


A = (4.9) 


with the kinetic inductance of the bottom layer Lg, estimated from the resistance 
Ry measured in situ as R(d) during the sample deposition, yields the geometry 


factor g(r) on the order of 1071, as given in Tab. [4.1] [Wil+22]. 


As can be seen in Fig./4.6 b); the Josephson length enlarges with increasing tempera- 


ture since both the critical current and the electrode’s dominating kinetic inductance 
strongly depends on the temperature. The data fit the empirical dependency 


3 21/4 
dere: (: = (=) (4.10) 


that is motivated by the temperature dependence of the two-fluid model with the 
4 a 
characteristic factor (1 = (E) ) . The temperature dependence of the main lobe’s 


maximal critical current is again in good agreement with the model of Ambegaokar 


and Baratoff of Eq. (2.59) [AB63] as depicted in Fig.|4.6 b) 


4.5 Microwave induced switching 


The Josephson plasma frequency is the small-amplitude oscillation frequency of 
a virtual phase particle in a potential well of an unbiased Josephson junction, as 
illustrated in Fig.|2.10 a)| Since the junction’s potential is tilted by a bias current, 


the oscillation frequency strongly depends on the bias as wo(Y) = wp (1 — yy" 


| Dah+68] (see Sec. |2.4.1). 


An off-resonant external AC bias merely yields weak oscillations of the potential 
tilt. It thus adds slightly to the DC bias so that smaller DC currents suffice that 
the junction escapes from the potential well. However, if an additional microwave 
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Figure 4.7: Resonance current at different external microwave drives of a squared junction with 
(10um)? of trilayer B at 20 mK. a) The junction’s phase oscillation frequency in a potential well 
strongly depends on the bias current. Forced oscillations by microwave irradiation can provoke a phase 
particle’s escape from such a potential well if its resonance frequency or n subharmonics matches 
the drive frequency. b) If so, the junction’s switching current becomes multivalued, and the secondary 
peak in the switching current distribution is identified as resonance current, where its full width half 
maximum is the error. (This figure is also published in Ref. [wil+22]. 


tone drives the oscillations resonantly, subharmonically or superharmonically, a 
significantly smaller DC bias current is necessary to escape from a potential well 
and reach the running state characterized by a nonzero voltage [Gro+04]. This 
shows up as a secondary peaks in the switching current distribution 
|BCG10], as exemplified in Fig. that is identified as resonant cur- 


rent, for which the fixed external drive frequency equals the junction’s resonance 
frequency wo, its integer multiples or fractions of wo. 


In the experiment with a 10 um x 10 pm squared junction of trilayer B as an example 
[Wil+22], switching current measurements with 10000 events are taken under 
microwave irradiation of a variable but fixed frequency. After finding a multivalued 
switching current distribution, the power of the microwave excitation is adjusted so 
that the additional peak is still clearly visible but does not yet shift towards smaller 


currents |Wal+03a]. This is the case if the integrals over the individual peaks are 
nearly equal. 
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4.6 Current steps 


The current value of the additional peak and the corresponding frequency are a 
single data point in Fig. and for many such measurements branches of the 
different harmonic n can be found, obeying 


2\ 1/4 
wn(y) = = € = (+) ) : (4.11) 


The extrapolation of the oscillation frequency with orthogonal distance regression 
to I = 0 yields the plasma frequency wp/27 = (13.28 + 0.05) GHz and the critical 
current I, = (8.36 + 0.08) pA. 


This fit enables the calculation of the tunnel barrier’s specific capacitance using 
the definition of the Josephson plasma frequency in Eq. (2.64). The value cy = 
(36.4 + 0.4) fFum~? is a typical number for tunnel barriers made from thermally 
oxidized aluminum and indicates that the additional kinetic inductance 
does not affect the plasma frequency to the first order. 


4.6 Current steps 


Apart from the static measurements on long junctions in applied magnetic fields 
to determine the Josephson length and the measures on short junctions to de- 
termine the plasma frequency, their product, the Swihart velocity © = Ajwp, is 
determined in independent measurements separately. As explained in Sec. [2.4.3] 
extended Josephson junction stacks can be viewed as nonlinear transmission lines 
through which different excitations, such as plasma waves or solitons in the form 
of Josephson vortices, can propagate with at most the Swihart velocity ¢ [Ped86]. 
The impedance mismatch at the two open ends of a one-dimensional long junction 
forms a cavity that allows standing wave excitations. 


One-dimensional linear long junctions with open boundaries have no topological 
constraints on the number of inherent vortices and antivortices. Therefore, fluxons 


can be trapped in such junctions even in zero magnetic fields?| where the overlap 

junction’s boundaries hold (gz + boz) | = (9: + Bos) | = 0. In this geome- 
= x= 

try, solitons are driven shuttle-like along the tunnel barrier by a Magnus-Lorentz 


force caused by a nearly homogeneously distributed bias current. At the junction 
ends, they are reflected, reversing not only their direction but also their polarity to 


? At first sight this may seem counterintuitive, but theoretical considerations show that in long junctions 
even the smallest perturbations are sufficient for a Josephson vortices to be energetically more 


favorable [Lik86]. 
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Figure 4.8: Current singularities due to junction cavity resonances of 100 um x 5 um and 110 um x 
5um long junctions of sample B in quasi-overlap geometry without and with magnetic fields, 
respectively. a) Zero-field steps arise only with sufficient damping in the junctions, which is realized by 
temperatures just below the critical temperature Te ~% 1.25 K. b) Fiske steps occur at different magnetic 
fields and their characteristic rounded shape originates from the increased damping at T = 1.0K. Both 
current steps feature a bias current depending characteristic shape according to Eq. 4.12), whose fits 
and errors as displayed as dark gray line and the gray shaded area. (These figures are also published in 


Ref. [wil+22].) 


satisfy the boundary condition. This process can be regarded as a virtual collision 
with an antisoliton. Thereby, soliton and antisoliton move towards the junction 
boundary with opposite velocities, collide at the edge and pass through each other 
if their velocities are large enough, such that the soliton leaves, whereas the anti- 
soliton enters the junction [Ust05]. Due to its reverse polarity, this antisoliton is 
driven in the opposite direction of the soliton. In doing so, the soliton is reflected as 
antisoliton at the one junction boundary and vice versa at the opposite boundary, 


yielding resonances in the junction cavity in zero magnetic fields [CFL71; |FD73 
FD74b]. 


Accordingly, two 27t-kinks pass each point of the junction per period At, giving rise 
to the time-averaged voltage (V) = 2@g/At across the junction. The period is given 
by the junction length and the propagation velocity from Eq. (2.98) as At = 2¢/v, 


so that the voltage drop is given for 6 = 0 by |[MS78}|BP82}|Lik86] 


Doe n 
v (4) = 


£ 1+ (8) 


Here, n denotes the number of participating solitons that cause equidistant current 
steps in the IV curves, which approach the voltage ®9C/£ in the relativistic limit 


(4.12a) 
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4.6 Current steps 


(v — č) achieved for large bias current. For their experimental observation, ho- 
mogeneous bias current distributions such as those obtained in the quasi-overlap 
geometry are advantageous to ensure a shuttle-like soliton motion and avoid their 
edge pinning as observed in inline junctions. 


To reliably access such current singularity features in the IV characteristic arising 
from the subgap resistance branch, it further helps to increase the dissipation of the 
underdamped junctions by increasing the sample temperature. Then, however, as 
shown in Fig.|4.8a)} the Stewart-McCumber branch cuts the lower part of the higher- 
order steps. Nevertheless, the current steps are fitted to Eq. to determine the 
Swihart velocity. 


In finite magnetic fields, this perception of shuttle-like vortex motion remains valid 
only in a very narrow parameter range of low magnetic fields H ~ Hq /2, 
intermediate junction lengths £ S 10, and comparable large damping [EFP83]. In 
this configuration, fluxons can travel the Josephson transmission line without being 
reflected at the junction edges, but being rather annihilated at a boundary. The 
released energy travels back as a slowly decaying plasmon, which in turn triggers 
a fluxon at the opposite boundary [Due+81}/Ust05]. Consequently, a single 27 per 
fluxon kink passes each point of the junction per full period, which results in half 


the voltage of the zero-field steps |Fis64}|Kul65;|Kul67 


_ Doe n 


(4.12b) 


Thus, such Fiske steps with even numbers of participating fluxons n coincide exactly 
with the former zero-field steps. Fiske steps can also be observed apart from this 
narrow parameter range with identical characteristics, although the underlying 
mechanism differs. In general, the Josephson frequency of a biased junction excites 
electromagnetic standing waves in the junction cavity, whose higher harmonics 
give rise to periodic current singularities, as shown in Fig. For underdamped 


junctions, cavity oscillations are unstable for w S wp [Cir+97 , which 


explains missing the first Fiske step in Fig.|4.8 b) 


Nevertheless, the Swihart velocity ¢ is determined by the periodicity of the current 
singularities and the well-defined junction length £ (see Tab. 4.1). Knowing the 


junction width w from fabrication and the specific capacitance cy, the impedance 
1 
wecy 


of the Josephson medium can be calculated by Zy = 


To analyze the impact of the electrode’s kinetic inductance on LJJs, their character- 
istic parameters, listed in Tab. are compared with estimations for conventional 
LJJs with equal tunnel barrier properties je and cy, but made from pure aluminum. 
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4 Fluxons in high-impedance long Josephson junctions 


Here, we assume that pure aluminum electrodes have negligible kinetic inductance. 
The result of this comparison is that the electrode’s kinetic inductance reduces both 
the Josephson length A; and the Swihart velocity ¢ by a factor of up to 40, while the 
Josephson plasma frequency wp remains nearly unchanged. Accordingly, the wave 
impedance of LJJs is increased by the same factor. The inductance contributing 
to Josephson plasma oscillations is dominated by the macroscopic stack tunnel 
barrier rather than the nanoscopic tunnel barriers in AlO, due to the much stronger 
intergrain coupling so that the increase in Lẹ can be neglected to the first order. 
The combination of the independently measured parameters corresponds to the 
conventional sine-Gordon model with © = Aywp. 


4.7 Concluding remarks 


The reported transport measurement results demonstrate that the novel high- 
impedance Josephson junctions made from disordered oxidized aluminum as a 
high-kinetic inductance superconductor are promising candidates for studying the 
quantum coherence of Josephson vortices. These long junctions feature high quality. 
The observation of zero-field steps proves the existence of Josephson vortices in 
the long junctions, and thus confirms the sine-Gordon model being applicable for 
high-impedance long junctions. 


Beyond, the experiments demonstrate a decrease in the vortex’s size and a reduction 
of its limiting (Swihart) velocity by about one order of magnitude compared with 
conventional LJJs, whereas the plasma frequency is unaffected to the first order. 
The measured Swihart velocities down to a small fraction of 3 x 1073 of the vacuum 
light velocity, in turn, correspond to an increase in the junction’s wave impedance 
up to 140) compared to 40) of conventional, similarly made LJJs. The enhancement 
of this impedance goes along with an increased quantumness Åy « Z. As hj plays 
the key role in experimentally reaching the quantum regime of Josephson vortex 
dynamics [Wal+03b], high-kinetic inductance electrodes also facilitate observing 
the quantum electrodynamics phenomena in long Josephson junctions. 


Apart from the goal in this work of quantum coherence of Josephson vortices, 
these high-impedance junctions are promising for the research field of active 
Josephson microwave components. The high-kinetic inductance electrodes enable 
tailoring the junction impedance and facilitate solving the long-standing problem 
of impedance matching LJJs to external circuits and 500 cables. Matching the 
impedance to external loads is crucial for increasing the efficiency of Josephson 
flux-flow oscillators used for microwave generation and amplification [Wil+22]. 
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5 Long junctions embedded in a 
high-impedance resonator 


Observations of quantum coherence are experimentally challenging using DC 
measurements since the sample is galvanically connected to the environment, 
representing a massive decoherence channel. However, excitations with alternating 
electromagnetic fields allow for manipulation and readout of the circuit with- 
out such DC lines. For this purpose, long junctions are coupled to high-quality 
microwave resonators that are completely galvanically decoupled from the envi- 
ronment and act as very narrow band-pass filters even for AC by coupling only on 
resonance. This readout scheme is well established for superconducting quantum 
bits [Wal+04], but not yet for long Josephson junctions. The following 
chapter elaborates on a method to couple high-impedance long Josephson junc- 
tions (LJJs) to such resonators and their fabrication process. Further, spectroscopy 
measurement data involving the anharmonicity and the magnetic field dependence 
are compared to a theoretical model of coupled junction and resonator. 


5.1 Sample considerations and design 


One possible way to couple a high-impedance long junction with a distributed 
A/2 high-kinetic inductance resonator is to expand the junction’s bottom electrode 
to a resonator. In other words, a section of the resonator is shared as the bottom 
electrode of the long junction, where a tunnel barrier and a top electrode are 


deposited on top (see Fig. 5.1). 


For infinitely long junctions, the currents in this assembly, driven by alternating 
electromagnetic fields, branch out in the two electrodes according to their induc- 
tance ratio as described by Kirchhoff’s current law. Because of the supercurrent 
tunneling through the junction, a resonator-junction interaction is achieved, whose 
coupling strength is controlled by the inductance ratio of the bottom to the top 
electrode as well as the widths of the junction and resonator. Strong coupling can, 
for instance, be realized by a high-kinetic inductance bottom electrode made of 
AIO, and a top electrode of pure aluminum, whose orders of magnitude smaller 
inductance is negligible (see Fig. 5-1). In this extreme case, the major part of the 
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5 Long junctions embedded in a high-impedance resonator 


a) side view Al 
sapphire 
b) top view 


Figure 5.1: Schematic sample design for long junctions embedded in a microwave resonator. a) An 
open-end microstrip LC transmission line forms a distributed A/2 resonator. The center part of this 
high-kinetic inductance resonator serves as bottom electrode of a high-impedance long Josephson 
junction, whose top electrode is made from a pure aluminum. b) The resonator is coupled capacitively 
to a microwave feedline in notch type configuration. 


current driven in the resonator tunnels through the junction and flows in the upper 
branch. 


The extent of the junction’s edge regions, where this current tunnels across the 
barrier, is governed by the ratio between Josephson and electrode inductance, and 
thus by the dimensionless screening parameter fr = Lo/Lj, which is related to the 
Josephson length. 


As distributed resonators exhibit spatially varying electromagnetic field distribu- 
tions, the choice of the LJJ’s position is decisive: Assuming a homogeneous kinetic 
inductivity, the spatial center corresponds to the symmetry point of the current 
distribution in the resonator. This means there are current nodes for even and 
antinodes for odd harmonics, respectively. 


For this reason, the LJJ’s impact on the resonator varies depending on the excited 
mode since the phase drop along the junction Ab « LI differs. 


FEM simulationg!| analyze the impact of varying inductance in the resonator’s 
center part (this is where the LJJ will be placed) on its resonant frequency. In 
the extreme case of zero inductance of the inner 10% of the resonator area even 
harmonics are affected insignificantly as the resonant frequency changes only by 
0.03% for the second harmonic and 0.8% for the fourth harmonic. In contrast, 
odd harmonics with current antinodes in the center of the microstrip resonator 


1 The resonant frequencies are estimated by employing Ansys HFSS, a multipurpose, full-wave 3D 
electromagnetic simulation software that utilizes a finite element method (FEM) eigenmode solver. 
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5.1 Sample considerations and design 


drastically change their resonant frequency by ~ 10.4%. Thus, spectroscopy mea- 
surements to determine higher-order resonance frequencies of such samples can 
be used to experimentally verify whether the proposed arrangement couples LJJs 
to resonators, making the former suitable for RF readout. 


As LJJs’ inductances are sensitive to magnetic in-plane fields, such fields should 
significantly impact the resonant frequencies. For this purpose, a model of phase- 
biased junctions was developed in this work, which is described in detail in Sec. 


To attribute the expected anharmonicity and modified magnetic field dependence 
to the LJJ, a reference resonator without any LJJ on top is fabricated in the same 
fabrication run. Frequency multiplexing of several notch-type resonators on a 
common feedline is utilized to measure multiple resonators with a small spatial 
distance to be as insensitive as possible to inhomogeneities in the trilayer. In 
addition, fewer RF lines and cooldowns are needed to measure and compare 
different resonators, but the measured resonances need to be unambiguously 
identified. 


Two-dimensional on-chip waveguides are usually realized in a coplanar structure, 
where center conductors are separated by an insulating gap from semi-infinite 
ground planes or by microstrip lines, whose ground electrodes are placed on the 
chip’s backside. Since homogeneous in-plane fields are applied in the present ex- 
periment, the microstrip variant with a normal instead of superconducting ground 
plane will be used. Normal metals prevent bending magnetic fields near the junc- 
tions due to Meissner screening in the ground plane. Besides, a normal conducting 
backside metallization avoids trapping Abrikosov vortices, which would interfere 
with the applied magnetic field. 


In the AC experiments, it is convenient to use microwave frequencies on the 
order of a few gigahertz since their corresponding energies are, on the one hand, 
sufficiently below the superconducting gap (e.g., 2A/h = 87 GHz for aluminum) 
to avoid quasiparticle excitations. On the other hand, gigahertz frequencies are 
adequately above the available temperatures (kgT/h ~% 200 MHz for conventional 
dilution refrigerators reaching T ~ 10mK) to dominate thermal excitations and 
enable deliberate manipulation. Moreover, the manifold commercial electronic 
equipment used for telecommunication operates in the kilohertz to gigahertz band. 
Access to this microwave technology required for the RF measurements, such as 
vector network analyzers, microwave sources, amplifiers suitable for cryogenic 
temperatures, circulators, mixers, etc., ease experimental implementations. 
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5 Long junctions embedded in a high-impedance resonator 
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Figure 5.2: Schematic patterning of a long junction embedded in a microwave resonator. The upper 
row depicts the circuit’s top view, and the lower row the side view along the resonator and long junction 
(dashed line) after each fabrication step. Starting from the deposited junction stack (i) the resonator 
area that includes the junction is patterned by an argon-chlorine-based inductively coupled plasma 
(ICP) etching process. The top electrode of the remaining trilayer (ii) is also removed by ICP etching at 
the outer parts of the resonator to define the long junction in the middle of the resonator only (iii). 


To achieve a large coupling strength between resonator and junction, granular 
aluminum oxide with RZ = 2 kO which corresponds to Lg = 1.5 nH, as bottom and 


pure aluminum as top electrode are desired. In order to exhibit fundamental modes 
in the frequency range 4.5 GHz-7.5 GHz, the resonators’ lengths are designed to 
l ~ 300 pm-500 um at a uniform width of w = 2m. The junctions extend over 
the entire width of the resonators w = w and their lengths are fixed to 10 % of the 
resonator length, so that £ ~ 30 um-50 um. While these values should correspond 
to long normalized junction lengths £ ~ 9-15, the Josephson length is targeted 
around A; ~ 3.5 pm. Therefore, the desired critical current density is estimated to 
be je ~ 100 A cm™?. 


5.2 Sample fabrication details 


The trilayer D, whose measurements are shown below, meets the desired parameters 
(see Tab. and is sputtered as described in the Sec. Resonators and 
long junctions are patterned each by photolithography and argon-chlorine-based 
inductively coupled plasma (ICP) reactive ion etching (RIE), as schematically 
illustrated in Fig. First, areas of the resonators, including the long junctions, are 
covered by a photoresist, and the rest of the trilayer is subsequently etched down 
to the sapphire substrate. Second, only the smaller junction area in the center of the 
resonator is covered by a photoresist, and the top electrode layer only is removed 
next to it so that the outer part of the resonator consists only of the AlO, layer. This 
step is delicate because it must be ensured that the junction, i.e., the low-inductance 
aluminum top electrode, is completely removed while the plasma does not etch too 
deeply into the bottom electrode. Unfortunately, since there is no selective etch that 
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Figure 5.3: Atomic force micrograph of an AlO, strip with an LJJ on top. The atomic force micrograph 
(a)) shows the long junction embedded in a microstrip resonator, and the schematic of the microstructure 
is illustrated in b). From the measured height profile, the width of the vertical trilayer strip is estimated 
to = 1.7 pm. In the lower part (y S 9.5 um), the long junction’s trilayer stack survives, whereas the top 
electrode is etched away in the upper region. For this purpose, a photoresist, which protects against 
etching, is patterned photolithographically with a horizontal excess length of ~ 21m, resulting in a 
slightly higher substrate there. 


etches only pure aluminum but not granular aluminum oxide, this etching result 
can only be achieved by choosing the appropriate etching time. See Appendix |B.2| 


for details of parameters concerning optical lithography (Tabs. and |B.4) 
and ICP etching (Tab. |B.5). 


Finally, an additional 100 nm thick layer of niobium is deposited onto the microstrip 
feedline in an optical liftoff process. Niobium features high critical temperature, 
high critical magnetic in-plane fields, and negligible kinetic inductance that allows 
for well-defined impedance matching to the standard microwave components with 
500. To confine the electromagnetic fields properly, a backside metallization is 
deposited in the form of a 100nm thick layer of silver with a 5nm thick coat of 
titanium beneath to improve its adhesion on the sapphire substrate. 


During patterning of sample D, the top electrode was slightly overetched, that 
is to say the granular aluminum layer was thinned down by a few nanometers. 
As the thickness of the resonator layer decreases, the sheet resistance or rather 
the kinetic inductance increases (compare the Rn(d) curve in Fig. [3.2 b)), which 
can be verified by test structures at room temperature. Four-point probes on test 
strips with 150-250 squares yield a sheet resistance of Rọ = (4.6 + 0.6) kQ, which 
corresponds to Ly = (4.9 + 0.6)nH according to Eq. (2.32). FEM simulations 
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5 Long junctions embedded in a high-impedance resonator 


with this estimated kinetic inductance value yield fundamental modes between 
2.6 GHz-4.3 GHz for the six fabricated resonators, whose dimensions are listed in 


Tab.|C.1lin Appendix 


The height profile from atomic force microscopy of a typical fabrication results 
of this process is depicted in Fig. 5.3} There, the structures can be measured in 
all three spatial directions so that the width of the resonators and junctions is 
determined to be about * 1.7 ym. From the height profile, it can be confirmed that 
a thin layer of the top electrode was removed in the second etch step. The window 
extending beyond the junction comes from the fact that the resist in the second etch 
step protrudes about 2 um beyond the functional structure onto the bare substrate 
to become insensitive to misalignment along this direction. 


5.3 Phase-biased Josephson junctions 


To determine the eigenfrequencies of the coupled system of the long junction and 
microwave resonator, the eigenmode solver of the electromagnetic FEM simulation 
suite Ansys HFSS is deployed. For this purpose, the inductances of the resonator’s 
parts with and without the junction are treated differently. In the outer regions of 
the bare AlO, layer, the dominant kinetic inductance is estimated from Eq. (2.28). 
The spatial phase distribution along the junction is analyzed in the resonator’s 
critical center part with the LJJ. The inductance, included in the FEM simulation, is 
then calculated from the phase drop between the two junction ends. 


For this analysis of the phase distribution, the inductance in the top electrode LA! 
is neglected since it is small compared to the inductance in the bottom electrode 
La > ENL As shown in Fig. the distributed LJJ is divided into N — co 
discrete segments and Kirchhoff’s current law is applied at each node j € [1, N]. In 
doing so, inside the junction, the j™ node flux is coupled by the nearest neighbors 
as 


a a a E a a a 
Lo/ (N —1) N Dy ! Lo/ (N —1) 
At the junction edges, the bias current I, driven by the resonator, are bought into 


the equation, so that the boundary conditions hold 


(5.1a) 


P- Pi- OOM (2m 
|) =I 1b 
Lo/ (N-11) No AG, r (d.1b) 
k . /27 Py — Pn- — OY + OY, 
fa ET by) = . 61 
yo (= N :) Lo/(N— 1) (5 c) 
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5.3 Phase-biased Josephson junctions 


QH 


Figure 5.4: Schematic circuit of a galvanically coupled long junction to a microwave resonator. The 
microstrip structure forms an open-end LC transmission line with the strip’s geometric and kinetic 
inductance Lo = L6/ and the capacitance against the backside ground electrode C = C’l. The resonator’s 
conductor is shared as the long junction’s bottom electrode, which is weakly coupled with the top 
electrode by an insulating tunnel barrier with the critical current I. = Jel. Since the bottom electrode’s 
kinetically dominated inductance exceeds the purely geometric inductance of the top electrode by 
orders of magnitude, this is dropped in the schematic circuit. 


The node fluxes can be translated into node phases by means of the junction’s 
well-known flux-phase relation ®; = Fo. Additionally, a homogeneous external 
magnetic flux causes a seas external phase drop between two adjacent nodes 


= ont = x2 ar. 1 = DAY Y j. By further identifying the Josephson induc- 


tance Ly = = and introducing the screening parameter By, = Im = en (for 
details see Appendix[A.2), the conditional equations can be formulated as 


= j ext | BL . BLY 
92 = 91 + Ag™ 4 NNI) sin (90) + Nov’ (5.2a) 
Pj+1 = 29) — Yj-1 NOT D sin (gj) , (5.2b) 
our = pn — pn-1 Ay 4 nt T) sin (pn). (5.20) 


Starting from the phase at one junction boundary 91, the phase develops according 
to these recursive formulas depending on the bias current I and the external 
phase Ag™*. 


For numerical calculations with various but ae I and Ag, we guess values 
of the n phase 9ı € [0,27), apply the Eqs. ( and (5.2b), and verify by 
Eq. (5.2c), whether the chosen initial phase an to the boundary conditions 
or not. z valid solution stands out due to current conservation, which means that 
the injected current at the one junction boundary equals the ejected current at the 
other one. 
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Figure 5.5: Numerically simulated anharmonicity and magnetic field dependence of a phase-biased 
junction. a) Odd harmonics of hybrid systems of resonators with LJJs (stars) feature significantly 
increased eigenfrequencies compared to reference resonators without any junction (dots). b) This 
fundamental distinction is also evident in the resonant frequencies’ magnetic field dependence, whose 
odd harmonics are affected way stronger by the LJJ than even harmonics. 


By running a nested interval algorithm, multiple mathematical solutions are found, 
where the energetically most favorable solution is picked as the physical solution. 
The phase drops along the junction gy_1 — Po for different bias currents yields 
a current-phase relation I(Ap), which is fitted to a fifth-degree polynomial. The 
coefficient of the linear term equals the reciprocal of the inductance, which is 
inserted in the FEM simulation to obtain the resonant frequencies. 


The simulated resonant frequencies of the first four harmonics, as well as their 
magnetic field dependence, are depicted in Fig. Both results show that the 
coupling strength between LJJ and resonator depends on the excited harmonic, 
and thus confirm the earlier discussed considerations. 


5.4 Spectroscopic characterization 


The transmission spectrum of the complex S21-parameter exhibits numerous res- 
onances, each as described in Sec. These are different harmonics of the 
frequency multiplexed resonators, which first have to be assigned to the corre- 
sponding circuits (see Fig. |5.6a)). The fundamental modes match the results of 
FEM simulations and are governed by the differing resonator length /, as width, 
capacitance per unit length, and sheet inductance are constant. 
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5.4 Spectroscopic characterization 


Table 5.1: Measured and simulated anharmonicity of resonators with and without an embedded 
long junction in zero magnetic fields. The anharmonicities of the first three modes from spectroscopy 
measurements fit exactly to the Ansys HFSS FEM simulations in the case of the bare reference resonator 
with a sheet inductance of L- = 5.649 nH. Assuming the same inductance for the AlOx layer of the 
resonator and the junction’s sheet inductance of L- = 0.554nH, the FEM simulations show good 
agreement with the measured resonance frequencies. 


circuit fol fi fs/fi 


experiment simulation experiment simulation 


bare resonator 1.915 +0.005 1.915 0.003 2.807 +0.005 2.809 + 0.012 
resonator & LJJ 1.768 +0.005 1.759 0.004 2.76 +0.02 2.807 + 0.014 


The higher harmonics are attributed to the fundamental modes by two-tone spec- 
troscopy measurements, whose schematic setup is explained in Sec. Thereby, 
the fundamental mode is probed, while the second tone’s frequency is swept with 
a way higher power. If the second tone is in resonance with any harmonic, the 
resonator is populated so much that its fundamental mode shifts to lower frequen- 
cies due to the nonlinear inductance. Using this Kerr effect and the well-known 
resonator geometries, the resonant frequencies can be attributed to the frequency 
multiplexed resonators, as shown in Fig. 


The reference resonator (#3) without any coupled LJJ can be unambiguously iden- 
tified since it shows almost no magnetic field dependence at fields far below He2 of 


AlO; (see Fig. [C-T]in Appendix |C) |Bor+20]| and as its harmonics match FEM sim- 
ulations as listed in Tab. 5.1] From the comparison of the measured and simulated 


resonance frequencies, a sheet inductance of L Alo, = 5.469 nH can be estimated. 
This value agrees with the resistivity measurements of the test structures at room 
temperature and originates in the already discussed overetching of the sample. 


5.4.1 Anharmonicity of the coupled system 


Their relative spatial current distributions for the other five coupled systems of 
resonators and LJJs are identical since each LJJ covers the inner 10 % of the resonator. 
This similarity is why the fundamental modes of these coupled systems show a 
uniform length dependence, as can be seen in Fig. Compared to the reference 
resonator, the fundamental resonant frequency is increased by about 7.3 %, which 
is in good agreement with the FEM simulations described in the model of phase- 
biased junctions in Sec. 


In contrast, both resonators with and without an embedded LJJ feature a uniform 
length dependence of their second harmonic frequency. This is because the low- 
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Figure 5.6: First three harmonic resonance frequencies of different resonators with and without LJJ 
in the zero-field and their magnetic field dependence. a) The measured resonant frequencies of the 
coupled systems of the resonator and long junction in zero magnetic field show a length dependence 
illustrated by straight lines to guide the eye. Unlike the second harmonic of the reference resonator 
(#3), the resonant frequency deviates from this course for odd harmonics. This different behavior of 
circuit #3 is apparent from the anharmonicity plotted as stars on the right ordinate. b) The magnetic 
field-dependent amplitude signal of a resonator with LJJ (exemplarily shown for circuit #2) exhibits a 
significantly larger shift to lower resonant frequencies for odd than for even harmonics. 


inductive shunt of the junction is placed in the current node, and thus modifies the 
phase drop along this segment only slightly. 


The third harmonics show an increased resonant frequency similar to the first 
harmonics because there is again a current antinode at the position of the LJJ. 
Unfortunately, not all third harmonics of the six circuits are measurable, but the 
remaining four values resemble the theory. The fourth and higher harmonics are 
regrettably not measurable. 


A comparison of the average measured and simulated anharmonicities of the 
resonators with and without LJJ show good agreement between theory and ex- 
periment, as listed in Tab. The slight deviation originates in the assumption 
for the simulation of the long junctions’ homogeneous sheet inductance that, in 
reality, spatially differs. Especially at the long junctions’ edge regions where the 
supercurrent mainly tunnels through the tunnel barrier, the inductance peaks as 
the vertical Josephson inductance adds to the lead inductance. 
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5.5 Concluding remarks 


5.4.2 Magnetic field dependence of the resonant frequencies 


An LJJ embedded in a resonator affects not only the anharmonicity but also 
the magnetic field dependence, which was theoretically discussed in Sec. 
Microwave spectroscopy data confirm a comparable dependency experimentally, as 
exemplarily depicted for the circuit #2 in Fig. Since the reference resonator 
does not feature such a field dependence, but is constant in its resonant frequency 
far below Ho (see Fig.|C.1|in Appendix|C}, it is reasonable to attribute this effect 
to the additional LJJ. 


5.5 Concluding remarks 


These proof of principle measurements show that LJJs can be galvanically cou- 
pled to microwave resonators by extending the bottom electrode to a microstrip 
resonator that biases the junction at microwave frequencies. The combination of 
LJJ position and current distribution of the excited resonator mode is crucial for 
the resonant frequency of the coupled system and its magnetic field dependence. 
Spectroscopy measurements of these two parameters show good agreement with 
FEM simulations, in which the phase drop along the junction is numerically cal- 
culated by the model of a phase-biased junction. These simulations further reveal 
that the coupling strength of the two systems varies with the inductance ratio of 
the bottom and top electrodes and the width ratio of the junction and resonator. 


This RF manipulation and readout scheme has been used for the first time for long 
junctions and offers tremendous advantages for experimentally investigating quan- 
tum coherence effects in long junctions. Since the resonator is a narrowband Purcell 
filter and the circuit is galvanically isolated from the environment, this strong 
decoherence channel is significantly reduced compared to DC measurements. 
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6 Josephson vortices as quantum 
particles 


In the previous chapters, we could demonstrate by transport measurements (see 
Sec. 4) that the junction impedance, and thus the quantumness of long junctions, 
can be increased with additional kinetic inductance in the junction electrodes. With 
the RF measurement scheme described in Sec.|5} an appropriate method was found 
for quantum coherence measurements. We are now ready to investigate the vortices’ 
quantum character with these prerequisites. For this purpose, very narrow junctions 
are patterned from high-impedance trilayers and coupled to a microwave resonator 
each. These samples are characterized by steady-state microwave spectroscopy and 
pulsed microwave measurements. 


6.1 Sample considerations and design 


In principle, a similar sample design is used as described in Sec. 5.1} but the LJJs are 
significantly narrowed to increase the inherent vortices’ quantumness. A junction 
width of w = 200nm is a good tradeoff between this aimed quantumness and 
reliably producible structure sizes without creating constrictions or broadenings 
that would cause pinning potentials for vortices [SBM97]. Moreover, tens of microns 
long and 200nm wide AlO; strips as junction electrodes are still wide enough to 
ensure a certain phase of the superconducting wave function and avoid phase slips 


[Vos+21]. 


To read out the system’s quantum state using the dispersive shift of the resonator, 
it is mandatory that the resonator’s eigenfrequency remains stable over time. The 
previous measurements of high-impedance resonators (see Sec. By, however, have 
shown fluctuations of the resonant frequency in the order of hundreds of kilohertz. 
Such strong variations are known from different works about superconducting films 
near the superconductor to insulator transition and appear if the sheet resistance 
exceeds the superconducting resistance quantum Rq = h/ (2e)? = 6.45 KQ 
Kri22]. In the circuits of trilayer D, the high sheet resistance Ry = (4.6 + 0.6) kQ 
results from the already large intended value of Rā = 2.3kQ and the overetching. 


Therefore, to reduce frequency fluctuations in new samples, AlO, bottom layers 


with a sheet resistance of Ry = 1kN at an enlarged thickness of da = 50nm 
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Figure 6.1: Schematic patterning of junction stacks for microwave measurements. The functional 
microstructure’s top and side view after each fabrication step are illustrated in the upper and lower 
row, respectively. Starting from the deposited junction stacks (i), the long junctions’ top electrodes are 
patterned in by electron-beam lithography and an argon-chlorine-based inductively coupled plasma 
(ICP) etching process (ii). The resonators (iii) are defined afterwards by optical lithography and ICP 
etching of the remaining bottom electrode, too. 


are targeted. Since the slope of R(d) (see Fig. decreases with larger film 
thicknesses, the sheet resistance of thicker films is less sensitive to overetching. 
According to Eq. (2.28), the desired sheet resistance of RẸ = 1 kQ corresponds to a 
sheet inductance of Lg = 0.75nH. 


Again, a 30nm thick layer of pure aluminum is used as top electrode to strongly 
couple junctions to resonators. Moreover, the coupling strength is varied by design- 
ing resonators with a width of 2 um, 4 um, and 8 um, while the junction width re- 
mains fixed at 200 nm. To obtain resonant frequencies in the range 4.5 GHz-7.5 GHz 
for the fundamental modes, a 750 um and a 900 um long resonator are designed for 
each resonator width. The junctions lengths are still 0.1-/, i.e., 75 um and 90 um. 


6.2 Sample fabrication details 


Due to the similar design, the fabrication process is related to the one of LJJs coupled 
to microwave resonators, which is presented in Sec. The main difference is the 
requirement of electron beam (e-beam) lithography to resolve the LJJs’ structure 
size of 200 nm. 


This process step will be performed before patterning isolated microstrip resonators 
to prevent damaging the LJJ by electrostatic discharge due to the deposited electrons 
during e-beam lithography. Due to the much smaller tunnel resistance over the 
extensive area of 2cm x 2cm, the risk of such damage is reduced significantly. 
The deposited charge can further be led away from the still completely connected 
trilayer much easier via clamps. 
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The photoresist ma-N 2403 is suitable as it is sensitive to both e-beam and deep UV, 
it provides a resolution of = 50 nm, and it is resilient against argon-chlorine-based 
inductively coupled plasma (ICP) reactive ion etching (RIE). Moreover, using this 
negative photoresist saves exposure time since its exposed areas persist in the 
development while the unexposed resist is lifted so that only the tiny junction areas 
need to be written by the electron beam. With this method, illustrated in Fig. 
long junctions are patterned from trilayer E, whose characteristic parameters are 
given in Tab. The resonators are defined in the second process step, using 
photolithography and ICP etching as before. The subsequent deposition of a 
backside metallization is analog to previous samples and consists of 5nm titanium 
and 100 nm silver. See Appendix [B.2| for detailed parameters concerning optical 


and e-beam lithography (Tabs. and {B.4) and ICP etching (Tab. [B.5p. 


6.3 Energy spectrum of vortex tunneling 


Spectroscopy measurements show resonances with internal and coupling quality 
factor Qi = (1.4 + 0.4) x 10° and Qe = 1828+ 4 in the single photon limit, whose 
resonance frequencies fluctuate insignificantly. The overall magnetic field depen- 


dence resembles previous spectroscopy measurements of an LJJ embedded in a 
resonator, reported in Sec. described with the model of phase-biased junctions. 


In addition, pairs of avoided level crossings are observed that appear at particular 
flux densities B,,, as can be seen in Fig. {6.2|and Fig. Such features indicate 
that the driven resonator exchanges energy with a coupled system at a similar bare 
resonant frequency. As long as the coupled system is sufficiently nonlinear, its 
mode spectrum vs. applied magnetic field can be tracked by two-tone spectroscopy, 
similar as explained in Sec. [2.3.3] For this purpose the resonance dip is determined 
at each field value, the resonant frequency is probed with the VNA for and the 
second tone’s drive frequency fy is swept with an additional microwave source 
(for the setup see Fig. 3.5a)). If this extra drive tone excites the second mode, 
the resonance curve shifts dispersively so that both phase and amplitude signals 
change drastically. As shown in Fig. such two-tone spectroscopy data exhibit 
hyperbolas, which are an unambiguous signature of quantum tunneling in a two- 
level system that is associated with a double-well potential. 


The ratios of the magnetic flux density values B,,, at which these pairs of avoided 
level crossings occur, match the intersections of adjacent fluxons states in the 
magnetic diffraction pattern as shown in Fig. Thus, the two states can be 
understood as |n — 1) and |n) vortices in the LJJ, whose eigenenergies are similar, 
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Figure 6.2: Magnetic field-dependent spectrum and two-tone spectroscopy data. a) The overall mag- 
netic field dependence of the samples of trilayer E resembles those of the samples D as expected since 
the basic design is the same, and only the dimensions of the long junction differ. Therefore, the system 
of a long junction embedded in a resonator can further be described by the model of a phase-biased 
junction. In addition, sets of two mirrored avoided level crossings are observed at particular magnetic 
field values. b) Two-tone spectroscopy on one of these anticrossings shows a mode that depends hyper- 
bolically on the applied magnetic flux density. The constant line at f = 8.015 GHz can be attributed to 
the resonator, which naturally also changes its resonant frequency because of the self-Kerr effect when 
strongly driven by the second tone. 


so that a single vortex can tunnel in and out of the junction. Within the scope of the 
standard tunneling model (STM), this tunneling between two energetically similar 


states is described by the Hamiltonian |FLS64b}|Leg+87| 


Hsm = ; i = = 50: i wade. (6.1) 
Here, A denotes the coupling energy between the lowest bound states in each 
potential well, so to say the tunneling amplitude, £ is the energy asymmetry of the 
two wells, and ôy, ô- the Pauli matrices. The double well potential is symmetric in 
the sweet spots B,, where the |n — 1) and |n) fluxons states are energetically equal 
and by overcoming the energy at the junction boundary, one fluxon can tunnel in 
and out of the LJJ. At applied magnetic flux densities slightly below the sweet spot 
Ba < Bn, the |n — 1) state becomes more favorable, whereas the |n) state is more 
likely for slightly larger flux densities Ba = B,. The off-diagonal elements in the 
STM-Hamiltonian of Eq. couple the |n — 1) and |n) states through tunneling, 
so that the first two eigenstates of the two wells hybridize as 


|p) = cos() In — 1) sin($) In), (6.2a) 
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Figure 6.3: Operating principle of a Josephson vortex qubit. a) The |S21 |7 -spectrum against applied 
magnetic flux shows three pairs of avoided level crossings, whose magnetic flux densities B,, coincide 
with the simulated intersections of two adjacent fluxon states for a normalized junction length of 2 ~ 7 
c) (indicated by gray arrows). For these values, the states of |n — 1) and |n) fluxons are energetically 
similar, so a single fluxon can tunnel into the junction. In doing so, the fluxon has to overcome the 
energy of the junction boundary, which can be modeled as a double well potential c). In the sweet spot, 
where the two energies are identical (E, = E„-ı), the potential is symmetric (¢ = 0) and the fluxon is 
maximally delocalized between the two minima, so to say in the junction inside and outside. Off the 
sweet spot |e] > A the fluxon is trapped in one of the two minima. d) The energy difference of the two 
hybridized eigenstates |) shows a hyperbola as a function of the asymmetry e and can be excited as 
shown by two-tone spectroscopy (see Fig. 6.2)). 


py) = sin( 2) |n —1)4 cos( 5) In), (6.2b) 


where the mixing angle 0 = arctan(A/e) is introduced. In the basis of these 
eigenstates, the STM-Hamiltonian reads Hstm = TE ö, with the energy difference 
between the eigenstates 


E = V4 +2. (6.3) 


The energy of Eq. exactly represents the previous mentioned hyperbola, 
measured as a function of the applied magnetic flux, and thus the asymmetry of 


the two potential wells (see Fig.|6.2 b)|and Fig. |6.3 d)). 
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6.4 Coherence properties of vortex tunneling 


Knowing the transition energy between the ground and excited state from two- 
tone spectroscopy measurements, the quantum system can be driven with pulsed 
microwave signals to explore its coherence properties. For these measurements, a 
magnetic field bias at the sweet spot is suitable since the influence of field noise 
on the qubit frequency is minimal there. For the following pulsed microwave 
measurements, a magnetic flux density of Ba ~ By ~ 48.44 mT is applied, which 
cannot be guaranteed to permanently correspond to the exact sweet spot due to 
temporal instability of the magnetic bias over the measuring period of several days. 


6.4.1 Rabi oscillations 


In the experiment, the system is repeatedly excited from the ground state by 
x-pulses of variable duration T, meaning rotated around the x-axis to different 
angles 0, and read out immediately afterwards. Including decoherence and energy 
relaxation, the projection of this rotation with constant Rabi frequency Qpr on the 
z-axis results in a damped cosinusoidal oscillation 

(62) = e7*/TRabi cos (=). (6.4) 

2t x 

Figure |6.4a)| depicts such Rabi oscillations, which are characterized by alternating 
excitations due to photon absorption and deexcitations due to stimulated emission. 
From these measured data, on the one hand, the time constant of the envelope 
Trabi = 2/ (TI Ti T, 1) can be determined, which is given by the relaxation time 
Tı and dephasing time T. On the other hand, the half-period duration tz can be 
determined from the Rabi frequency, which is as the pulse length to excite the 
qubit from |0) to |1), i.e., to perform a 7t-rotation on the Bloch sphere (compare to 


Sec.2.3.4). 


The quantum Rabi frequency depends both on the drive frequency and amplitude, 
as explained in Sec. [2.3.4] For Rabi measurements with different drive frequency fa, 
Chevron patterns are expected according to Eq. (see Fig. [2.9 a)|for theoretical 
simulations). Measured data, as exemplarily depicted in Fig. are in good 
accordance to theoretical predictions since the Rabi frequency speeds up while the 
amplitude decreases with increasing detuning. These effects can be explained by 
the fact that the driven rotation is no longer around the pure x-axis, but around a 


tilted axis in the xz-plane [SN17 Sch20]. As a consequence, the rotations 


still start at |0), but |1) cannot be reached anymore due to the tilt. 
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Figure 6.4: Quantum mechanical Rabi oscillations in a Josephson vortex qubit. a) Repeatedly applied 
manipulation pulses of variable duration T rotate the qubit state by different azimuth angles on the Bloch 
sphere, so that the expectation value (6,) features Rabi oscillations. b) The Rabi frequency depends on 
the drive frequency or rather the detuning from the qubit frequency. Eq. reproduces the measured 
data and exhibits the typical Chevron pattern. c) For quantum mechanical Rabi oscillations, their 
frequency increases linearly with the amplitude of the drive tone. The good agreement of measured 
data and the theory proves the quantum nature of the Josephson vortex and confirms the device to be a 
Josephson vortex qubit. 


In Fig. 6.4c)|the quantum Rabi oscillation’s linear dependence on the drive ampli- 
tude V is shown, which is in excellent agreement with the theoretical expectations. 
This power dependence is characteristic of quantum mechanical Rabi oscillations, 
whose bare Rabi frequency without detuning equals the drive amplitude V. In 
other words, at higher powers, there are, on average more photons in the system, 
which drive the qubit more strongly, thus faster. The linear slope of the generalized 
Rabi frequencies Qpr against the drive amplitude V finally proves the Joseph- 
son vortex’ quantum character. Together with the hyperbolic spectrum, shown in 
Fig. these data manifest that a Josephson vortex tunnels coherently in the 
system under consideration, which makes the device a Josephson vortex qubit. 
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6.4.2 Relaxation and dephasing times 


The envelope of Rabi oscillations shows an exponential decay caused by interactions 
with the environment. This coupling to the environment is indispensable since 
qubits are intended to be manipulated and read out from outside. Two processes 
can be fundamentally distinguished: 


First, real qubits relax from the excited state to the ground state under spontaneous 
emission into the resonator, the environment, or internal modes in the solids such 
as defects. The time constant of this exponential decay is usually called energy 
relaxation time T4. In the Bloch sphere illustration, the energy relaxation manifests 
as a decay along the z-axis into |0), so that the qubit state is no longer bound to 
the sphere’s surface but may be in its interior. 


Second, fluctuations in the transition frequency wo; lead to the dephasing of the 
qubit. Averaging multiple measurement values with random dephasing reveals a 
decay towards zero within the xy-plane. For example, magnetic field noise leads to 
a fluctuating phase, and thus to another exponential decay with the characteristic 
dephasing time Tp. 


Relaxation time 


To measure the energy relaxation time T}, the qubit is brought into the excited state 
|1) by means of a t„-pulse. After a delay time At, the probability of the excited 
state (ö,) is determined using the resonator’s dispersive shift. During the time At 
the qubit is free to evolve and decays to the ground state, which is described by 


(02) =1—2e 4/0, (6.5) 


Figure [6.5 a)| depicts the schematic measurement protocol, the illustration in the 
Boch sphere, and the measured (6) data. By fitting these data to Eq. (6.5), the 
energy relaxation time T; can be determined, whose statistical distribution is shown 
in Fig. [6.5 d)| For this purpose, the T, measurement is repeated 1465 times within 
~ 86h and yields on average Tı = (45 +5) ps. 


Ramsey dephasing time 


To determine the dephasing time Tz, a “Ramsey” measurement is performed, 
where the qubit state is brought with a t/2-pulse around the x-axis to the state 
A (|0) + |1)) on the equator. Arriving in the xy-plane, the qubit evolves freely for 
a variable duration At, which is nothing but an effective Lamor precession at the 
detuning frequency Aw in the rotating frame. The quantum state can only be read 
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Figure 6.5: Energy relaxing and dephasing time measurements on a Josephson vortex qubit. a) To 
determine the energy relaxation time Tı experimentally, the qubit is excited to |1) with a t-pulse and 
after a waiting time Af the qubit state is read out. During the free evolution Af the qubit decays expo- 
nentially along the z-axis, where the half-life is the energy relaxation time T4. b) Since dephasing occurs 
perpendicular to z, but only the projection on z is measurable, in Ramsey dephasing measurements 
the qubit is rotated to and from the equator with a t/2-pulse at the beginning and just before the 
measurement, respectively. In between, the qubit evolves freely in the xy-plane for time At. Thereby 
detuning leads to an effective Lamor precession and dephasing to a decay towards the origin of the 
Bloch sphere, which in combination results in a damped oscillation. c) Low-frequency dephasing can be 
circumvented in Hahn-echo measurements by a t-pulse around y after half of the free evolution. The 
individual pulses are illustrated on the Bloch sphere in the same color. d) Long-term measurements for 
= 86h of Tj, Tr and I yield Gaussian distributions. 
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out by projection on the z-axis so that a second t,,/2-pulse around the x-axis is 
necessary. Without dephasing (as it is the case for At = 0), this pulse sequence 
effectively corresponds to a t,-pulse, split in two parts, so that the qubit ends 
up in |1). With dephasing, however, the state rotates around the z-axis, and thus 
cosinusoidal in y, which is rotated by the second t,/2-pulse in the z-direction|!|In 
addition to the effective Lamor precession due to detuning, dephasing causes the 
qubit to decay towards the origin of the Bloch sphere. In consequence, the overall 
result is a damped oscillation 


(êz) = -—e AT cos(AwAt) (6.6) 


with half-life TÈ. Since for very small detunings Aw, the first quarter wave of the 
oscillation would be difficult to distinguish from the superimposed exponential 
decay, and thus challenging to fit, one usually deliberately chooses a detuning 
Aw ~ 1MHz to see a few oscillations during the decay. The Ramsey pulse sequence, 
its illustration on the Bloch sphere, and example measurement data are depicted 
in Fig. Long-term measurements for ~ 86h with 568 repetitions give the 
distribution shown in Fig. and a mean Ramsey dephasing time of TR = 
(0.9 =e 0.4) HS. 


Hahn-echo dephasing time 


To become insensitive to dephasing with slower frequencies than the inverse 
measurement duration, a tz-pulse around y can be introduced at half of the free 
evolution time, i.e., after At/2 (see Fig. [6.5 c)). This pulse mirrors the state on the y- 
axis and the two rotations around the z-axis during At/2 partially compensate each 
other. In these so-called “Hahn-echo” or “spin-echo” measurements, an exponential 
decay 


yagi (6.7) 
is obtained with the characteristic Hahn-echo dephasing time TË > TÈ (see 


Fig. 6.5)). The significantly longer Hahn-echo dephasing time TF = (7.2 + 1.4) ps 
(for 1441 measurement repetitions) indicates that the Ramsey dephasing time 


strongly suffers from low-frequency noise of the qubit transition frequency. One 
major noise source is current noise of the magnet coil. Persistent current coils, for 
example, could provide a remedy for this. 


1 This is why, the measured (ô+) is plotted in Fig.|6.5b)| although the qubit decays in the xy-plane. 
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6.5 Concluding remarks 


In summary, we have demonstrated the, to our knowledge, first experimental real- 
ization of a Josephson vortex qubit using a quasi-one-dimensional high-impedance 
long junction. At particular applied magnetic flux densities, the device features 
coherent quantum tunneling of a single vortex in and out of the junction, form- 
ing a two-level system. Such quantum systems have a characteristic hyperbolic 
mode spectrum, observed experimentally by two-tone spectroscopy and strongly 
indicates the increased quantumness of the inherent Josephson vortices. 


Moreover, pulsed microwave measurements reveal quantum Rabi oscillations, 
whose frequency increases proportionally to the applied drive frequency. This ob- 
servation again proves the Josephson vortex’ quantum nature. The average energy 
relaxation time of T} ~ 451s is comparable to state-of-the-art superconducting 
qubits. This result is impressive against the background that the measured qubit 
is the very first of its kind, focusing on the principle implementation but not on 
optimized relaxation and coherence times. With that said, the average Hahn-echo 
dephasing time of TE ~ 7 ps is appealing, too. 


T should be theoretically limited by 2T,, which is why much longer T, times can 
be expected with an optimized setup regarding magnetic field noise. Apart from 
these limitations, which could be similarly solved for other superconducting qubit 
types, the vortex qubits devised in this work are a promising development for the 
research field of superconducting quantum circuits and the long-standing dream 
of quantum computing. 
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The goal of this work was to experimentally investigate the quantum coherence of 
fluxons in quasi-one-dimensional long Josephson junctions. Entering the Josephson 
vortices’ quantum regime is experimentally challenging to reach since 
such excitations in conventional junctions occur on an energy scale that exceeds 
the system’s quantized eigenenergies by orders of magnitude. In consequence, the 
Josephson medium and its inherent excitations are pretty classical. 


The long junction’s crucial energy scale decreases with its increasing impedance 
A= {L7 cw?) and is defined by the electrodes’ sheet inductance L-, the 
specific junction capacitance cy, and its width w. The approach followed in this work 


was to replace bulk junction electrodes with thin films of high-kinetic inductance 
superconductors in addition to very narrow junctions. The kinetic inductance in 
strips of the superconductor granular aluminum oxide exceeds the geometrical 
part by several orders of magnitude. It thus facilitates an increased quantumness 
of the long junction, as outlined in Sec. 2.4.4] 


Within the scope of this work, we developed fabrication techniques for impedance- 
tailored junctions featuring good controllability of the kinetic inductance during 
sputter deposition and, if desired, high impedances. 


With transport characterization measurements at millikelvin temperatures, we 
demonstrate high quality factors of the junctions. Long junctions also show typical 
magnetic diffraction patterns of their critical currents and zero-field and Fiske 
steps, which we use to determine the Josephson length and the Swihart velocity, 
respectively. These zero-field steps prove the existence of Josephson vortex exci- 
tations in the high-impedance long junctions. Furthermore, we determined the 
Josephson plasma frequency by switching current measurements under microwave 
irradiation. 


The results of these three independent measurements demonstrate that the sine- 
Gordon model is valid without restrictions for high-impedance long junctions 
and © = Aywp is still satisfied. Here, Ay and © are reduced by about one order of 
magnitude due to the additional kinetic inductance, while wp itis nearly unaffected. 
The reduced Swihart velocity is accompanied by an impedance increase up to 140, 
which corresponds to the largest value ever measured in long junctions to our 
knowledge [wil+22]. From these measurement results, we conclude that the chosen 
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approach is suitable to increase Josephson vortices’ quantumness so that they could 
behave like quantum particles. 


A suitable measurement scheme for quantum coherence properties is possible in 
the framework of circuit quantum electrodynamics, namely a Josephson medium 
connected to a microwave resonator. For this purpose, we galvanically couple high- 
impedance long junctions to a distributed microstrip resonator each, which we 
manipulate and read out by microwave photons. By measuring the anharmonicity 
and the magnetic field dependence of the coupled system’s resonant frequency, 
we attribute the junction’s influence on the resonator unambiguously. This result 
confirms the coupling of the two systems and is understood by the theoretical 
treatment of phase-biased junctions. 


We combine our findings gathered from previous experiments and focus on the in- 
trinsic quantumness of Josephson vortices by significantly increasing the impedance 
of the long junction through high kinetic electrode inductances and narrow junction 
widths. This junction is embedded in a microwave resonator for manipulation and 
readout with single photons. In microwave spectroscopy measurements versus ap- 
plied magnetic flux density, we observe hyperbolic modes at particular fields. This 
mode spectrum is characteristic of quantum tunneling in a two-level system. We 
attribute this signature to a tunneling process of a single fluxon in and out of the 
long junction since measured magnetic field values correspond to the intersections 
between |n — 1) and |n) fluxon states from theoretical predictions. 


By driving the transition of this two-level system by pulsed microwave signals, we 
observe quantum Rabi oscillations, whose frequency is proportional to the drive 
amplitude. This measurement demonstrates the quantum nature of the vortices 
for the first time to our knowledge. The vortex qubit features a measured mean 


energy relaxation time of Ti = (45 + 5) us and an average Hahn-echo dephasing 
time of TH = (7.2 + 1.4) ps. 


The experimental methods developed in this work could also be used in various 
future experiments. For instance, the coupling of the Josephson junction to a 
microwave resonator suggests using these circuits to study the internal mode 
spectrum of the long junction and to directly excite, e.g., the Josephson plasma 
frequency, zero-field, and Fiske resonances. 


Besides, the high-kinetic inductance electrodes enable impedance-tailored junctions 
to match external circuits and 50 © cables. This impedance matching is crucial for 
increasing the efficiency of Josephson flux-flow oscillators applied for microwave 
generation and amplification. 
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The discovery of Josephson quantum vortices is a promising starting point for 
further experimental and theoretical work. One could study the dynamic properties 
of such quantum particles in the nonlinear Josephson medium, i.e., quantum 
solitons within the scope of the quantum sine-Gordon model [Roy+21]. 


Beyond these uniform media, one could design more sophisticated junction geome- 
tries, such as spatially periodic structures, including potential barriers and wells. 
They can be realized by broadenings (“microshorts”) or constrictions (“microre- 
sistors”) of the long junction, where the quantum vortices are repelled or pinned. 
Quantum coherent entities in such metamaterials should lead to a band structure 
for propagating collective linear and nonlinear excitations, similar to phonons in 
crystals. 


The last point, we want to outline here, is optimizing the properties of the Joseph- 
son vortex qubit. The priority should be to reduce its dephasing by filtering the 
magnetic field noise or, ideally, using persistent current coils. Furthermore, alter- 
native vortex qubits are also possible. For instance, microshort- and mircoresistor- 
impurities can be used to design two-level systems whose pinning and repulsion 
potentials allow for coherent tunneling. 


111 


Bibliography 


[ACC66] 


[Aff81] 


[AP95] 


[AB63] 


[AR63] 


[ASB89] 


[Ann+10] 


[Ann10] 


[ADK18] 


[Bar50] 


B. Abeles, R. W. Cohen, and G. W. Cullen: Enhancement of Supercon- 
ductivity in Metal Films. Physical Review Letters 17 (1966), 632-634. 


por:/10.1103/PhysRevLett . 17 .632|(cit. on p.[56}. 


I. Affleck: Quantum-Statistical Metastability. Physical Review Letters 46 


(1981), 388-391. por:|10.1103/PhysRevLett . 46.388] cit. on p.{75). 


G. L. Alfimov and A. F. Popkov: Magnetic vortices in a distributed 
Josephson junction with electrodes of finite thickness. Physical Review B 


52 (1995), 4503-4510. vor. [10.1103/PhysRevB .52 .4503|(cit. on p. [40). 


V. Ambegaokar and A. Baratoff: Tunneling Between Superconduc- 
tors. Physical Review Letters 10 (1963), 486-489. Dor: |10 . 1103 / 


PhysRevLett 10.486) (cit. on pp. B1]78]79) 
P. W. Anderson and J. M. Rowell: Probable Observation of the Josephson 
Superconducting Tunneling Effect. Physical Review Letters 10 (1963), 


230-232. nox 10. 1108/PhyskevLett 10.230] (cit. on pp BO) 
S. M. Anlage, H. J. Snortland, and M. R. Beasley: A current controlled 
variable delay superconducting transmission line. IEEE Transactions on 
Magnetics 25 (1989), 1388-1391. por:|10 . 1109/20. 92554 | (cit. on pp.|16] 
17). 

A. J. Annunziata, D. F. Santavicca, L. Frunzio, G. Catelani, M. J. Rooks, 
A. Frydman, and D. E. Prober: Tunable superconducting nanoinductors. 


Nanotechnology 21 (2010), 445202. DOT: |10 . 1088/0957 -4484/21/44/ 
[445202] (cit. on p- 6). 


A. J. Annunziata: Single-Photon Detection, Kinetic Inductance, and Non- 
Equilibrium Dynamics in Niobium and Niobium Nitride Superconduct- 
ing Nanowires. PhD thesis. Yale University, 2010. urL: [https : //| 
[Thesis. pdf|(cit. on p.|16). 

P. W. Atkins, J. De Paula, and J. Keeler: Physical chemistry. 11th ed. 
Oxford, New York: Oxford University Press, 2018. URL: 


academia. edu/51098021|(cit. on p.|62). 


J. Bardeen: Zero-Point Vibrations and Superconductivity. Physical Review 


79 (1950), 167-168. vor. [10.1103/PhysRev. 79. 167.3|(cit. on p.[12). 


113 


Bibliography 


[Bar+10] 


[BJV75] 


[BP82] 


[BL64] 


[BK03] 


[BK04] 


[Bes06] 


[BCG10] 


[Bla+04] 


[Blo46] 


[BTS58] 


114 


R. Barends, N. Vercruyssen, A. Endo, P. J. de Visser, T. Zijlstra, T. M. 
Klapwijk, and J. J. A. Baselmans: Reduced frequency noise in super- 
conducting resonators. Applied Physics Letters 97 (2010), 033507. Dot: 


10. 1063/1 . 3467052] (cit. on p.|17). 


A. Barone, W. J. Johnson, and R. Vaglio: Current flow in large Josephson 
junctions. Journal of Applied Physics 46 (1975), 3628-3632. por: [10.] 


[106371322089] (cit. on pp. 2) Z4) 25478). 


A. Barone and G. Paterno: Physics and Applications of the Josephson 


Effect. John Wiley & Sons, Ltd, 1982. por: [10 . 1002/352760278X|(cit. on 
pp- (31) 8267 4) 7) 76} 82136). 

C. P. Bean and J. D. Livingston: Surface Barrier in Type-II Supercon- 
ductors. Physical Review Letters 12 (1964), 14-16. DOI: 
FaysRevbett. 12.14) (cit. on p-|) 

K. H. Bennemann and J. B. Ketterson, eds.: The Physics of Superconduc- 
tors. Vol. 1. 2 vols. Berlin, Heidelberg: Springer, 2003. pot: 
[578-3-642-68675-9](cit. on p. [12 

K. H. Bennemann and J. B. Ketterson, eds.: The Physics of Superconduc- 
tors. Vol. 2. 2 vols. Berlin, Heidelberg: Springer, 2004. vor: 
978-3-822-38917-2(cit. on p12) 

M. Bestehorn: Hydrodynamik und Strukturbildung. Springer Berlin, Hei- 
delberg, 2006. 392 pp. por:[10. 1007/3-540-33797-0] (cit. on p.|1). 

J. A. Blackburn, M. Cirillo, and N. Grenbech-Jensen: On the classical 
model for microwave induced escape from a Josephson washboard potential. 


Physics Letters A 374 (2010), 2827-2830. D01:[10.1016/j.physleta.] 
2910.04.077)cit. on p.B0) 

A. Blais, R.-S. Huang, A. Wallraff, S. M. Girvin, and R. J. Schoelkopf: 
Cavity quantum electrodynamics for superconducting electrical circuits: 
An architecture for quantum computation. Physical Review A 69 (2004), 


062320. vor: [f0.1708/PhysRevA..69.062320| (cit. on pp. T8} 2985) 
F. Bloch: Nuclear Induction. Physical Review 70 (1946), 460-474. Dot: 
[10.4108/PhysRev.70.460|(cit. on p. 27}. 


N. N. Bogoljubov, V. V. Tolmachov, and D. V. Sirkov: A New Method in 
the Theory of Superconductivity. Fortschritte der Physik 6 (1958), 605- 


682. por: [10 . 1002/prop . 19580061102] (cit. on p./13). 


Bibliography 


[Bol+08] 


[Bor+20] 


[Bra18] 


[BH54] 


[BHL83] 


[CF91] 


[CFL71] 


[CHS71] 


[Chu+81a] 


[Chu+81b] 


A. T. Bollinger, R. C. Dinsmore, A. Rogachev, and A. Bezryadin: 
Determination of the Superconductor-Insulator Phase Diagram for One- 
Dimensional Wires. Physical Review Letters 101 (2008), 227003. pot: 


10.1103/PhysRevLett . 101.227003|(cit. on p.{17). 


K. Borisov, D. Rieger, P. Winkel, F. Henriques, F. Valenti, et al.: Super- 
conducting granular aluminum resonators resilient to magnetic fields up to 


1 Tesla. Applied Physics Letters 117 (2020), 120502. por:[10. 1063/5 .| 
poteoa] cit. on pp. B3) P3) 


J. Braumüller: Quantum simulation experiments with superconducting 


circuits. 2018. vor: (cit. on pp. 
102). 

W. Buckel and R. Hilsch: Einfluß der Kondensation bei tiefen Temperaturen 
auf den elektrischen Widerstand und die Supraleitung für verschiedene 
Metalle. Zeitschrift für Physik 138 (1954), 109-120. DOI: 


BF01337903| (cit. on p.|56). 


M. Büttiker, E. P. Harris, and R. Landauer: Thermal activation in ex- 
tremely underdamped Josephson-junction circuits. Physical Review B 28 


(1983), 1268-1275. vor:|10.1103/PhysRevB. 26 .1268|(cit. on p.|75). 


J. G. Caputo and N. Flytzanis: Kink-antikink collisions in sine-Gordon 
and $ \varphi*4 $ models: Problems in the variational approach. Physical 
Review A 44 (1991), 6219-6225. Dot: 
(cit. on pp. 47} |136). 

J. T. Chen, T. F. Finnegan, and D. N. Langenberg: Anomalous dc current 
singularities in Josephson tunnel junctions. Physica 55 (1971), 413-420. 


DOT: |10. 1016/0031-8914 (71) 90282-5] (cit. on pp. {1}|82). 


P. V. Christiansen, E. B. Hansen, and C. J. Sjöström: Negative self- 
inductance in superconducting thin wires and weak links. Journal of Low 


Temperature Physics 4 (1971), 349-389. pot: 
(cit. on p.|16). 

T. Chui, G. Deutscher, P. Lindenfeld, and W. L. McLean: Conduction in 
granular aluminum near the metal-insulator transition. Physical Review B 
23 (1981), 6172-6175. vor. |10. 1103/PhysRevB.23.6172|(cit. on p.56). 
T. Chui, P. Lindenfeld, W. L. McLean, and K. Mui: Coupling and isola- 
tion: Critical field and transition temperature of superconducting granular 


aluminum. Physical Review B 24 (1981), 6728-6731. vor: 
PiystevB. 26,676) (cit.on p.53} 


115 


Bibliography 


[Cir+97] 


[Cir+98] 


[Cla66] 


[CW08] 


[CA68] 


[Col15] 


[Czy17] 


[Dah+68] 


[Dau+21] 


[DP19] 


[DF61] 


116 


M. Cirillo, T. Doderer, S. G. Lachenmann, F. Santucci, and N. Grønbech- 
Jensen: Dynamical evidence of critical fields in Josephson junctions. Physical 
Review B 56 (1997), 11889-11896. por:|10.1103/PhysRevB.56. 11889] 
(cit. on p.|83). 

M. Cirillo, N. Grønbech-Jensen, M. R. Samuelsen, M. Salerno, and 
G. V. Rinati: Fiske modes and Eck steps in long Josephson junctions: Theory 
and experiments. Physical Review B 58 (1998), 12377-12384. vor: |10.] 


1103/PhysRevB.58.12377|(cit. on p. 83). 


J. Clarke: A superconducting galvanometer employing Josephson tun- 
nelling. The Philosophical Magazine: A Journal of Theoretical Ex- 


perimental and Applied Physics 13 (1966), 115-127. pot: 
Garesasssosarımi|(c. on p. B) 


J. Clarke and F. K. Wilhelm: Superconducting quantum bits. Nature 453 


(2008), 1031-1042. por: [10.1038 /nature07126] cit. on p.[30). 


R. W. Cohen and B. Abeles: Superconductivity in Granular Aluminum 


Films. Physical Review 168 (1968), 444-450. por: |10.1103/PhysRev.| 
68.424] (it. on pp. B58} 


P. Coleman: Introduction to Many-Body Physics. Cambridge: Cambridge 


University Press, 2015. Do1:]10.1017/CB09781139020916|(cit. on p.[13). 


G. Czycholl: Theoretische Festkörperphysik Band 2: Anwendungen: Nicht- 
gleichgewicht, Verhalten in äußeren Feldern, kollektive Phänomene. 2017. 
por: 10. 1007/878-8-663-53701-8|(cit. on pp.12][3) 

A. J. Dahm, A. Denenstein, T. F. Finnegan, D. N. Langenberg, and D. J. 
Scalapino: Study of the Josephson Plasma Resonance. Physical Review 
Letters 20 (1968), 859-863. por: [10.1103/PhysRevLett.20 859] (cit. on 
pp. B4] 79). 

H. Dausy, L. Nulens, B. Raes, M. J. Van Bael, and J. Van de Von- 
del: Impact of Kinetic Inductance on the Critical-Current Oscillations of 
Nanobridge SQUIDs. Physical Review Applied 16 (2021), 024013. por: 


10.1103/PhysRevApplied . 16.024013](cit. on p.|17). 


P. G. De Gennes and P. A. Pincus: Superconductivity of Metals and Alloys. 


Boca Raton: CRC Press, 2019. 292 pp. por:[10 . 1201/9780429497032| 
(cit. on pp. [11} 12}. 

B. S. Deaver and W. M. Fairbank: Experimental Evidence for Quantized 
Flux in Superconducting Cylinders. Physical Review Letters 7 (1961), 


43-46. por:(10.1103/PhysRevLett .7 .43|(cit. on pp. [1}[8}. 


Bibliography 


[DMO08] 


[Deu+73a] 


[Deu+73b] 


[DMC85] 


[DN61] 


[Due+81] 


[EGK15] 


[EH05] 


[EFP83] 


[F118] 


D. Depla and S. Mahieu, eds.: Reactive Sputter Deposition. Red. by R. 
Hull, R. M. Osgood, J. Parisi, and H. Warlimont. Vol. 109. Springer 
Series in Materials Science. Berlin, Heidelberg: Springer, 2008. Dot: 


19.10077978-3-5%0-76884-3|ct. on p. 7}. 

G. Deutscher, H. Fenichel, M. Gershenson, E. Griinbaum, and Z. 
Ovadyahu: Transition to zero dimensionality in granular aluminum su- 
perconducting films. Journal of Low Temperature Physics 10 (1973), 


231-243. por: [10.1007 /BF00655256|(cit. on p.[56). 


G. Deutscher, M. Gershenson, E. Grünbaum, and Y. Imry: Granular 
Superconducting Films. Journal of Vacuum Science and Technology 10 


(1973), 697-701. por:|10. 1116/1 .1318416|(cit. on p.56). 

M. H. Devoret, J. M. Martinis, and J. Clarke: Measurements of Macro- 
scopic Quantum Tunneling out of the Zero-Voltage State of a Current-Biased 
Josephson Junction. Physical Review Letters 55 (1985), 1908-1911. pot: 


[10.1108/PhysRevLett 65. 1908]{cit. on pp. 74} 73) 
R. Doll and M. Nabauer: Experimental Proof of Magnetic Flux Quantiza- 
tion in a Superconducting Ring. Physical Review Letters 7 (1961), 51-52. 


vor. 1108/PhysRevLett. 7.51] (cit. on pp. [}8). 
B. Dueholm, E. Joergensen, O. A. Levring, J. Mygind, N. F. Pedersen, 
M. R. Samuelsen, O. H. Olsen, and M. Cirillo: Dynamic fluxon model 
for Fiske steps in long Josephson junctions. Physica B+C 108 (1981), 1303- 


1305. por: (cit. on p. 3). 

E. G. Ekomasov, A. M. Gumerov, and R. V. Kudryavtsev: On the 
possibility of the observation of the resonance interaction between kinks of 
the sine-Gordon equation and localized waves in real physical systems. JETP 


Letters 101 (2015), 835-839. por:|10.1134/80021364016120061|(cit. on 
PER 


C. Enss and S. Hunklinger: Low-Temperature Physics. Berlin, Heidelberg: 


Springer, 2005. por: [10 .1007/b137878| (cit. on pp. |1} 62). 


S. N. Erné, A. Ferrigno, and R. D. Parmentier: Fluxon propagation and 
Fiske steps in long Josephson tunnel junctions. Physical Review B 27 


(1983), 5440-5446. vor:|10.1103/PhysRevB 27.5440] (cit. on p.|83). 


M. V. Feigel’man and L. B. Ioffe: Microwave Properties of Superconductors 
Close to the Superconductor-Insulator Transition. Physical Review Letters 


120 (2018), 037004. vor: [10 .1103/PhysRevLett . 120.037004] (cit. on 
p-P?}. 


117 


Bibliography 


[FP63] 


[FLS64a] 


[FLS64b] 


[FFT73] 


[Fis64] 


[FU03] 


[Fri44] 


[Frö50] 


[FD74a] 


[FD74b] 


[FD73] 


118 


R. A. Ferrell and R. E. Prange: Self-Field Limiting of Josephson Tunneling 
of Superconducting Electron Pairs. Physical Review Letters 10 (1963), 


479-481. Do1:|10.1103/PhysRevLett.10.479|(cit. on pp.|38}|39} |78). 
PP 


R. P. Feynman, R. B. Leighton, and M. L. Sands: The Feynman Lectures 
on Physics: Electromagnetism and Matter. Vol. 2. Addison-Wesley Pub. 


Co., 1964. URE: 
[Eoc Rem] (it. on p.B0) 

R. P. Feynman, R. B. Leighton, and M. L. Sands: The Feynman Lectures 
on Physics: Quantum mechanics. Vol. 3. Addison-Wesley Pub. Co., 1964. 


URE: 
(cit. on p.|100). 

B. F. Field, T. F. Finnegan, and J. Toots: Volt Maintenance at NBS via 
2e/h: A New Definition of the NBS Volt. Metrologia 9 (1973), 155-166. 


pot: [10 .1088/0026- 1394/9/4/003](cit. on p. 30). 


M. D. Fiske: Temperature and Magnetic Field Dependences of the Josephson 
Tunneling Current. Reviews of Modern Physics 36 (1964), 221-222. pot: 


10.1103/RevModPhys . 36. 221|(cit. on pp.|1}|83). 
PP 


M. V. Fistul and A. V. Ustinov: Josephson vortex interaction mediated 
by cavity modes: Tunable coupling for superconducting qubits. Physical 


Review B 68 (2003), 132509. por: |10.1108/PhysRevB.68. 132509|(cit. 
on p.|2). 

H. T. Friis: Noise Figures of Radio Receivers. Proceedings of the IRE 32 
(1944), 419-422, vor {0.11087 JRPROG. 1944, 252049] (cit. on p. BA) 


H. Frohlich: Theory of the Superconducting State. I. The Ground State at 
the Absolute Zero of Temperature. Physical Review 79 (1950), 845-856. 


por:{10.1103/PhysRev .79.845]|(cit. on p.{12). 


T. A. Fulton and L. N. Dunkleberger: Lifetime of the zero-voltage state in 
Josephson tunnel junctions. Physical Review B 9 (1974), 4760-4768. Dot: 


10.1103/PhysRevB.9.4760|(cit. on p.|73). 


T. A. Fulton and L. N. Dunkleberger: Vortex propagation and radiation 
emission in Josephson tunnel junctions. Rev.Phys.Appl. 9 (1974), 299-303. 


por:|10.1051/rphysap : 0197400901029900} (cit. on pp. |1}/82). 
pnoysap PP 


T. A. Fulton and R. C. Dynes: Single vortex propagation in Josephson 
tunnel junctions. Solid State Communications 12 (1973), 57-61. Dor: 


10. 1016/0038- 1098 (73) 90345- 1|(cit. on pp. [1}[82). 


Bibliography 


[GC08] 


[Geb22] 


[Geb+20a] 


[Geb+20b] 


[Gia+06] 


[GL09] 


[GT57] 


[GPS14] 


[Gor59] 


[GC34] 


J. Garrison and R. Chiao: Quantum Optics. Oxford Graduate Texts. 


Oxford: Oxford University Press, 2008. 730 pp. Dot: 
[5078750188508861 .001.0001] (cit. on p. 25), 


R. Gebauer: A Flexible FPGA-based Control Platform for Superconduct- 


ing Multi-Qubit Experiments. PhD thesis. 2022. DOI: 
[0001aT695](cit. on p.68) 

R. Gebauer, N. Karcher, D. Gusenkova, M. Spiecker, L. Grünhaupt, 
et al.: State preparation of a fluxonium qubit with feedback from a custom 
FPGA-based platform. AIP Conference Proceedings 2241 (2020), 020015. 


por: 0-108375 -00 11721] (cit. on p68) 
R. Gebauer, N. Karcher, J. Hurst, M. Weber, and O. Sander: Accelerating 
complex control schemes on a heterogeneous MPSoC platform for quantum 


computing (2020). Do1:|10.48550/arXiv.2004.07755|(cit. on p./68). 

F. Giazotto, T. T. Heikkilä, A. Luukanen, A. M. Savin, and J. P. Pekola: 
Opportunities for mesoscopics in thermometry and refrigeration: Physics 
and applications. Reviews of Modern Physics 78 (2006), 217-274. DOT: 


10.1103/RevModPhys.78.217|(cit. on p.|75). 


V. L. Ginzburg and L. D. Landau: “On the Theory of Superconductiv- 
ity”. On Superconductivity and Superfluidity: A Scientific Autobiography. 
Ed. by V. L. Ginzburg. Berlin, Heidelberg: Springer, 2009, pp. 113-137. 


DOI:|10. 1007/978-3-540-68008-6_4|(cit. on p. 8). 


R. E. Glover and M. Tinkham: Conductivity of Superconducting Films for 
Photon Energies between 0.3 and $ 40kT_c $. Physical Review 108 (1957), 


243-256. por: [10. 1108/PhysRev. 108.249] (cit. on p. 56) 
H. Goldstein, C. P. Poole, and J. L. Safko: Classical Mechanics. 3rd ed. 
Pearson, 2014. URL: [https : / / pearson . com/ en - us / subject - 
catalog/p/classical-mechanics/P200000006871/9780201657029 
(cit. on p. 48). 

L. P. Gor’Kov: Microscopic derivation of the Ginzburg-Landau equations 
in the theory of superconductivity. Soviet Physics JETP 36 (9) (1959), 
1364-1367. URL: |http://www.jetp.ras.ru/cgi-bin/dn/e_09_06_ 
E362:pa2|(cit. on p.[10} 

C. J. Gorter and H. Casimir: On supraconductivity I. Physica 1 (1934), 


306-320. por: (cit. on pp.B][7). 


119 


Bibliography 


[Gra+08] 


[Gre+04] 


[GMD16] 


[Grü+18] 


[Grü+19] 


[Gu+17] 


[Ham00] 


[HR06] 


[Haz+19] 


120 


G. J. Grabovskij, L. J. Swenson, O. Buisson, C. Hoffmann, A. Monfar- 
dini, and J.-C. Villégier: In situ measurement of the permittivity of helium 
using microwave NbN resonators. Applied Physics Letters 93 (2008), 


134102. por: 10. 1063/1. 2996263] (cit. on p.{I7). 

N. Gronbech-Jensen, M. G. Castellano, F. Chiarello, M. Cirillo, C. 
Cosmelli, L. V. Filippenko, R. Russo, and G. Torrioli: Microwave-Induced 
Thermal Escape in Josephson Junctions. Physical Review Letters 93 (2004), 


107002. por:[10.1103/PhysRevLett .93. 107002 (cit. on p.|80). 


R. Gross, A. Marx, and F. Deppe: Applied Superconductivity: Josephson 
Effect and Superconducting Electronics. De Gruyter Textbook Series. 
Walter De Gruyter Incorporated, 2016 (cit. on p.2). 


L. Grünhaupt, N. Maleeva, S. T. Skacel, M. Calvo, F. Levy-Bertrand, 
A. V. Ustinov, H. Rotzinger, A. Monfardini, G. Catelani, and I. M. 
Pop: Loss Mechanisms and Quasiparticle Dynamics in Superconducting 
Microwave Resonators Made of Thin-Film Granular Aluminum. Physical 


Review Letters 121 (2018), 117001. pot: |10.1103/PhysRevLett .121.| 
[117001] (cit. on p- 55). 

L. Grünhaupt, M. Spiecker, D. Gusenkova, N. Maleeva, S. T. Skacel, et 
al.: Granular aluminium as a superconducting material for high-impedance 
quantum circuits. Nature Materials 18 (2019), 816-819. pot: [10. 1038/] 


s41563-019-0350-3](cit. on p.|55). 
X. Gu, A. F. Kockum, A. Miranowicz, Y.-x. Liu, and F. Nori: Microwave 


photonics with superconducting quantum circuits. Physics Reports. Mi- 
crowave photonics with superconducting quantum circuits 718-719 


(2017), 1-102. por: [£0. 404675 .paysrep -2017-10-00 (cit. on p. 3) 


C. A. Hamilton: Josephson voltage standards. Review of Scientific In- 


struments 71 (2000), 3611-3623. pot: (cit. on 
p.B0). 

S. Haroche and J.-M. Raimond: Exploring the Quantum: Atoms, Cavities, 
and Photons. Oxford Graduate Texts. Oxford: Oxford University Press, 


2006. 616 pp. vor: |10. 10837acprof:080/9780198509121..001.0001] 
(cit. on p.[29). 

T. M. Hazard, A. Gyenis, A. Di Paolo, A. T. Asfaw, S. A. Lyon, A. Blais, 
and A. A. Houck: Nanowire Superinductance Fluxonium Qubit. Physical 


Review Letters 122 (2019), 010504. por:(10.1103/PhysRevLett .122.| 
[010504] (cit. on p.|17). 


Bibliography 


[HSB94] 


[1590] 


Dak+64] 


DC 63] 


[JV09] 


[Jos62] 


[Jos64] 


[Jos65] 


[Kam08] 


[Kam11] 


[K196] 


Z. Hermon, A. Stern, and E. Ben-Jacob: Quantum dynamics of a fluxon 
in a long circular Josephson junction. Physical Review B 49 (1994), 9757- 


9762. Do1:|10.1103/PhysRevB.49.9757|(cit. on p.[52). 


Y. M. Ivanchenko and T. K. Soboleva: Nonlocal interaction in Josephson 


junctions. Physics Letters A 147 (1990), 65-69. Dot: 
[801009900156] cit. on p. 0. 

R. C. Jaklevic, J. Lambe, A. H. Silver, and J. E. Mercereau: Quantum 
Interference Effects in Josephson Tunneling. Physical Review Letters 12 


(1964), 159-160. por: T0. 11037Physkeviett. 12. 189](cit. on p. B0) 


E. T. Jaynes and F. W. Cummings: Comparison of quantum and semiclas- 
sical radiation theories with application to the beam maser. Proceedings 


of the IEEE 51 (1963), 89-109. por: 10. 1109/PROC 1963.1664] (cit. on 
p- 29). 

D. Jeltsema and A. J. Van Der Schaft: Lagrangian and Hamiltonian 
formulation of transmission line systems with boundary energy flow. Reports 


on Mathematical Physics 63 (2009), 55-74. DOI: 
fS77COSYO0OGH-3| (cit. on p. B7) 


B. D. Josephson: Possible new effects in superconductive tunnelling. 


Physics Letters 1 (1962), 251-253. DOI: 
[81369-0|cit. on pp. 0) 

B. D. Josephson: Coupled Superconductors. Reviews of Modern Physics 
36 (1964), 216-220. vor: [10.1103/RevModPhys .36. 216|(cit. on pp. [I] 
32] [381] Bo) 

B. D. Josephson: Supercurrents through barriers. Advances in Physics 14 
(1965), 419-451. por: 10. 1080700018736500101081] (cit. on pp. I) B2 
47). 

H. Kamerlingh Onnes: The liquefaction of helium (1908), 168. URL: https :| 
/ / www . dwc . knaw . nl / DL / publications /PU00013525 . pdf| (cit. on 
p.|5). 

H. Kamerlingh Onnes: The resistance of pure mercury at helium temper- 
atures. Commun. Phys. Lab. Univ. Leiden, b 120 (1911). urr: [https] 


//physics.ucf.edu/~rep/EDII/Onnes1911.pdf}(cit. on p- B). 


T. Kato and M. Imada: Macroscopic Quantum Tunneling of a Fluxon in 
a Long Josephson Junction. Journal of the Physical Society of Japan 65 


(1996), 2963-2975. vor. |10. 1143/3P8J.68.2963](cit. on pp. 2} 5). 


121 


Bibliography 


[Kem06] 


[KWUO02] 


[Kha+12] 


[KD11] 


[KDP06] 


[Kit05] 


[Kw19] 


[KMM95] 


[Kos+01] 


122 


A. Kemp: Quantum and thermal phase escape in extended Josephson systems 


(2006). URL:|https://www.osti. gov/etdeweb/biblio/20865471| (cit. 
on p.|53). 

A. Kemp, A. Wallraff, and A. V. Ustinov: Josephson Vortex Qubit: Design, 
Preparation and Read-Out. physica status solidi (b) 233 (2002), 472-481. 
DOI: 
[6052-3] (cit. on pp. 253) 

M. S. Khalil, M. J. A. Stoutimore, F. C. Wellstood, and K. D. Osborn: 
An analysis method for asymmetric resonator transmission applied to super- 
conducting devices. Journal of Applied Physics 111 (2012), 054510. pot: 


10. 1063/1 . 3692073] (cit. on p. 22). 


J. H. Kim and R. P. Dhungana: Effects of a resonant cavity on macroscopic 
quantum tunneling of fluxons in long Josephson junctions. Physical Review 


B 83 (2011), 064503. por: 10. 1103/PhysRevB .83. 064503|(cit. on p.P). 
J. H. Kim, R. P. Dhungana, and K.-S. Park: Decoherence in Josephson 
vortex quantum bits: Long-Josephson-junction approach to a two-state sys- 


tem. Physical Review B 73 (2006), 214506. DOI: 
[F3-214506](cit. on pP). 


C. Kittel: Introduction to Solid State Physics. 8th ed. New York: Wiley, 


2005. URL: [https : //www. wiley .com/en-us/Introduction+to+ 
Solid+State+Physics+8th+Edition-p-9780471415268J|(cit. on pp.|5} 
14). 


R. Kleiner and H. Wang: “Intrinsic Josephson Junctions in High Tem- 
perature Superconductors”. Fundamentals and Frontiers of the Josephson 
Effect. Ed. by F. Tafuri. Springer Series in Materials Science. Cham: 
Springer International Publishing, 2019, pp. 367—454. DOI: 
[578-3-080-20726-7_10|(cit. on pp. [15] 87) 

A. W. Kleinsasser, R. E. Miller, and W. H. Mallison: Dependence of 
critical current density on oxygen exposure in Nb-AlO/sub x/-Nb tunnel 
junctions. IEEE Transactions on Applied Superconductivity 5 (1995), 
26-30. DOT: [10 . 1109/77 . 384565] (cit. on pp. 59} (141). 

V. P. Koshelets, P. N. Dmitriev, A. B. Ermakov, A. S. Sobolev, A. M. 
Baryshev, P. R. Wesselius, and J. Mygind: Radiation linewidth of flux-flow 
oscillators. Superconductor Science and Technology 14 (2001), 1040- 


1043. Do1:|10.1088/0953-2048/14/12/312]|(cit. on p.[1p. 


Bibliography 


[Kra+19] 


[Kri+19] 


[Kri22] 
[Kul65] 


[Kul67] 


[Lan37] 


[LG50] 


[Lau+95] 


[LS67] 


[Led+10] 


[Leg+87] 


[Leg87] 


P. Krantz, M. Kjaergaard, F. Yan, T. P. Orlando, S. Gustavsson, and 
W. D. Oliver: A quantum engineer’s guide to superconducting qubits. 
Applied Physics Reviews 6 (2019), 021318. Dot: 
(cit. on pp. [2] 27). 

S. Krinner, S. Storz, P. Kurpiers, P. Magnard, J. Heinsoo, R. Keller, J. 
Lütolf, C. Eichler, and A. Wallraff: Engineering cryogenic setups for 100- 
qubit scale superconducting circuit systems. EPJ Quantum Technology 6 


(2019), 1-29. pot-f0.1140/epjqt/540507-018-0072-O|(cit. on p. (6) 


M. Kristen. In preparation (2022) (cit. on p.|97). 


I. O. Kulik: Theory of “Steps” of Voltage-current Characteristic of the 
Josephson Tunnel Current. ZhETF Pisma Redaktsiiu 2 (1965), 134. URL: 


http: //jetpletters.ru/ps/1597/article_ 24482. shtml (cit. on 
PP. 83) 

I. O. Kulik: Theory of the tunnel effect in a superconductor. Soviet Physics 
Journal 10 (1967), 12-14. D01:|10.1007/BF00819974|(cit. on pp. 1183). 


L. D. Landau: On the theory of phase transitions. I. Phys. Z. Sowjet. 11 


(1937). por: 0. 101678878-0-08-010688-4.50034-1|(it. on p. B) 
L. D. Landau and V. L. Ginzburg: On the theory of superconductivity. Zh. 
Eksp. Teor. Fiz.; (USSR) 20:12 (1950), 1064-1082. Do1:|10.1016/B978-| 
0-08-010586-4.50035-3|(cit. on p-B). 

A. Laub, T. Doderer, S. G. Lachenmann, R. P. Huebener, and V. A. 
Oboznov: Lorentz Contraction of Flux Quanta Observed in Experiments 
with Annular Josephson Tunnel Junctions. Physical Review Letters 75 


(1995), 1372-1375. por:|i0,1103/PhysRevLett.. 75. 1372|(cit. on p.). 


P. Lebwohl and M. J. Stephen: Properties of Vortex Lines in Supercon- 
ducting Barriers. Physical Review 163 (1967), 376-379. pot:{10. 1103/| 
Phyakev. 163.379] (cit. on p. 59). 

H. G. Leduc, B. Bumble, P. K. Day, B. H. Eom, J. Gao, et al.: Titanium 
nitride films for ultrasensitive microresonator detectors. Applied Physics 
Letters 97 (2010), 102509. DOT: [10 . 1063/1. 3480420] (cit. on p.[17). 

A. J. Leggett, S. Chakravarty, A. T. Dorsey, M. P. A. Fisher, A. Garg, 
and W. Zwerger: Dynamics of the dissipative two-state system. Reviews 
of Modern Physics 59 (1987), 1-85. por: [10.1103/RevModPhys . 59. 1] 
(cit. on p.|100). 

O. Legrand: Kink-antikink dissociation and annihilation: A collective- 
coordinate description. Physical Review A 36 (1987), 5068-5073. Dot: 


[10.1108/PhysRevA. 38.8083] (cit. on pp. 47 136), 


123 


Bibliography 


[LPS83] 


[Lik86] 


[Lon48] 


[Lon50] 


[LL35] 


[LU08] 


[Luo+14] 


[Mae+95] 


[MW 18] 


[MK83] 


[MB58] 


124 


O. A. Levring, N. F. Pedersen, and M. R. Samuelsen: Perturbation 
calculation of magnetic field dependence of fluxon dynamics in long inline 
and overlap Josephson junctions. Journal of Applied Physics 54 (1983), 


987-991. por:|i0, 1063/1 .332025|(cit. on p.|48). 

K. K. Likharev: Dynamics of Josephson junctions and circuits. Philadel- 
phia, Pa.: Gordon and Breach Science Publishers, 1986. URL: 
//www.taylorfrancis.com/books/mono/10.1201/9781315141572/ 


(cit. on pp. [1} 89} [81] 82). 


F. London: On the Problem of the Molecular Theory of Superconductivity. 


Physical Review 74 (1948), 562-573. DOI: 
(cit. on pp. [7} [10} [13). 


F. London: Macroscopic Theory of Superconductivty. John Wiley & Sons. 


Superfluids 1 (1950). URL: http: //w2agz.com/Library/Classic% 
20Papers%20in%20Superconductivity/London , 420Superf luids 7% 
20 (Wiley, %20New%20York, 4201950) . pdf] (cit. on pp. {1} 8). 


F. London and H. London: Supraleitung und Diamagnetismus. Physica 2 


(1935), 341-354. vor: (cit. on p.f). 


A. Lukashenko and A. V. Ustinov: Improved powder filters for qubit 
measurements. Review of Scientific Instruments 79 (2008), 014701. DOI: 


10. 1063/1 . 2827515] (cit. on p.|65). 


J. Luomahaara, V. Vesterinen, L. Grönberg, and J. Hassel: Kinetic in- 
ductance magnetometer. Nature Communications 5 (2014), 4872. DOT: 


10. 1038/ncomms5872|(cit. on p.|17). 


M. Maezawa, M. Aoyagi, H. Nakagawa, I. Kurosawa, and S. Takada: 
Specific capacitance of Nb/AlOx/Nb Josephson junctions with critical current 
densities in the range of 0.1-18 kA/cm?. Applied Physics Letters 66 (1995), 


2134-2136. vor. [10.1063/1 113927] (cit. on p. BT). 


W. Martienssen and H. Warlimont: Springer handbook of condensed mat- 
ter and materials data. 2nd ed. Springer Handbooks. Cham: Springer 


International Publishing, 2018. pot: 
(cit. on p. 11). 

A. Matsuda and T. Kawakami: Fluxon Propagation on a Josephson Trans- 
mission Line. Physical Review Letters 51 (1983), 694-697. DOT: 


FaysRevLett 61.694] (it. on p.B7) 
D. C. Mattis and J. Bardeen: Theory of the Anomalous Skin Effect in 
Normal and Superconducting Metals. Physical Review 111 (1958), 412- 


417. vot [10.1108/PhysRev 111 412|cit. on pp-(18]56), 


Bibliography 


[Max50] 


[MSJ69] 


[McC68] 


[MS78] 


[MO33] 


[Mer70] 


[MT69] 


[Mon+13] 


[Nag+83] 


[Nag+85] 


[NPT99] 


E. Maxwell: Isotope Effect in the Superconductivity of Mercury. Physical 


Review 78 (1950), 477-477. vor: 10. 1103/PhysRev. 78.477] (cit. on 
p- 53. 

A. F. Mayadas, M. Shatzkes, and J. F. Janak: Electrical resistivity model 
for polycrystalline films: the case of specular reflection at external surfaces. 
Applied Physics Letters 14 (1969), 345-347. por: 


(cit. on pp. [58}[59)}. 
D. E. McCumber: Effect of ac Impedance on dc Voltage-Current Character- 
istics of Superconductor Weak-Link Junctions. Journal of Applied Physics 


39 (1968), 3113-3118. por: 10.106371. 1656743] (cit. on pp. B2} 33). 


D. W. McLaughlin and A. C. Scott: Perturbation analysis of fluxon dynam- 


ics. Physical Review A 18 (1978), 1652-1680. vor: 
18.1852](cit. on pp. [1} 8) 51) 52] 82). 

W. Meissner and R. Ochsenfeld: Ein neuer Effekt bei Eintritt der Supraleit- 
fähigkeit. Naturwissenschaften 21 (1933), 787-788. DOI: 
[5F01504259| (cit. on p.) 


J. E. Mercereau: Superconducting magnetometers. Revue de Physique 


Appliquée 5 (1970), 13-20. por: 10. 1051/physap:019700050101300] 
(cit. on p. BO). 

R. Meservey and P. M. Tedrow: Measurements of the Kinetic Inductance of 
Superconducting Linear Structures. Journal of Applied Physics 40 (1969), 


2028-2034. por: |10.1063/1. 1657905|(cit. on p.[T6}. 
R. Monaco, V. P. Koshelets, A. Mukhortova, and J. Mygind: Self-field 
effects in window-type Josephson tunnel junctions. Superconductor Science 


and Technology 26 (2013), 055021. Do1:|10. 1088/0953- 2048/26/5/ 
[55023] cit. on pp. 28). 


T. Nagatsuma, K. Enpuku, F. Irie, and K. Yoshida: Flux-flow type Joseph- 
son oscillator for millimeter and submillimeter wave region. Journal of Ap- 


plied Physics 54 (1983), 3302-3309. pot: [10 . 1063/1 . 332443} (cit. on 
p. [). 

T. Nagatsuma, K. Enpuku, H. Iwakura, and K. Yoshida: Flux-Flow- 
Type Josephson Linear Amplifier with Large Gain and Wide Linear Range. 
Japanese Journal of Applied Physics 24 (1985), L599. vor: 
[aP -27-1599](cit. on p.f} 

Y. Nakamura, Y. A. Pashkin, and J. S. Tsai: Coherent control of macroscopic 
quantum states in a single-Cooper-pair box. Nature 398 (1999), 786-788. 


DOr: |10. 1038/19718] (cit. on p. BO}. 


125 


Bibliography 


[NSS60] 


[NBB19] 


[Nol14] 


[Nol18] 


[Nor95] 


[Ols+86] 


[OD91] 


[OS67] 


[PRS91] 


[Par78] 


[Ped83] 


126 


J. Nicol, S. Shapiro, and P. H. Smith: Direct Measurement of the Super- 
conducting Energy Gap. Physical Review Letters 5 (1960), 461-464. DOI: 


[10.1108/PhysRevbett 5. 461|(cit. on pp.13]70) 


D. Niepce, J. Burnett, and J. Bylander: High Kinetic Inductance NbN 
Nanowire Superinductors. Physical Review Applied 11 (2019), 044014. 


D01:[10.1103/PhysRevApplied.11.044014|(cit. on p.|17). 


W. Nolting: Grundkurs Theoretische Physik 2. Springer-Lehrbuch. Berlin, 


Heidelberg: Springer, 2014. por:[10.1007/978-3-642-41980-5|(cit. on 
p.#8). 

W. Nolting: Grundkurs Theoretische Physik 1. Ed. by W. Nolting. Berlin, 
Heidelberg: Springer, 2018. por:[10.1007/978-3-662-57584-0|(cit. on 
p. (25). 

J. E. Nordman: Superconductive amplifying devices using fluxon dynamics. 
Superconductor Science and Technology 8 (1995), 681-699. Dot: [o] 


1088/0953-2048/8/9/001|(cit. on p.|1). 
O. H. Olsen, N. F. Pedersen, M. R. Samuelsen, H. Svensmark, and D. 


Welner: Perturbation treatment of boundary conditions for fluxon motion 
in long Josephson junctions. Physical Review B 33 (1986), 168-173. Dot: 


10.1103/PhysRevB. 33.168] (cit. on p.|48). 


T. P. Orlando and K. A. Delin: Foundations of applied superconductivity. 
Addison-Wesley Pub. Co., 1991 (cit. on pp. 32} 50}. 


C. S. Owen and D. J. Scalapino: Vortex Structure and Critical Currents 
in Josephson Junctions. Physical Review 164 (1967), 538-544. DOI: [10.] 


11037PhysRev.184.538](cit. on pp. B9} P22). 


S. Pagano, B. Ruggiero, and E. Sarnelli: Magnetic-field dependence of the 
critical current in long Josephson junctions. Physical Review B 43 (1991), 


5364-5369. por: 10. 1103/PhysRev8 43-5384] (cit. on pp. #2} E3). 


R. D. Parmentier: “Fluxons in Long Josephson Junctions”. Solutions 
in Action. Ed. by K. Lonngren and A. Scott. Academic Press, 1978, 


pp. 173-199. por: 10. 1016/B978-0- 12-455580-8. 50014-7|(cit. on 
p- (i). 

N. F. Pedersen: “Solitons in Long Josephson Junctions”. Advances in 
Superconductivity. Ed. by B. Deaver and J. Ruvalds. NATO Advanced 
Science Institutes Series. Boston, MA: Springer US, 1983, pp. 149-181. 


DO1:|10.1007/978-1-4613-9954-4_5|(cit. on p.[ip. 


Bibliography 


[Ped86] 


[PS73] 


[Poo+14] 


[Poz11] 


[Pro+15] 


[qki] 


[Rey+50] 


[Rot+16] 


[Roy+21] 


[SN17] 


N. F. Pedersen: “Chapter 9 - Solitons in Josephson Transmission 
Lines”. Modern Problems in Condensed Matter Sciences. Ed. by S. E. 
Trullinger, V. E. Zakharov, and V. L. Pokrovsky. Vol. 17. Solitons. Else- 


vier, 1986, pp. 469-501. nor: 10. 1018/8978-0- 444-87002-8.50015-1] 
(cit. on p.|81). 

N. F. Pedersen and K. Saermark: Analytical solution for a Josephson- 
Junction model with capacitance. Physica 69 (1973), 572-578. Dot: [10 .] 


1016/0031-8914(73) 90089- X| (cit. on p.|37). 


C. Poole, H. Farach, R. Creswick, and R. Prozorov: Superconductivity. 


3rd ed. Elsevier Insights. Elsevier Science, 2014. DOI: 
(= O757S-(cit. on p. 6). 

D. M. Pozar: Microwave Engineering. 4th. Hoboken: Wiley, 2011. URL: 
https: //www.wiley.com/en-us/Microwave+Engineering+4tht+ 
[sdition-p- 9760470631663] (cit. on p.[19) 

S. Probst, F. B. Song, P. A. Bushev, A. V. Ustinov, and M. P. Weides: 
Efficient and robust analysis of complex scattering data under noise in mi- 
crowave resonators. Review of Scientific Instruments 86 (2015), 024706. 


por [f0.1068/1 4807935 (it. on pp. 22) 23). 


qkitgroup: Qkit — a quantum measurement and evaluation suite in python. 
uRL: {https ://github . com/qkitgroup/qkit|(cit. on p.|66). 

C. A. Reynolds, B. Serin, W. H. Wright, and L. B. Nesbitt: Superconduc- 
tivity of Isotopes of Mercury. Physical Review 78 (1950), 487-487. Dot: 


[0 1108/PhysRev.78.487|(cit. on p12). 
H. Rotzinger, S. T. Skacel, M. Pfirrmann, J. N. Voss, J. Miinzberg, S. 
Probst, P. A. Bushev, M. P. Weides, A. V. Ustinov, and J. E. Mooij: 
Aluminium-oxide wires for superconducting high kinetic inductance circuits. 
Superconductor Science and Technology 30 (2016), 025002. pox: [10.] 


1088/0953-2048/30/2/025002|(cit. on pp. {17} 55H57}. 
A. Roy, D. Schuricht, J. Hauschild, F. Pollmann, and H. Saleur: The 
quantum sine-Gordon model with quantum circuits. Nuclear Physics B 


968 (2021), 115445. por:[10. 1016/j .nuclphysb 2021. 115445] (cit. on 
pp- 2111). 


J. J. Sakurai and J. Napolitano: Modern quantum mechanics. 2nd ed. 


Cambridge University Press, 2017. DOI: [10 .1017/9781108499996] (cit. 
on pp. 102). 


127 


Bibliography 


[Sam+16] 


[Sar+15] 


[Sar+91] 


[SMU97] 


[Sch20] 


[Sch+20] 


[SB07] 


[Sch19] 


[Sch70] 


[Sco69] 


[Sco70] 


128 


N. Samkharadze, A. Bruno, P. Scarlino, G. Zheng, D. P. DiVincenzo, L. 
DiCarlo, and L. M. K. Vandersypen: High-Kinetic-Inductance Supercon- 
ducting Nanowire Resonators for Circuit QED in a Magnetic Field. Physical 


Review Applied 5 (2016), 044004. Do1:[10.1103/PhysRevApplied.5.] 
[644004|(cit. on p.[17). 

B. Sarabi, A. N. Ramanayaka, A. L. Burin, F. C. Wellstood, and K. D. 
Osborn: Cavity quantum electrodynamics using a near-resonance two-level 
system: Emergence of the Glauber state. Applied Physics Letters 106 


(2015), 172601. vor: |10.1063/1. 4918775] (cit. on p.P6). 


E. Sarnelli, S. Pagano, B. Ruggiero, and M. Russo: Magnetic field de- 
pendence of the critical current in long 'quasi-overlap’ Josephson junctions. 
IEEE Transactions on Magnetics 27 (1991), 2716-2719. DOI: 


[20133773] (it. on p.. 


V. V. Schmidt, P. Müller, and A. V. Ustinov: The physics of superconduc- 


tors. Berlin, Heidelberg: Springer, 1997. vor: 
[S5OF=8) it. on p. T6} 


A. Schneider: Quantum Sensing Experiments with Superconducting 


Qubits. 2020. vor (it. on pp. 
27) 28 B0 4 (102) 

Y. Schon, J. N. Voss, M. Wildermuth, A. Schneider, S. T. Skacel, M. P. 
Weides, J. H. Cole, H. Rotzinger, and A. V. Ustinov: Rabi oscillations 
in a superconducting nanowire circuit. npj Quantum Materials 5 (2020), 


1-5. por:|10. 1038/s41535-020-0220-x|(cit. on p.55}. 


J. R. Schrieffer and J. S. Brooks: Handbook of High-Temperature Supercon- 


ductivity. New York, NY: Springer, 2007. Dot: 
58734-6](cit. on p. IA} 


J. R. Schrieffer: Theory of Superconductivity. Boca Raton: CRC Press, 


2019. 352 pp. nor: [0 1201/8780429496709|(cit. on pp. 1213). 


K. Schwidtal: Type-I and Type-II Superconductivity in Wide Josephson Junc- 


tions. Physical Review B 2 (1970), 2526-2532. vot: 
2526). on pp. FA TB) 


A. C. Scott: A Nonlinear Klein-Gordon Equation. American Journal of 


Physics 37 (1969), 52-61. por:[i0. 111971. 1976404|(cit. on p.). 


W. C. Scott: Hysteresis in the DC switching characteristics of Josephson 
junctions. Applied Physics Letters 17 (1970), 166-169. pot: 


1653950] (cit. on p.87) 


Bibliography 


[SK05] 


[She+18] 


[SBM97] 


[Ste68] 


[Swe+13] 


[Swi61] 


[Tay+67] 


[Tie+21] 


[Tin04] 


P. D. Shaju and V. C. Kuriakose: Double-well potential for a Josephson 
vortex qubit. Physica C: Superconductivity and its Applications 424 


(2005), 125-132. por:[10.101673 -physe .2005. 05 .008](cit. on p.P). 

A. Shearrow, G. Koolstra, S. J. Whiteley, N. Earnest, P. S. Barry, F. J. 
Heremans, D. D. Awschalom, E. Shirokoff, and D. I. Schuster: Atomic 
layer deposition of titanium nitride for quantum circuits. Applied Physics 


Letters 113 (2018), 212601. por:|10 . 1063/1 . 5053461) (cit. on p.|17). 


A. Shnirman, E. Ben-Jacob, and B. Malomed: Tunneling and resonant 
tunneling of fluxons in a long Josephson junction. Physical Review B 56 


(1997), 14677-14685. por: 10.1103/PhysRevB 56. 14677|(cit. on pp.[52| 
97). 
W. C. Stewart: Current-voltage characteristics of josephson junctions. Ap- 


plied Physics Letters 12 (1968), 277-280. DOI: 
(cit. on pp. 32/33). 

L. J. Swenson, P. K. Day, B. H. Eom, H. G. Leduc, N. Llombart, C. M. 
McKenney, O. Noroozian, and J. Zmuidzinas: Operation of a titanium 
nitride superconducting microresonator detector in the nonlinear regime. 


Journal of Applied Physics 113 (2013), 104501. pot: 
[£754609) (cit. on pp. [71 23). 

J. C. Swihart: Field Solution for a Thin-Film Superconducting Strip Trans- 
mission Line. Journal of Applied Physics 32 (1961), 461-469. por:|10.| 


198371.1736028](it. on pp. [} B9). 
B. N. Taylor, W. H. Parker, D. N. Langenberg, and A. Denenstein: On 
the Use of the AC Josephson Effect to Maintain Standards of Electromotive 


Force. Metrologia 3 (1967), 89-98. DOT: |10 . 1088/0026- 1394/3/4/001 
(cit. on p. B0). 

E. Tiesinga, P. J. Mohr, D. B. Newell, and B. N. Taylor: CODATA 
recommended values of the fundamental physical constants: 2018. Reviews 


of Modern Physics 93 (2021), 025010. pox: 10. 1108/ReviiodPhys.98 | 
[025010] (cit. on p.|8). 

M. Tinkham: Introduction to superconductivity. 2nd ed. Dover books on 
physics. Mineola, New York: Dover Publications, 2004. URL: 


// store doverpublications. con/ 0486436032. html] (cit. on pp. [I 
Teje E 


129 


Bibliography 


[Ust98] 


[Ust05] 


[Val58] 


[Vis+10] 


[VD17] 


[Vos+21] 


[Wal+00a] 


[Wal76] 


[Wal+03a] 


[Wal+00b] 


130 


A. V. Ustinov: Solitons in Josephson junctions. Physica D: Nonlinear Phe- 
nomena. Annual International Conference of the Center for Nonlinear 


Studies 123 (1998), 315-329. vor: 
(cit. on pp.|1}50). 

A. V. Ustinov: Solitons in Josephson Junctions: Physics of Magnetic Fluxons 
in Superconducting Junctions and Arrays. John Wiley & Sons, 2005 (cit. 


on pp. [82] [83). 

J. G. Valatin: Comments on the theory of superconductivity. Il Nuovo 
Cimento (1955-1965) 7 (1958), 843-857. por: [10 . 1007/BF02745589| (cit. 
on p.|13). 

M. R. Vissers, J. Gao, D. S. Wisbey, D. A. Hite, C. C. Tsuei, A. D. 
Corcoles, M. Steffen, and D. P. Pappas: Low loss superconducting titanium 
nitride coplanar waveguide resonators. Applied Physics Letters 97 (2010), 


232509. por: 16. 1063/1. 3517269] cit. on p.[17). 


U. Vool and M. Devoret: Introduction to quantum electromagnetic circuits. 
International Journal of Circuit Theory and Applications 45 (2017), 


897-934. por:|10.1002/eta, 2359](cit. on pp. 2} 5). 

J. N. Voss, Y. Schön, M. Wildermuth, D. Dorer, J. H. Cole, H. Rotzinger, 
and A. V. Ustinov: Eliminating Quantum Phase Slips in Superconducting 
Nanowires. ACS Nano 15 (2021), 4108-4114. por:[10. 1021/acsnano .| 
[6.08721] (cit. on pp. 53,97. 

C. H. van der Wal, A. C. J. ter Haar, F. K. Wilhelm, R. N. Schouten, 
C. J. P. M. Harmans, T. P. Orlando, S. Lloyd, and J. E. Mooij: Quantum 
Superposition of Macroscopic Persistent-Current States. Science 290 (2000), 


773-777. vor: |f0.1126/science .280.5492.779 (cit. on p. 0}. 


J. R. Waldram: The Josephson effects in weakly coupled superconductors. 


Reports on Progress in Physics 39 (1976), 751-827. por:[10 .1088/0034-| 
(cit. on pp. ZT). 

A. Wallraff, T. Duty, A. Lukashenko, and A. V. Ustinov: Multiphoton 
Transitions between Energy Levels in a Current-Biased Josephson Tunnel 


Junction. Physical Review Letters 90 (2003), 037003. Dor: 
PhysievLett. 90 037003] cit. on p80). 


A. Wallraff, Y. Koval, M. Levitchev, M. V. Fistul, and A. V. Ustinov: 


Annular Long Josephson Junctions in a Magnetic Field: Engineering and 
Probing the Fluxon Interaction Potential. Journal of Low Temperature 


Physics 118 (2000), 543-553. por: [10 . 1023/A: 1004674908169] (cit. on 
p- B). 


Bibliography 


[Wal+03b] 


[Wal+04] 


[Wal01] 


[WM08] 


[Wei69] 


[Wil+22] 


[Wol11] 


[ZHB78] 


[ZB22] 


[ZTH70] 


A. Wallraff, A. Lukashenko, J. Lisenfeld, A. Kemp, M. V. Fistul, Y. 
Koval, and A. V. Ustinov: Quantum dynamics of a single vortex. Nature 


425 (2003), 155-158. por:[i0. 1038/nature01826|(cit. on pp. }84][109). 
A. Wallraff, D. I. Schuster, A. Blais, L. Frunzio, R.-S. Huang, J. Majer, 
S. Kumar, S. M. Girvin, and R. J. Schoelkopf: Strong coupling of a single 
photon to a superconducting qubit using circuit quantum electrodynamics. 


Nature 431 (2004), 162-167. por:|10. 1038/nature02851|(cit. on pp.{18] 
85). 

A. Wallraff: Fluxon dynamics in annular Josephson junctions: from rela- 
tivistic strings to quantum particles (2001) (cit. on pp. {53} (74). 


D. F. Walls and G. J. Milburn: Quantum optics. Ed. by D. F. Walls and 


G. J. Milburn. Berlin, Heidelberg: Springer, 2008. pox: [10 .1007/978-| 
B B40-2957Ea cit. on p 5). 


M. Weihnacht: Influence of Film Thickness on D. C. Josephson Current. 


physica status solidi (b) 32 (1969), K169-K172. por: |10. 1002/pssb.| 
[19690320359] cit. on pp. 40] 79). 

M. Wildermuth, L. Powalla, J. N. Voss, Y. Schön, A. Schneider, M. V. 
Fistul, H. Rotzinger, and A. V. Ustinov: Fluxons in high-impedance long 
Josephson junctions. Applied Physics Letters 120 (2022), 112601. por: 


10 . 1063/5 . 0082197] (cit. on pp. {41} [58} 59] |69} [72] [74] |78}180} 182] [84} 


109). 
E. L. Wolf: Principles of Electron Tunneling Spectroscopy: Second Edition. 


OUP Oxford, 2011. 618 pp. por 10, 1088/acprof 020/9780198E85154] 
[04.0603] (cit. on p. BT). 


P. Ziemann, G. Heim, and W. Buckel: Oxygen content and oxide barrier 
thickness in granular aluminum films. Solid State Communications 27 


(1978), 1131-1135. vor: (cit. on 
p. 56). 


J. Zierep and K. Bühler: Principles of Fluid Mechanics. 1st ed. Springer 


Wiesbaden, 2022. 237 pp. vor: |10. 1007/978-3-658-34812-0|(cit. on 
p. [). 
J. E. Zimmerman, P. Thiene, and J. T. Harding: Design and Operation 


of Stable rf-Biased Superconducting Point-Contact Quantum Devices, and 
a Note on the Properties of Perfectly Clean Metal Contacts. Journal of 


Applied Physics 41 (1970), 1572-1580. Dor: [10 . 1063/1. 1659074] (cit. 
on p. B0}. 


131 


List of Publications 


[Wil+22] M. Wildermuth, L. Powalla, J. N. Voss, Y. Schön, A. Schneider, M. V. 
Fistul, H. Rotzinger, and A. V. Ustinov: Fluxons in high-impedance long 
Josephson junctions. Applied Physics Letters 120 (2022), 112601. por: 
10.1063/5.0082197 


[Sch+20] Y. Schön, J. N. Voss, M. Wildermuth, A. Schneider, S. T. Skacel, M. P. 
Weides, J. H. Cole, H. Rotzinger, and A. V. Ustinov: Rabi oscillations 
in a superconducting nanowire circuit. npj Quantum Materials 5 (2020), 


1-5. D01:|10.1038/s41535-020-0220-x 


[Vos+21] J. N. Voss, Y. Schön, M. Wildermuth, D. Dorer, J. H. Cole, H. Rotzinger, 
and A. V. Ustinov: Eliminating Quantum Phase Slips in Superconducting 


Nanowires. ACS Nano 15 (2021), 4108-4114. vor:|10.1021/acsnano | 


0c08721 


133 


Appendix 


A Circuit quantization of long Josephson junctions 


A.1 sine-Gordon model of long Josephson junctions 
including high-kinetic inductance 


Long Josephson junctions can be modeled by N — œ infinitesimally small lumped 
elements of resistively and capacitively shunted junctions in z-direction, which 
are extended along the x-axis, and thus connected via inductive bottom and top 
electrodes (see Fig. [A.1a)). Using common circuit quantization techniques with N 
node fluxes in the top and bottom electrodes D; and oP, the Lagrangian holds 


N-1 
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=w | E (910° aa) zul (1 cos( g (o æ) )) (A.1b) 
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2L; 2L 


Here, cy = C}*/w£ denotes the specific junction capacitance, Ly = LpP'w// and 
Lg = Li?'w/l the top and bottom electrode’s sheet inductance that sums up both 
the geometric and kinetic parts and je = 1{'/w£ describe the critical and bias 
current density. 
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Figure A.1: Equivalent circuit of long Josephson junctions. a) A long junctions is modeled by lumped 
element junctions that are connected by inductive top and bottom electrodes, whose inductance may 
differ. b) In this circuit, the two inductance are equivalent to only one effective inductance in one 
electrode. 


Applying the Euler-Lagrange equation |BP82} |Leg87; |CF91] yields the coupled 


equations of motion 
d aL d aL o£ 
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Calculating (A.2b)/L5 — / Ly yields 
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(A.3a) 


Ind A 2 
= jdub — ke sin (20). (A.3b) 
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In the last step, the node fluxes 6'— 6’ = © are substituted by one single node flux, 
describing the difference of top and bottom electrode and the vector potentials A! — 
AP = A are treated analogously. Furthermore, an effective electrode inductance 


is identified as LPLP | > +% iP +10 =L 


b t 
electrode’s inductance. Hence, the circuit can be simplified to having only one total 


inductance as shown in Fig.|A.1b) 


į the sum of top and bottom 


to 
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A.2 Model of phase biased long Josephson junctions 


In order to analyze the phase distribution of a distributed long junction is divided 
into N — oo discrete segments and Kirchhoff’s current law is applied at each node 
j € [1, N]. In comparison to the high-kinetic inductance in the bottom electrode 
poes > L}", the top electrode’s purely geometric inductance L! is negligible, and 
the Ji node flux is coupled by the nearest neighbors as 


ey ee na 
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including the boundary conditions at the junction edges, defined by the bias 

current I. The node fluxes can be translated into node phases by means of the 

junction’s well-known flux-phase relation $j = $2 9. Besides, a homogeneous 

external magnetic flux causes a constant external phase drop between two adjacent 
xt 

nodes Bex, — pet = Pu = Agh y j. 
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By further identifying the Josephson inductance Ly = si and introducing the 
screening parameter BL = m = me, the conditional equations result in 
ext BL : BLY 
= 914 ; A. 
pt 
Pj+1 = 29) — Pj-1 NN (9), (A.6b) 
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B Experimental details 


B.1 Cryogenic measurement setup 
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Figure B.1: Spectra of RF and DC measurement lines. a) Spectrum of RF signal lines at room tem- 
perature. The attenuation of the stainless steel coaxial cables shows the typical 1/,/f behavior if the 
wavelength A is smaller than the diameter of the coaxial line ~ 0.3mm. b) Spectrum of DC measurement 
lines at cryogenic temperatures. The spectrum shows characteristics of two low-pass filters with ~ 30 dB 
attenuation each. The lower cutoff frequency © 1 kHz is attributed to an RCR filter at the second pulse 
tube stage, made from temperature insensitive thin-film resistors with R = 180 Q and multilayer ceramic 
(NPO) capacitors with C = 1nF. The upper cutoff frequency 300 kHz corresponds to the CLC low-pass 
filter at room temperature, consisting of shielded wire-wound SMD inductor with a ferrite core with 
L = 100 pH and multilayer ceramic (NPO) capacitors with C = 1nF. This type of LC circuit shows a 
different spectrum in resonance depending on whether the complex or real-valued signal is measured. 
The data in the low-frequency range from 1 Hz to 5 MHz are acquired with a Lock-in amplifier, whereas 
the microwave range from 14 kHz to 20 GHz is measured with a VNA. 
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B.2 Fabrication parameters 


Table B.1: Resist application. To coat a sample with a resist of thickness d, the spin coater accelerates 
with a, where the resist is distributed of the whole chip in a first ramping step with speed vr for the 
duration tr. In a second step, the targeted thickness is achieved by spinning with speed vs for the time 
ts. Afterwards the resist is softbaked at the temperature Tg for the time tg to vapor solvents and set the 
resist. 


resist d a Ur tr Us ts Tg tg 
(um) (rpm/s) (rpm) (s) (rpm) (s) (C) (s) 
S1805 ~0.5 6000 500 1 6000 60 115 60 
(optical) 
AZ5214E ~1.2 7500 500 5 6000 60 110 50 
(optical) 
ma-N 2403 ~0.3 1000 3000 65 90 60 
(e-beam) 
S1818 ~1.8 2000 4500 60 80 300 
(protective) 


Table B.2: Optical resist exposure parameters. A XeHg lamp emitting the wavelength A = 365nm 
is used to expose the photoresist under a mask for the time tex. Thereby a constant power mode, 
controlling the power Pex, or constant intensity mode, adjusting the dose Dex, is available. 


resist process mask Pex Dex fex 
(W) (mWem) (s) 

S1805 positive soda lime 370 1.8 
AZ5214E positive soda lime 13.0 5.0 


Table B.3: Electron-beam resist exposure parameters. Depending on the required resolution and 
accuracy of the desired structures, two different beam currents I, step sizes Ax are used to pattern the 
ma-N 2403 resist, whereas the average area dose D is kept constant. The small beam current and step 
size is used for writing the 200 nm wide long junction structures. The microstrip feedline and DC leads 
of test structures, however, are written with higher beam currents and step sizes in order to save writing 
time. 


resist D I Ax 
(Ccm?) (nA) (nm) 
ma-N 2403 240 2 10 


ma-N 2403 240 30 50 
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Table B.4: Resist development parameters. Each exposed resist is developed in a selective solvent 
(developer) for a time ty, which is stopped in a stop bath. 


resist developer tg stop bath 
(s) 
$1805 AZ developer : H20 (3:2) 33 H20 
AZ5214E AZ developer : H20 (1:1) 33 H,O 
ma-N 2403 MF-319 50 H20 


Table B.5: ICP trilayer etching parameters. Depending on requirements of (an-) isotropy and etching 
rate, different ICP processes are used that differ in participating gases, RF, and ICP power. The etching 
time is chosen to etch the desired layers as properly as possible. 


trilayer layers Tay Ico, Io, Pre Pıcr t 

(sccm) (sccm) (sccm) (W) (W) (s) 
A trilayer 15 15 0 50 100 205 
A top 15 15 0 50 100 90 
B trilayer 15 15 0 50 100 180 
B top 15 15 0 50 100 90 
C trilayer 15 15 0 50 100 215 
C top 15 15 0 50 100 110 
D trilayer 2 12 0 100 200 70 
D top 15 3 1 100 200 32 
E top 15 3 1 100 100 5 
E bottom 2 12 0 100 200 120 


Table B.6: Anodic oxidation parameters. The sample is facing a platinum anode in a fixed distance of 
2cm in an electrolyte of 156g ammonium pentaborate in 1120 ml ethylene glycol and 760 ml demineral- 
ized water. The oxide layer is defined by the applied voltage V and the measured current saturated 
over the time t. The influence of different voltages was studied at trilayer A, where three 5mm x 5mm 
in a row were dipped successive deeper and deeper in the solution by decreasing the voltage. 


trilayer chip V t 
(V) (min) 
A 1-3 9 15 
A 46 11 15 
A 7-9 13 15 
B 1-9 11 15 
C 1-9 11 15 
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Figure B.2: Static oxidation parameters for tunnel barrier creation. The critical current densities of 
trilayers A, B, and C are determined from short squared junctions measured at millikelvin temperatures. 
According to Kleinsasser et al. |JKMM95], je depends exponentially on the product of oxygen (partial) 


pressure Pox and oxidation time fox. 


C Additional measurement data 


Table C.1: Design parameters of Trilayer D for long junctions, each embedded in a resonator. A long 
junction of length £ and width w is placed in the center of a rectangular microstrip resonator of length 


l and width w. 


device l w £ w 
(um) (pm) (pm) (pm) 
Res #1 503.0 2.0 50.3 2.0 
Res #2 440.4 2.0 44.0 2.0 
Res #3 390.9 2.0 0.0 0.0 
Res #4 350.8 2.0 35.1 2.0 
Res #5 317.7 2.0 31.8 2.0 
Res #6 289.9 2.0 29.0 2.0 
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Figure C.1: Magnetic field dependence of coupled systems of a long junction embedded in a res- 
onator. The amplitude’s frequency spectrum |S21| of Res #1 in the lowest row to the Res #6 in the 
uppermost row is depicted as color depending on the applied magnetic field. The magnetic field 
dependence of the system’s resonant frequency (dark lines) differs significantly for the first (left), 
second (center), and third (right) harmonics caused by the LJJ’s magnetic field dependent inductance. 
(Note the different frequency scales for different harmonics in the three columns.) The eigenfrequencies 
of the reference resonator (Res #3) feature no magnetic field dependence. 
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In the last decades, superconducting devices have emerged as a promising 
platform for quantum technologies, including quantum sensing and quantum 
computing. Their key elements are Josephson junctions, which allow for coherent 
supercurrent tunneling between two weakly linked superconductors. If such a 
junction is extended in one direction to a long junction, the superconducting 
phase difference can vary in space and time and may allow for quantized phase 
windings that drive supercurrent vortices. The physics of such Josephson vorti- 
ces spans from nonlinear soliton dynamics with relativistic effects to technical 
applications like microwave generation and amplification. In allthese cases the 
vortices behave as non-quantum particles. This restriction to classical behavior 
originates in the long junction’s limited geometrical properties, in particular its 
electrode inductance. 

The advent of superconducting high-kinetic inductance circuits relaxes these 
constraints and enables an increased junction impedance, which goes along 
with the vortices’ quantumness. In this work it is demonstrated that the junction 
impedance can be enhanced beyond the geometric limitations, which facilitates 
various applications. A hybrid system is studied particularly, that consists of a 
high-impedance long junction embedded in a microwave resonator. This galvanic 
coupling enables the dispersive readout scheme to determine the quantum states 
of different vortex configurations. In the vortices’ quantum regime, coherent 
tunneling of single vortices in a two-level system is observed, forming a Joseph- 
son vortex quantum bit. Quantum Rabi oscillations with energy relaxation and 
dephasing times in the microsecond range are measured, making the system 
promising for future quantum technologies. 
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